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This review discusses how low-energy, valence excitations created by swift electrons can render 
information on the optical response of structured materials with unmatched spatial resolution. 
Electron microscopes are capable of focusing electron beams on sub-nanometer spots and probing 
the target response either by analyzing electron energy losses or by detecting emitted radia- 
tion. Theoretical frameworks suited to calculate the probability of energy loss and light emission 
(cathodoluminescence) are revisited and compared with experimental results. More precisely, a 
quantum-mechanical description of the interaction between the electrons and the sample is dis- 
cussed, followed by a powerful classical dielectric approach that can be in practice applied to 
more complex systems. We assess the conditions under which classical and quantum-mechanical 
formulations are equivalent. The excitation of collective modes such as plasmons is studied in 
bulk materials, planar surfaces, and nanopartielcs. Light emission induced by the electrons is 
shown to constitute an excellent probe of plasmons, combining sub-nanometer resolution in the 
position of the electron beam with nanometer resolution in the emitted wavelength. Both electron 
energy-loss and cathodoluminescence spectroscopics performed in a scanning mode of operation 
yield snap shots of plasmon modes in nanostructures with fine spatial detail as compared to other 
existing imaging techniques, thus providing an ideal tool for nanophotonics studies. 
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I. INTRODUCTION 

Investigating minute physical processes has been es- 
sential for advancing knowledge and generating the right 
questions all the way from the beginning of modern sci- 
ence up to recent developments in nanotechnology. Far- 
field optical microscopes have contributed to this end, 
but they are limited by diffraction to a spatial resolution 
of about half the light wavelength in practice. Actually, 
the best resolution (< lA) is currently achieved using 
electrons rather than light (Nellist et ai, 2004; Ruska, 
1987), for which the de Broglie wavelength (de Broglie, 
1925) is well below 0.1 A at energies above 20 keV. 

Besides acquiring static images of the nanoworld, we 
are interested in finding out about the excitations that 
small objects can sustain, which inform us on their dy- 
namical evolution and are relevant for encoding and ma- 
nipulating information and for exploring a myriad of ap- 
plications in fields such as molecular biology. Thus, our 
ultimate goal is to perform spectroscopy at the shortest 
possible length scale with the highest energy resolution. 
With this focus in mind, we can classify the available 
experimental techniques as shown in Fig. 1, clearly indi- 
cating that electron-based spectroscopies offer the best 
choice for optimizing spatial resolution. 

Two routes have been devised so far for dealing with 
the diffraction limit :^ (1) to reduce the wavelength of 
the probe and (2) to employ near-field detection tech- 
niques. (1) Moving to shorter photon wavelengths is 
difficult due to the lack of versatile lenses and mirrors 
beyond the ultraviolet. Nonetheless, recent advances in 
x-ray microscopy (XRM) have allowed 15nm imaging 
resolution (Chao et ai, 2005) at energies above 250 eV. 
Alternatively, shifting from photons to electrons of the 
same energy encompasses a significant decrease in wave- 
length and improved spatial resolution.^ (2) Exploiting 
the near field is another option, particularly when rely- 
ing on evanescent components, the fast decay of which 
provides an extra handle for enhancing resolution using 
localized probes. For instance, near-field scanning optical 
microscopy (NSOM, based on a subwavelength tip at the 



^ The existence of a diffraction Umit has been questioned by recent 
investigations of super-oscillating functions (Zheludev, 2008). 

^ Heavier particles such as protons have been argued to pro- 
vide good spatial resolution compared to electrons (Demkov and 
Meyer, 2004). Actually, a helium ion microscope has been re- 
cently released (Ward et ai, 2006) and achieved 2.4 A resolution 
(Carl Zeiss 2008 announcement). 
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FIG. 1 (Color online) Atlas of spatially-resolved spectroscopy 
techniques, organized according to their space and energy 
resolution (see Appendix D for a list of acronyms). The 
relation between wavelength and energy is represented by 
dashed curves for photons and electrons. The universal elec- 
tron inelastic-mean-free-path (e-imfp) is given as a function 
of electron energy. The skin depth of Ag is calculated from 
optical data (Palik, 1985). 



end of a fiber that brings, collects, or scatters light) and 
tip-enhanced NSOM can push spatial resolution down 
to tens of nanometers (Betzig et al., 1992; Hartschuh, 
2008). Even better detail in the sub-A domain is rou- 
tinely achieved using the formerly developed scanning 
tunneling microscope (STM), on which scanning tunnel- 
ing spectroscopy (STS) yields 0.1 eV energy resolution 
(Hormandinger, 1994). However, STS only works with 
metals, and its extreme spatial accuracy (~ 0.01 A in the 
vertical direction) requires close proximity between the 
tip and the sample surface, so that single-electron exci- 
tations dominate the spectra and mask collective modes 
relevant to optics, such as plasmons. 

Electron microscopes are thus the best option for re- 
solving both localized and extended excitations with sub- 
nanometer spatial detail and < 0.1 eV energy resolution 
in any type of material (Lazar et al., 2006). These in- 
struments are sufficiently versatile to be surface-sensitive 
and to simultaneously procure information on bulk prop- 
erties. Their performance has considerably improved in 
recent years due to extraordinary advances in energetic 
electron sources and optics. A number of studies in dif- 
ferent fields benefit from the unprecedented spatial reso- 
lution of electron microscopy, which no other technique 
can currently match. 

Electron microscopy has different ffavors depending 
on the type of signal that we measure [inelastic losses, 
cathodoluminescence (CL) light emission, secondary elec- 
trons (SEs), etc.]. Several widely available types of se- 
tups are schematically represented in Fig. 2. The trans- 
mission electron microscope (TEM) provides by far the 
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most powerful combination of spectral and energy reso- 
lution, mainly via analysis of loss events. This type of 
microscope can operate like a conventional optical mi- 
croscope, projecting a bright image of an object on a 
phosphor screen or a CCD camera to picture the mag- 
nified sample. Energy-filtered TEM (EFTEM) allows 
constructing images out of those electrons that have do- 
nated a determined amount of energy to the specimen, 
which is useful to visualize selected losses [for instance, 
core excitations identifying the chemical environment of 
atomic species (Batson, 1993)]. Excellent spatial resolu- 
tion is obtained using TEMs operated in scanning mode 
(STEM) , in which the beam is focused and scanned over 
the sampled area, either to form images out of large- 
momentum transfers recorded by annular dark-field elec- 
tron detectors or to collect loss spectra at specific loca- 
tions. 

Most TEMs use energetic electrons in the range of 50- 
400 keV and require very thin specimens 100 nm) 
that are relatively transparent to these electrons. Sam- 
ple preparation is thus an important issue for achieving 
optimum spatial resolution. 

In contrast, scanning electron microscopes (SEMs) can 
work with any sample that is covered by a thin metal 
layer (~ 1 — 2 nm), forming images by collecting SEs 
upon bombardment with a primary beam of ^ 2 nm in 
diameter. Furthermore, one can perform spectroscopy 
on SEMs through wavelength-resolved CL emission. Al- 
though CL spectroscopy can be also collected on a TEM 
(Yamamoto et ai, 1996a), it is most commonly available 
in SEMs equipped with a light collector system (e.g., an 
ellipsoidal mirror) and a photon spectrometer. 

Low-energy electron microscopy (LEEM) (Rocca, 
1995) relies on the use of 1-100 eV electrons and can 
achieve < 50 nm spatial resolution accompanied by sub- 
electronvolt energy resolution. As an example of appli- 
cation, LEEM has been recently used to observe acoustic 
plasmons (Diaconescu et ai, 2007). In contrast, pho- 
toemission electron microscopy (PEEM) (Bauer, 1994) 
is particularly useful to map the light intensity in illumi- 
nated nanostructures, as recently demonstrated for Utho- 
graphically patterned Ag nanoparticle arrays (Aeschli- 
mann et al., 2007; Cinchetti et ai, 2005). The femtosec- 
ond dynamics of plasmon excitations has been also re- 
solved using multiphoton PEEM with pump-probe illu- 
mination (Kubo et ai, 2005, 2007). 

Finally, there are other configurations that make use 
of electrons to perform spectroscopy. For instance, re- 
flection electron energy-loss spectroscopy (REELS) has 
been recently used to determine optical properties of no- 
ble metals (Went et ai, 2008; Werner, 2006; Werner et ai, 
2007) and graphite (Calliari et ai, 2008) after careful 
data analysis. 



(a) STEM 

80-400 keV 
electron 
beam 



(b) SEM 

1-50 keV 
electron 
beam 




thin 
sample 



ADF 



SEE 



(c) PEEM 

electron 
microscope 

<50-eV 




sample 



scanning imaging 
& EELS 



1 



V, If Pl^otoelectrons 



50 nm 

sample 



FIG. 2 (Color online) Schematic representation of various 
types of electron microscopes equipped to perform spec- 
troscopy, (a) Scanning transmission electron microscopes 
(STEMs) allows analyzing the energy distribution of elec- 
trons transmitted through a thin specimen (< 100 nm thick), 
on which the incident beam is focused down to a < 0.1 nm 
spot. Large-momentum transfers are collected through annu- 
lar dark-field detectors (ADFs). Cathodoluminescence (CL) 
light detection is also possible on a STEM, (b) Scanning elec- 
tron microscopes (SEMs) collect secondary electrons (SEs) to 
form images or CL emission to perform spectroscopy, with 
spatial resolution down to ~ 2nm. (c) Photoemission elec- 
tron microscopes (PEEMs) can render spectroscopic informa- 
tion on the photo-excitation process and electron-hole (e-h) 
dynamics with < 50 nm resolution. 



A. Spectroscopy using electron microscopes 



A swift electron impinging on a microscope specimen 
gives rise to secondary electron emission (SEE) and CL 
that can be energy-analyzed, thus yielding information 
on the excitation spectrum of the sample. Electron 
energy-loss spectroscopy (EELS) performed in STEMs 
adds up to this suite of probes. The main advantage 
of electron microscopes is that these types of spectro- 
scopies can be performed with truly nanometer spatial 
resolution. 

Unlike freely propagating light, the bare field of an 
electron is evanescent, as we show in Sec. II. A. This 
is advantageous to investigate localized excitations, in- 
volving wave-vector components of the electromagnetic 
field that lie outside the light cone. For example, elec- 
tron beams were instrumental in discovering and char- 
acterizing collective excitations of conduction electrons 
in metals (Chen and Silcox, 1975b; Powell and Swan, 
1959; Raether, 1980; Ritchie, 1957; Ruthermann, 1948; 
Watanabe, 1956), known as plasmons because they are 
quasi-particles associated to oscillations of the electron 
gas plasma. Specifically, bulk plasmons in a homogenous 
metal are conspicuous in EELS, since they are electro- 
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static and longitudinal in nature,^ and for this reason 
they couple very efficiently to moving charges. 

Similarly, plasmons can be confined at the surface of 
a metal, and they are actually the source of interesting 
phenomena and applications that configure the field of 
plasmonics (see Sec. VI). Surface plasmons (SPs) are ver- 
satile entities: they can be trapped in metal particles that 
are much smaller than the wavelength (Myroshnychenko 
et al., 2008b), and they can hybridize with light extend- 
ing over larger metallic structures (Coyle et ai, 2001). 
More precisely, surface plasmon polaritons (SPPs), which 
are a subset of SPs capable of propagating on planar sur- 
faces (Barnes et ai, 2003; Ozbay, 2006; Zia et al., 2006) 
or along ID waveguides (Bozhevolnyi et al., 2006), hold 
great promise to become the natural link between current 
nanoelectronics and future integrated nanophotonics, op- 
erating at frequencies that are ^ 10^ times higher than 
microchip clocks. 

As illustrated in Fig. 3, the passage of a fast electron 
can excite localized plasmons in metallic nanoparticles, 
but it can also launch SPPs in planar metallic surfaces 
(Bashevoy et ai, 2006; Heitmann, 1977; van Wijngaarden 
et ai, 2006) or in metal nanowires (Vesseur et ai, 2007). 
Quite different from SPPs, localized plasmons can decay 
radiatively, thus contributing to CL, although all sorts of 
plasmons can partially decay via inelastic channels that 
involve electronic excitations, including electron-hole (e- 
h) pair creation and SEE if the electron is near the surface 
and its energy above the vacuum level. We discuss this 
matter more thoroughly in Sec. IV. A. 

The electron can directly excite e-h pairs too, the decay 
of which gives rise to the emission of CL and SEs. This 
produces in general a complex cascade of SEs, which a 
transport equation approach is suited to model reliably 
(Rosier and Brauer, 1991). The cascade includes ener- 
getic electrons that generate further SEs and CL, adding 
complexity to this scenario. However an element of sim- 
plicity comes from the separation between coherent and 
incoherent processes, from the point of view of the emit- 
ted light, as we explain in more detail in Sec. IV. In 
particular, the coherent CL signal is dominant in metals. 

We are interested in low-energy excitations, which in- 
volve electrons and holes in the valence and conduction 
bands, as well as collective modes (e.g., localized plas- 
mons and SPPs). This is by far the most intense part of 
the loss spectrum (see Fig. 4). Its analysis yields infor- 
mation on the materials optical response with the kind 
of spatial resolution that is currently dreamed of in the 



Bulk plasmons in a source-free metal are characterized by zero 
magnetic field and longitudinal electric field (V X E = 0), so 
that they trivially satisfy Maxwell's equations under the condi- 
tion of vanishing permittivity. In contrast, surface plasmons are 
confined to metal-dielectric interfaces, they involve nonzero mag- 
netic fields, and they have transversal character (V ■ E = in 
each homogeneous region of space separated by the interface on 
which plasmons are defined). 




FIG. 3 (Color online) Schematic representation of some of 
the excitation processes triggered in a solid by a swift passing 
electron, and their connection to measured signals in elec- 
tron microscopes: CL, secondary electron emission (SEE), 
and electron energy-loss spectroscopy (EELS). The electron 
can launch propagating surface plasmon polaritons (SPPs; see 
Sec. V.D). It can also produce surface plasmons localized in 
nanoparticles or in surface features (Sec. VI). Localized sur- 
face plasmons typically decay by coupling to radiation, thus 
giving rise to contribution (1) to CL (Sec. IV. A. 1). Elec- 
tronic excitations in the sample (e.g., e-h pairs) can decay 
radiatively, back to the initial state [coherent CL emission 
(2), automatically included in the random-phase approxima- 
tion (RPA) dielectric function as a bubble diagram (Lindhard, 
1954); see Sec. IV.A.l], or to a different excited state (incoher- 
ent CL emission; Sec. IV. A. 2). In metals, secondary electronic 
excitations constitute the dominant decay channel, producing 
Auger electrons (AEs), sometimes above the vacuum level, so 
that they contribute to the detected SEE. Direct excitation 
of electrons from the target is also possible to produce true 
SEs. 

context of nanophotonics. It is our purpose to review 
historical developments, to summarize recent advances in 
this area, to present an overview of theoretical methods, 
and to point out some opportunities opened by electron 
microscopy in order to expand and complement nanopho- 
tonics studies in a way that can be particularly beneficial 
for emerging areas such as plasmonics. 



II. INTERACTION OF SWIFT ELECTRONS WITH 
MATTER 

We first consider some fundamental aspects of the in- 
teraction between fast electrons a matter, as well as the- 
oretical approaches suitable to simulate EELS. 



A. An evanescent source of light In matter 

The electromagnetic field that accompanies a point 
charge moving in vacuum can be regarded as an evanes- 
cent source of radiation which permits exploring regions 
of momentum-energy space that lie outside the light cone. 
This has interesting consequences: fast electrons gener- 
ate SPPs when passing near a metal surface, as we show 
in Sec. III.C, but they can also excite non-dipolar modes 
in small particles (Chu et ai, 2009), which would be dif- 
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FIG. 4 (Color online) Typical electron energy-loss spectrum 
showing a narrow zero-loss peak (ZLP; FWHM~0.2eV), col- 
lective modes in the valence loss region (e.g., the 3.7 eV plas- 
mon of Ag in the inset), and much weaker core excitations 
at higher lost energies. The spectrum has been calculated 
for 200 keV electrons traversing a 20 nm Ag film using optical 
data for Ag (Palik, 1985) and assuming a collection angle of 
5 mrad. 



ficult to resolve using external light instead. 

It is helpful to examine first the electric field produced 
by a swift charged particle moving inside a homogeneous 
medium. The field can be conveniently decomposed into 
different frequency contributions using the Fourier trans- 
form 



E(r,t) 



2^ 



E(r,w)e" 



(1) 



and also into momentum components with exp(iq-r) spa- 
tial dependence. We assume that the electron describes 
a straight-line trajectory with constant velocity vector v, 
crossing the origin at time i = 0. This is a reasonable as- 
sumption for the swift electrons and relatively thin speci- 
mens typically examined with TEMs. Then, the electron 
charge density becomes —^-KeSiio — q • v) in q — space. 
Direct solution of Maxwell's equations yields 



E(r,c^) = 



q/e — fcv/c 



e"i"-5(c.-q.v), (2) 



where k — lu/c is the light wavenumber in free space and 
e is the dielectric function of the homogeneous medium. 
It should be noted that e can depend on both q and 
uj. The consequences of the wavevector dependence are 
discussed later in Sec. III.B, but they are ignored in the 
remainder of this section. 

Equation (2) contains some interesting elements. Re- 
tardation effects show up both through k^e in the de- 
nominator, which modifies the strength of the field 
(producing, for example, weaker interaction in metals 
and stronger interaction in low-index dielectrics), and 
through a term proportional to the velocity vector in 
the numerator, which adds transversal components to the 
field. 
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FIG. 5 (Color online) Wavevector-frequency diagram showing 
the light dispersion relation — k^e in a polaritonic material 
described by Eq. (A3) (thick solid curves) and its intersection 
with allowed transfers coming from an electron moving with 
velocity v (shaded region). Point C signals the frequency 
threshold of Cherenkov radiation (CR) emission. 



The delta function inside the integral of Eq. (2) ex- 
presses energy conservation for transfers of frequency uj 
and wavevector q from the electron to the material. Ne- 
glecting relativistic corrections, the electron energy is re- 
duced from meW^/2 to |r7ieV — 7iqp/2me during the trans- 
fer, so that the energy difference Ticu leads to 



q • V — hq^/2me ~ q • v. 



(3) 



The right-most expression is the non-recoil approxima- 
tion, which works extremely well under the usual condi- 
tion hq ^ m^v.'^ This approximation remains valid if 
the energy transfer is computed using relativistic expres- 
sions. The kinematically allowed transfers span a solid 
area for all possible relative orientations between v and 
q, as shown by the shaded region in Fig. 5 with upper 
boundary to — qv. 

The zeros of the denominator in the integrand of Eq. 
(2) signal the dispersion relation of light in the medium, 
q = kyft. This has been represented in Fig. 5 for a polari- 
tonic material described by the dielectric function of Eq. 
(A3) (thick solid curves). The figure shows two different 
frequency domains separated by the condition v^t — c? 
(point C in Fig. 5 is a graphical solution of this equation). 
At lower frequencies with respect to C, the electron does 
not couple to excitations in the medium and the spec- 
tral components of the electric field decay exponentially 
away from the trajectory, as explicitly shown by solving 



* Electrons with typical TEM energies above 80 keV have veloc- 
ities V > 0.5c. When they undergo valence losses Tiuj < 50 eV, 
the momentum transfer is Tiq-^ ~ hio/v < 0.0004 met) along the 
direction of the trajectory. 
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FIG. 6 (Color online) Transverse spatial-direction depen- 
dence of the exp(— itjt) contribution to the electromagnetic 
field set up by an electron moving in vacuum with velocity 
V — 0.7c (7 « 1.4 and kinetic energy « 200 keV) along the 
positive z axis. The only non- vanishing components {Eu, Ez, 
and H^) decay exponentially at large distance R from the tra- 
jectory. The inset shows the orientation of these components 
relative to the electron velocity vector. The small R limit is 
dominated by the 1/R divergence of Er and H^p. 



electromagnetic field extends up to distances of the or- 
der of ~ t'7e/w, the Bohr cutoff. Notice however that one 
cannot assign this value to a characteristic Coulomb de- 
localization distance, since the field diverges at the origin 
as ~ so that large interaction contrast is expected 
across small distances in the region close to the trajec- 
tory. 

When the electron moves faster than light in the 
medium, under the condition 

V > c/yfe, (7) 

there is overlap between the photon dispersion relation 
and the shaded region in Fig. 5, so that the electron can 
emit Cherenkov radiation (CR) (see Sec. IV. D). The field 
displays oscillatory behavior and decays as 1/ VR away 
from the trajectory [this stems from the modified Bessel 
functions Kg and Ki for imaginary argument in Eq. (5)]. 

It should be stressed that Eqs. (4)- (6) and Fig. 6 re- 
fer to each monochromatic component of the electromag- 
netic field, evolving with time as exp{—iujt) [see Eq. (1)]. 
An electron moving in vacuum can be actually regarded 
as an evanescent source of super-continuum light with 
the spectral dependence shown in Fig. 6. 



the integral in Eq. (2). We find (Jackson, 1999) 



where 

g(r) = e''^^/" 



1a-„ ^ Z-A-, ii^lR 

7€ \VleJ \V^. 



(4) 



(5) 



7e — 1/ y/l — ev'^/c^ is the Lorentz contraction factor, 
and the notation r — (R, z) with R = {x,y) has been 
employed (see inset of Fig. 6), with the velocity taken 
along z. Similarly, the magnetic field reduces to 



H(r,c.) 



uc7j \vjf^ 



(6) 



where tp is the azimuthal unit vector (see Fig. 6). The 
announced exponential decay with R of both E(r,w) 
and H(r,a;) arises from the asymptotic behavior of 
the modified Bessel functions Km for large arguments 
(Abramowitz and Stegun, 1972). The non- vanishing 
components of these fields are represented in Fig. 6. The 



^ Interestingly, the Lorentz factor 7^ involves the velocity of light 
in the material, cj^fi, which is in turn frequency-dependent. 

® The square roots are chosen to yield positive real parts in this 
work. Notice that Im{e} is always positive in the retarded re- 
sponse formalism followed here, and it becomes a positive in- 
finitesimal in non- lossy dielectrics. 



B. Classical dielectric formalism 

The pioneering work of Fermi (1940) on the stopping 
of fast charged particles in dielectric materials opened up 
the application of classical electrodynamics to describe 
the interaction of swift electrons with matter. Following 
this useful tradition, we now extend the dielectric formal- 
ism of the previous section to cope with inhomogeneous 
samples and discuss in particular the spectral loss prob- 
ability, which is relevant to EELS experiments. 

The energy loss suffered by a fast electron moving with 
constant velocity v along a straight line trajectory r = 
re(i) can be related to the force exerted by the induced 
electric field E'"*^ acting back on the electron as (Ritchie, 
1957) 

AE = ej divE'"'^[re(t),i] 
hLuduj Feels (w). 



where the — e electron charge has been included (i.e., 
AE > 0) and 

TeelsH = ^ JdtRc {e-'"*v • E^^^ [re(t),c^]} (8) 

is the so-called loss probability, which is given per unit of 
transferred frequency uj. The problem of calculating the 
loss probability reduces then to solving the electric field 
set up by the electron. A large deal of work has been de- 
voted to obtaining the electric field for many geometries, 
including planar surfaces, isolated spheres, neighboring 
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spheres, circular cylinders, wedges, and more complex 
shapes, using both analytical and fully-numerical meth- 
ods, either within the non-retarded approximation, based 
upon solutions of Poisson's equation, or with full inclu- 
sion of retardation effects by solving Maxwell's equations 
(see Sec. Ill and referenc;es therein). Next, we outline the 
general features of this formalism. 



in terms of a frequency- and space-dependent dielectric 
function e(r,a;) provides often reliable results. Detailed 
expressions of W are offered in Appendix B for planar, 
spherical, and cylindrical geometries. Nevertheless, Eq. 
(10) is valid beyond the local response approximation, 
and subtle effects such as local-field corrections can be 
incorporated via more realistic quantum-mechanical rep- 
resentations of the screened interaction. 



1. Non-retarded approximation 



In the non-retarded approximation, we neglect the de- 
lay experienced by the electromagnetic signal that medi- 
ates the electron-sample interaction. Then, the electric 
field admits the form E(r, w) = —V<j){r,u)), and we can 
disregard H in the absence of magnetic response. 

It is useful to express the electric potential (j) in terms 
of the screened interaction W{r, r' , to), defined as the po- 
tential created at r by a unit point charge located at r' [an 
implicit exp(— iojt) time dependence is understood]. This 
quantity has to be combined with the charge density cor- 
responding to the moving electron. Considering a non- 
recoiled straight line trajectory, and assuming without 
loss of generahty that the velocity vector is directed along 
the positive z axis [i.e., the trajectory is r,,{t) = tq + vt 
and V = vz], the electron charge density in frequency 
space CO reduces to 



e j dte"^' (5(r-ro- vi) 
'^(R - Ro) e^'^(^-^°'/^ 



(9) 



where the notation shown in the inset of Fig. 6 has been 
used. From here, the potential reads 



(r,u;) 



dz' VF(r,Ro,z',a;) e' 



\w{z' — zo)lv 



Finally, plugging these expressions into Eq. (8), the non- 
retarded (NR) loss probabihty is found to be 



rgfLs(Ro,a;) = ^ j dzdz'c 
X Im{-VF(Ro,2:,Ro,-2',w)} 



i0{z — z') 



(10) 



where the dependence of r^gm on the electron impact 
parameter Rq is explicitly shown. 

The loss probability can be thus derived from the 
knowledge of the frequency-dependent screened interac- 
tion W. A local description of the sampled materials 



It should be noted that the reciprocity theorem \W{r,r',LS) = 
VF(r',r, a;)] and the fact that the bare Coulomb interaction is a 
real function have been utilized in the derivation of Eq. (10). 
® The induced field in Eq. (8) can be safely replaced by the total 
field, because the bare field of the moving charge does not pro- 
duce stopping. We have aecordingly dropped the superscript ind 
in Eq. (10). 



2. Retardation effects 

In high-voltage TEMs, retardation may become im- 
portant because the speed of the charged projectiles is 
a sizeable fraction of the speed of light. For instance, 
the Lorentz contraction factor at 200 keV takes already 
a value significantly different from one, 7 = 1.4. Retar- 
dation has two noticeable consequences for spectroscopy: 
(1) it increases the range of interaction of the electron 
probe in directions normal to the trajectory, as clearly 
illustrated by Eq. (5) and Fig. 6; and (2) it produces 
redshifts in excitation mode energies (Myroshnychenko 
et ai, 2008b). The latter are discussed below in further 
detail, but we can anticipate than this effect turns out to 
be important when the excitations extend over specimen 
distances that cannot be neglected in front of the corre- 
sponding light wavelength. A proper description of these 
effects requires calculating the electric field of Eq. (8) 
from Maxwell's equations in the presence of the moving 
electron and the sample under consideration. 

The electromagnetic response of a structured material 
is fully captured in its electric Green tensor. In partic- 
ular, the electric field produced by an external current 
density j{r,Lo) in an inhomogeneous medium of permit- 
tivity e(r, w) can be written in frequency space oj as 

E(r,a;) = -Amu) J rfr'G(r,r',w) ^(r',^) (11) 

in terms of G, the electric Green tensor of Maxwell's 
equations in Gaussian units, satisfying 



V X V X G(r, r', w) - k'^ e(r, lo) G(r, r', w) 
= ^5(r-r') 



(12) 



and vanishing far away from the sources (|r — r'| — > 00 
limit) . 

For the electron charge density of Eq. (9), the external 
current density reduces to j = vp, which upon insertion 
into Eq. (11), and this in turn into Eq. (8), allows us to 
write the loss probability as 



4e^ r 

Feels (Ro,w) = / dzdz' cos 
X Im{-Gzz[Ro, z,Iio, z' ,oj]} , 



co{z — z') 



(13) 
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where Gzz = z ■ G ■ z.^ Interestingly, this expression 
works for any sign of u, and therefore the loss probability 
is independent on whether the electron moves towards 
positive or negative z's. 



C. Quantum approach 

The quantum nature of both the electron probe and the 
excitations sustained by the targeted materials permeate 
many aspects of the electron-sample interaction. How- 
ever, Ritchie and Howie (1988) showed that a quantum 
mechanical description of EELS yields the same results 
as a semiclassical formalism if all the inelastic signal is 
collected, and therefore readers that are more interested 
in applications and in a general overview of the subject 
can overlook this section, unless they feel motivated by 
Sec. V.B.2. 

We work for simplicity in the non-retarded limit, in 
which the validity of Eq. (10) is demonstrated next un- 
der very common experimental conditions (Ritchie and 
Howie, 1988). 



1. Quantum description of the target 

The last century has witnessed important develop- 
ments in the field of interaction of fast charges with solids, 
particularly in the community of electronic and atomic 
collisions (see, for instance, Palmer and Rous, 1992, 
Ziegler, 1999, Winter, 2002, and references therein), 
from which our theoretical understanding of electron mi- 
croscopy has benefited considerably. In that context, ion 
and electron stopping was of great importance to un- 
derstand dynamical screening in solids, with interesting 
developments such as the description of plasmon excita- 
tions using second quantization schemes in planar (Lu- 
cas and Sunjic, 1971) and spherical surfaces (Ashley and 
Ferrell, 1976). Moreover, a non-retarded, fully quantum- 
mechanical expression for the loss probability was de- 
rived using a self-energy formalism (Echenique et al, 
1987, 1990). A quantum treatment of the target has been 
shown to be necessary for a correct assessment of delocal- 
ization in the excitation of core levels (Allen et al., 2003; 
Oxley and Allen, 1998). We here obtain a general ex- 
pression for the loss probability, starting from the more 
widely used Fermi golden rule, and assess the conditions 
under which it agrees with the semi-classical formalism 



In the retarded case, the reciprocity theorem states that 
G{r, r',aj) = G^(r',r,a;), a fact that we have used to recast 
exponential factors involving z and z' into a cosine function in 
the derivation of Eq. (13). Moreover, the free-space Green func- 
tion is entirely made of plane-wave components lying inside the 
light cone, so that it cannot contain wavevectors ui/v > k. This 
guarantees that the integral in Eq. (13) yields zero for an elec- 
tron moving in vacuum, and that we are allowed to utilize the 
total rather than the induced field to obtain this equation. 



final state 




initial state 



FIG. 7 (Color online) Schematic representation of the 
Coulomb interaction between a swift electron and a speci- 
men, showing the elements involved in Eq. (14), including 
the target electron-density operator p at point r'. 



presented above, as previously discussed by Ritchie and 
Howie (1988). 

During the interaction of a fast electron with a tar- 
get, the latter can undergo transitions from its ground 
state |0) of energy Uujo to excited states \n) of energy 
?iujm while the incoming electron of energy hsi and wave 
function '0i(r) acquires components ■(/'/(r) of lower en- 
ergy TiEf. Since the interaction with very energetic elec- 
trons is generally small, the transition rate is well de- 
scribed within first-order perturbation theory (Fermi's 
golden rule): 



rfF^R 
dt 



2^e^ 



E 



dhd^r 



(14) 



where we have used the target-probe Coulomb interac- 
tion and p is the target electron-density operator (see 
Fig. 7). We can now recast this expression into a more 
convenient form by relating the target matrix elements 
to the linear-response susceptibility (Pines and Nozieres, 
1966) 

Im{xir,r',uj)} 

= X E(0|/5(r)l") {n\p{r'm 5{^n -coo- cj), 

n 

valid for lu > 0, and this in turn to the screened interac- 
tion 



r - ri||r' - 



Then, the sum of Eq. (14) can be separated into specific 
values of the frequency transfer uj = Si — e j as 



duj 



dt 



dt ' 
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where 



dt 



J2 / d'rd^r' V/(r)V*(r)VKr')V'^(r') 



X lm{-W{r,r',w)}S{sf-Si+u}). (15) 



This expression is general for incident electrons of well- 
defined energy, and we remind the reader that W con- 
tains all quantum-mechanical details of the sample re- 
sponse, although expressions obtained from dielectric 
theory such as those offered in Appendix B yield reliable 
results in most situations encountered in practice. 



2. Quantum effects in the fast electrons 

Several authors have analyzed the formation of im- 
ages in electron microscopes from a quantum-mechanical 
viewpoint, considering the influence of instrumental pa- 
rameters, as for example the beam aperture and the col- 
lection angle (Batson, 1985; Kohl, 1983; MuUer and Sil- 
cox, 1995; Ritchie and Howie, 1988). Particular emphasis 
has been made on taming delocalization, which is rele- 
vant to devise ways of improving the spatial resolution of 
the inelastic signal (Allen et al, 2003; Oxley and Allen, 
1998). We are not going to enter into details of the in- 
struments here, but it is instructive to consider effects 
related to partial detection of the transmitted electrons. 

We can regard each electron in a TEM as consisting of 
a coherent superposition of plane waves, yielding for in- 
stance a narrow focus close to the sample in STEM mode. 
A very accurate approximation can be adopted if we keep 
in mind that the angular aperture of the beam has re- 
mained in the range of a few milliradians since the early 
days of electron microscopy (Ruska, 1987). The compo- 
nents of the incident electron wavevector perpendicular 
to the beam direction z are thus pi± < 10~^pi, where 
Pi is the total wavevector. The parallel components be- 
come = {p^ ~ Pi±.y^^ ~ Pit where we are neglecting 
~ pf^f^Pi ~ 10~^Pi in front of pj. We conclude that the 
incident charge can be reliably described near the sample 
by the wave function 



We are now prepared to recast Eq. (15) into a more 
practical formula. Using delta- function normalization for 
the final states, multiplying the transition rate by the in- 
teraction time L/v, and making use of the above consid- 
erations, we find 

2 dT^^iu) 



fNRc, a — f dT (c 



where 



^ gip/^.(R-R') ^iu,{z'-z)/v ini{-l^(r,r',w)} 

is the loss probability resolved in p/j_ (the lateral 
wavevector of the transmitted electron ± v), or equiv- 
alently, the outgoing direction. This formula indicates 
that STEMs can be used to retrieve the full nonlocal 
dependence of Im {—W{r, r' ,uj)}, but we defer a detailed 
discussion of this point to Sec. V.B.2, in which we explore 
vicinage effects. 

In practice, electron beams are polychromatic. How- 
ever, a trivial extension of the above formalism leads to 
the intuitive result that the loss probability is then the 
average of Eq. (17) over the incident-beam spectrum, 
thus precluding interference between components of dif- 
ferent primary energy. 

Equation (17) leads to a powerful result, established 
by Ritchie and Howie (1988), regarding the validity of 
the classical dielectric formalism employed in Sec. II.B.l. 
The unrestricted integral over all possible values of p / j_ 
yields ^(R — R'), so that the loss probability reduces to 



^ EELS 



M = / 



d^R |V>i^(R) 



rEELs(R'>'^)' 



where rggLg(R, w) is given by Eq. (10). In other words, 

the EELS probability is well described by Poisson's equa- 
tion if all the inelastic signal is collected (i.e., using a wide 
acceptance angle in the spectrometer), but it needs to be 
averaged over electron impact parameters weighed by the 
spot intensity |'i/'j_L(R)p. 



Mr) 



Li/2 



e'f*^Vi±(R), 



(16) 



where L is the quantization length of the trajectory and 
the R dependence refiects the variation perpendicular 
to z. This should be valid for typical sample thickness 
Az ^ 50 nm, along which the lateral divergence of the 
beam can be quantified as AR < 0.5 nm, a small value 
compared to typical dimensions of common photonic and 
plasmonic structures. 

The transmitted electrons can be again described by 
plane waves of wavevector pf. Then, the frequency trans- 
fer reduces to w « q • v under the approximation of Eq. 
(3), where q = — py is the wavevector transfer and 
V = (Ji/me)piZ is the incident electron velocity. 



III. ELECTRON ENERGY-LOSS SPECTROSCOPY 

Hillier and Baker (1944) were the first to propose and 
demonstrate EELS in TEMs, although earlier pioneering 
experiments reported energy losses of transmitted elec- 
trons in thin films (Leithauser, 1904). This technique 
has become standard in the electron microscopy commu- 
nity and is capable of providing information on electronic 
band structures and plasmons in the low-energy loss re- 
gion, as well as atomically-resolved chemical-identity en- 
coded in core losses (Browning et al, 1993). During its 
prolific existence, valence EELS has contributed to fields 
as varied as biochemistry (for example, in the study of 
excitations sustained by nucleic acid bases, reported by 
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Crcwc et al, 1971), interplanetary science (for instance, 
in the explanation of a 5 eV strong absorption feature in 
cosmic dust found by Bradley et al, 2005), and micro- 
electronics (in particular, in the investigation of the re- 
sistivity of CMOS elements performed by Pokrant et al, 
2006). 

Since the early days of EELS, transmission electron 
microscopes have undergone a tremendous series of im- 
provements that permit currently achieving ~ 0.1 eV en- 
ergy resolution for a sub-nanonictcr-sizc elec;tron beam. 
This opens up new vistas in the low-energy loss region, 
such as addressing the optical properties of photonic 
structures with unprecedented spatial detail. 

The excitation spectrum of a typical material is illus- 
trated in Fig. 4 for bulk Ag. The EELS intensity for 
a given energy loss hu) directly reflects the strength of 
specimen modes corresponding to that energy (Egerton, 
1996). The loss spectrum has been traditionally sepa- 
rated into valence- and core-loss regions, with the divi- 
sion between them arbitrarily established at ~ 50 eV. 
With the customary use of thin samples, the main fea- 
ture in an EELS spectrum is the zero-loss peak (ZLP) 
of unscattered electrons, in which unresolved, very-low- 
energy excitations (e.g., phonons) are buried. The inten- 
sity of valence features is over two orders of magnitude 
lower than the ZLP in the 20 nm Ag film considered in 
Fig. 4, while higher-energy core excitations are even less 
probable. This presents a typical scenario in which the 
inelastic signal is superimposed into the tail of the ZLP. 
Removal of the ZLP is thus important, and various so- 
phisticated procedures have been devised that produce 
very reliable results (Lazar et al., 2006; Nelayah et al, 
2007a). 

Valence losses are generally more intense than core 
losses and allow collecting any given amount of inelastic 
signal with less electron irradiation, therefore producing 
minimum damage to the sample. However, the interpre- 
tation of low-energy-loss images is less direct because it 
involves the excitation of delocalized modes (Erni et al., 
2008; Howie, 2003). In this respect, it is useful to rely on 
theoretical simulations, which are now advanced to the 
point of dealing with virtually any sample geometry in a 
predictive way. 

The books of Raether (1980) and Egerton (1996) are 
classic references of EELS that extensively cover these 
subjects and provide detailed considerations on the op- 
eration and principles of TEMs. Several review papers 
present shorter summaries (Brydson, 2001) and collec- 
tions of experimental spectra (Ahn, 2004). Comprehen- 
sive reports on valence EELS (Wang, 1996) and its theo- 
retical analysis (Rivacoba et al., 2000) are also available. 
The present article intends to supplement this field by 
providing a more extended overview of low-energy losses, 
their characteristics and theoretical understanding, and 
several examples of novel application to nanophotonics, 
for which they might become particularly useful. 



A. Space, momentum, and energy resolution 

Nowadays, some commercial TEMs incorporate the 

latest achievements in space and energy resolution. The 
latter has been tremendously improved over the last 
decade by the development of aberration correctors in 
the electron optics, particularly the spherical aberration 
(Uhlemann and Haider, 1998). This has reduced the 
size of the beam spots from the 1-2 nanometers com- 
monly achieved without such correctors down to sub-A 
dimensions (Krivanek et al., 1999), which allows for in- 
stance imaging closely spaced atoms in crystal samples 
using core losses (Batson et al, 2002; Nellist et al, 2004; 
Varela et al., 2004, 2007), and even identifying details of 
chemical bonding (MuUer et al., 2008). However, the in- 
terpretation of atomically-resolved images is sometimes 
diflicult and requires to rely on simulations (Oxley and 
Pennycook, 2008; Wang et al, 2008a). The search for im- 
proved spatial resolution through aberration-correction is 
actually a work in progress, for example through theoret- 
ical analysis for improved incoherent imaging (Intarapra- 
sonk et al., 2008) and through new sub- Angstrom beam 
characterization techniques (Dwyer et al, 2008). 

The factors that limit the spatial resolution of a mi- 
croscope are typically known as delocalization effects 
(Egerton, 2003, 2007; Erni and Browning, 2005). There 
are several of them, of very different nature, that can be 
grouped into three distinct categories: (1) instrumental 
or lens delocalization, which as noted above can be re- 
duced below 1 A; (2) Coulomb delocalization associated 
to the finite range of the field that accompanies a fast 
electron (see Fig. 6 and discussion in Sec. ILA); and (3) 
the extended nature of the excitations that are probed, 
varying from macroscopic distances in the case of low- 
energy SPPs and CR losses, to nanometers in particle 
plasmons, and less than 1 A for core losses. For example, 
sub-atomic resolution is currently limited by the size of 
the probe rather than by the extension of the ionizing 
interaction that is employed to resolve core levels (Allen 
et al, 2003). 

Additionally, spectral resolution has dramatically im- 
proved with the arrival of new spectrometers yielding 
^ 50meV accuracy (Brink et al, 2003), and electron 
monochromators (Su et al., 2003), which roughly consist 
in filtering out incident electrons outside a narrow energy 
window at the expense of reducing the beam current. An 
energy resolution of 0.1 eV has been achieved. Measuring 
excitations down to less than fiuj = 0.5 eV is now possible 
thanks to the limited extension of the monochromatized 
ZLP (Lazar et al., 2006; Terauchi et al., 1999). This 
spectral resolution should be sufficient for studying most 
collective excitations supported by metallic systems, in 
which the intrinsic width produced by absorption is gen- 
erally larger than a tenth of an electronvolt. 

Deconvolution techniques, utilized to eliminate the 
ZLP and enhance spectral resolution, have progressed 
tremendously in recent years (Lazar et al, 2006; van 
Benthem et al., 2001). Furthermore, dynamical insta- 
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bilities in the electron beam energy, which directly dam- 
age spectral details, are now corrected by resorting to 
fast spectra acquisition, well below the 50 — 60 Hz of the 
omnipresent electrical network, and by subsequent addi- 
tion of the collected spectra after repositioning the ZLP 
maximum (Nelayah et al., 2007a). 

A reasonable degree of momentum-transfer resolution 
is possible in TEMs by varying the divergence half-angle 
of the incident beam tpin and the collection half-angle of 
the spectrometer tpout [see Fig. 2(a)]. These parameters 
are typically in the range of a few milliradians. Of course, 
the radius of the electron spot AR, which controls spa- 
tial resolution in STEMs, is related to (pin through the 
uncertainty principle, (nie/tijARipin > l/2v. Actual op- 
eration conditions in TEMs are very close to this limit. 

These impressive achievements are the result of a his- 
torical quest, still in progress, for ever better spatial 
and spectral resolution (Egerton, 2003, 2007; Erni and 
Browning, 2005). However, many factors are remaining 
that limit the performance of electron microscopes, as 
recently reviewed by Egerton (2007). 



B. Bulk losses and determination of bulk dielectric 
functions 

Although the volume of homogeneous materials offers 

little chances for performing microscopy, bulk losses are a 
splendid source of information on optical properties (Pa- 
lik, 1985; Raether, 1980) and they configure an area of 
research in which the electron microscopy community has 
made extremely important contributions. In a homoge- 
neous medium, we can combine Eqs. (2) and (8) to write 
the EELS probabiUty as 



■ bulk 



(^) 



(18) 



X Im 



- - In 



where L is the length of the trajectory and e has been 

assumed to be independent of wavevector q. This is the 
local response approximation, which applies when low- 
enough momentum transfers are collected below a certain 
cutoff 

flQc « V {rrieViponty + {hoj/v)^ 

that is determined by the half-aperture collection angle 
of the microscope spectrometer ipo^t (assuming <^in <C 

font)- 

The non-retarded version (c — > oo limit) of Eq. (18), 

Tbu'lkH = ^^^{^} HQcvH, (19) 

is commonly employed to perform Kramers-Kronig (KK) 

analysis and to retrieve bulk optical constants from mea- 
sured EELS data after absolute calibration of the loss 



function Im{— 1/e} (Egerton, 1996; French et al, 1998; 
Pfliiger and Fink, 1991; Shiles et al, 1980; Zhang et al, 
2008). However, retardation corrections like CR losses 
(see Sec. IV. D) can be very damaging in the determina- 
tion of electronic band gaps and they must be carefully 
removed from the measured spectra in some cases (Erni 
and Browning, 2008; Jouffrey et al., 2004; Stoger-Pollach 
et al., 2006; Stoger-Pollach and Schattschneider, 2007). 
In addition, multiple inelastic scattering events and fac- 
tors depending on beam divergence and spectrometer ac- 
ceptance angles, as well as the precision in the align- 
ment of the sample, have to be properly addressed be- 
fore using KK to retrieve e (Bertoni and Verbeeck, 2008; 
Stoger-Pollach, 2008). Separating volume and surface 
losses adds another complication, which can be solved by 
comparing spectra acquired from specimens of different 
thicknesses (Mkhoyan et al., 2007). 

One is often interested in opening the collection an- 
gle to increase the inelastic signal and to reduce sample 
damage. Nonlocal effects are then apparent, for exam- 
ple through e-h pair excitations, as represented in Fig. 
8(a) for a free-electron gas of parameters corresponding 
to Al. For instance, this region is accessible when using 
200 keV electrons and LpoMt = lOmrad, a combination 
that results in the qc value indicated by the heavy arrow 
of Fig. 8(a). Under such circumstances, we need to in- 
clude spatial dispersion in the dielectric function, which 
has different forms for longitudinal and transversal fields, 
eion(9, w) and etr(g, w), respectively, as pointed out in Ap- 
pendix A. The momentum-resolved loss probability is 
then found to be 



<^rbuik('^) 



X Im 



1 



(20) 



q^ 



1 



1 



etr(g,w) 

1 

etr(?,w) 



i(a,w) 
1 



^2 - k'^etT{q,uj) 



where q_L is the 2D wavevector transfer perpendicular to 
the electron trajectory, and the total wavevector transfer 
satisfies q"^ = q±+ co'^/v'^ in virtue of Eq. (3). 

The momentum dependence of the transversal permit- 
tivity can be safely neglected in this discussion, since it 
constitutes a correction of order vp/c to the response, 
where vp is the Fermi velocity (e.g., vp = 1.39 x 
10* cms^^ in gold). Besides, the Lindhard formula (Lind- 
hard, 1954) for the response of a free-electron gas (see ex- 
plicit expression in Sec. A.l) provides a fair description of 
the longitudinal dielectric function in good metals such 
as Al, although local-field corrections can be also relevant 
(Vast et al., 2002). 

In this context, the response of gold has increas- 
ing technological importance, but detailed experimental 
analysis of the momentum-dependent optical constants is 
still missing. The local dielectric function determined by 
Johnson and Christy (1972) from ellipsometry measure- 
ments is routinely employed in plasmonic studies, and it 
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FIG. 8 (Color online) Bulk excitations in Al. (a) The con- 
duction band of Al is well described by a degenerate free- 
electron gas of Fermi wavevector fcp — 1.75 A~^. The e- 
h pair excitations supported by this system are contained 
in the shaded region within the plotted wavevector-energy 
{q — hui) diagram. The volume plasmon dispersion rela- 
tion measured by Batson and Silcox (1983) (symbols) is 
compared to the longitudinal plasmon-pole approximation 
Lo = [uj'^ + {h/m^f{/3^q^ + g*/!)]^''^ (thick solid curve) de- 
rived from the condition eion{q,u!) = 0, as obtained from Eq. 
(A6) with P = (3/5) ^^^^F and et = 1. (b) Logarithmic-scale 
representation of some the curves in (a) , compared to the light 
line (tj = qc) and the transversal bulk plasmon dispersion re- 
lation {uj = [ijj'p + q^c^]^^^). The latter, which is obtained by 
inserting Eq. (A7) into q^ = k^etr, gives rise to light emission 
in thin films irradiated by electrons (Ferrell, 1958; Vincent 
and Silcox, 1973). 



can be well described by the Drude-like formula of Eq. 
(A7) in the visible and NIR part of the spectrum (see Ta- 
ble II). In a recent report (Garcia de Abajo, 2008), the 
author has proposed a nonlocal eion for gold constructed 
as 



eion(g,w) 



\Lo)^e''iu;) + e'^{q,Lo), (21) 



where e and e are the Drude and Mermin dielectric 
functions given by Eqs. (Al) and (A4), respectively, and 
corresponding to the values of and rj listed in Table II, 
whereas e°^P is the measured local permittivity. Further 
experimental effort is still required to contrast this model. 

The spectrum in bulk metals is dominated by plasmon 
losses, first observed by Ruthermann (1948) and Lang 



(1948), and later identified as collective oscillations of 
the conduction-electron gas by Pines and Bohm (1952). 
These modes are signaled by the vanishing of eion, as 
deduced from the non-retarded (NR) limit of Eq. (20), 



dqj_ 



Im 



eion('7, w) 



(22) 



The dispersion relation of bulk plasmons has been ex- 
tensively studied in the past using EELS (Batson et al, 
1976; Chen et al., 1980, 1976; Raether, 1980; Watanabe, 
1956) and is shown in Fig. 8(a) for Al, as measured by 
Batson and Silcox (1983) (symbols). When broad col- 
lection angles are considered, one can reliably approxi- 
mate £1011(97 w) by a plasmon-pole analytical expression 
(Ritchie, 1957), reproduced in Eq. (A6). The plasmon 
dispersion relation derived from this formula provides a 
reasonable interpolation between the measured plasmon 
and the Bethe ridge lu — Tiq^ /2mi,, shown in Fig. 8(a) as 
a solid curve. 

Spatially-resolved bulk-plasmon spectroscopy has been 
used to identify the presence of metals in a specimen [e.g., 
Na in either Na20 or silicate glasses (Jiang et al., 2008)] 
and to demonstrate the effect of quantum confinement in 
carbon nanostructures by tracing local variations in the 
plasmon energy (Stolojan et al., 2006), which are related 
to changes in the conduction-electron density through 
Eq. (A2). 

Finally, magnetic circular dichroism has been mea- 
sured in the response of iron at x-ray energies 
(Schattschneider et al., 2006), suggesting that further 
work could eventually explore optical activity at visi- 
ble and NIR frequencies via EELS, for instance in chiral 
(Rogacheva et al., 2006) and left-handed (Boltasseva and 
Shalaev, 2008) metamaterials. 



C. Planar surfaces 

Surfaces and interfaces can host trapped modes that 
modify the local band structure and the optical response 
with respect to bulk materials. For instance, we find 
electronic surface states confined by directional gaps in 
noble-metal interfaces (Smith, 1985), but also specific 
collective modes of conduction electrons (e.g., surface 
plasmons). As we show later on, these features can be 
characterized by EELS performed in STEMs with high 
degree of spatial resolution. 

A particularly instructive situation is presented when 
the beam is directed parallel to a planar interface sep- 
arating two different media (Fig. 9, left). This configu- 
ration has been extensively studied both experimentally 
(Lecante et al., 1977; Moreau et al., 1997; Powell, 1968; 
Walls and Howie, 1989) and theoretically (Echenique and 
Pendry, 1975; Garci'a-Molina et al., 1985a; Lucas and 
Sunjic, 1971). The analysis of this geometry in the non- 
retarded limit can be readily made from the knowledge 
of the screened interaction W{y,y' ,uj), defined as the po- 
tential created at point r by a charge oscillating with 
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FIG. 9 (Color online) Planar interface separating two media 
of dielectric functions ei and e2, respectively. An electron 
is shown moving parallel to the boundary inside medium 1 
(left), with impact parameter b relative to the interface. An- 
other electron is represented crossing the interface at normal 
incidence (right). 



frequency lj and placed at r'. Expressions of W for sim- 
ple geometries, including that of Fig. 9, are given in Ap- 
pendix B. A closed- form expression for the loss probabil- 
ity of the electron following the parallel trajectory of Fig. 
9 (left) is then obtained from Eqs. (10) and Eqs. (B3). 
One finds 
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(23) 



£1 + £2 



Im 



£1 



where L is the length of the trajectory. The impact pa- 
rameter b and the dielectric functions ei and £2 are de- 
fined in Fig. 9. The first term inside the curly brackets is 
the bulk loss probability of Eq. (19). The square brackets 
in Eq. (23) contain the effect of the interface, which dies 
off exponentially with b as Kq{9) w exjp(—9) ^/tt~/29 for 
large = 2ujb/v (Abramowitz and Stegun, 1972). The 
first interface term describes the excitation of intrinsic 
boundary modes, signalled by the condition £1 + £2 = 0. 
It was first derived by Echenique and Pendry (1975), who 
obtained the expression 
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for an electron moving in vacuum near a medium of di- 
electric function e. The second interface term in Eq. (23) 
accounts for a reduction of bulk losses (i.e., the trans- 
fer of oscillator strength from volume to surface modes) . 
This is the so-called begrenzungseffekt. 

Equation (23) has been successfully applied to explain 
surface losses suffered by electrons passing near MgO 
cubes (Cowley, 1982a,b; Marks, 1982; Walls and Howie, 
1989; Wang and Cowley, 1988) in the so-called aloof con- 
figuration (see below), and also to EELS in a Si/Si02 
interface and in GaAs surfaces (Howie and Milne, 1985). 

Inclusion of retardation effects is critical to account for 
CR losses and to correctly assess the weight of interface 
effects, as demonstrated by Moreau et al. (1997) for the 



Si/Si02 interface, and later by Couillard et al. (2007) 
for multilayer structures containing HfO. These studies 
have been recently extended by Yurtsever et al. (2008) 
and Couillard et al. (2008), who have measured EELS 
in slabs and layered stacks formed by Si and Si02 using 
electron beams parallel to the interfaces and have demon- 
strated the important role of retardation corrections for 
a quantitative comparison with theory. 

The retarded counterpart of Eq. (23) has been reported 
by several authors (Garcfa-Molina et al., 1985a; Otto, 
1967). A particularly simple derivation consists in ex- 
panding the integrand of Eq. (2) into p (TM) and s 
(TE) electromagnetic plane- wave components, which are 
later reflected at the interface to act back on the electron 
(Forstmann et al., 1991). One finds 



EELS 



(lo) = rbulk(w) + Frcf (w), 



where Fbuik is the bulk loss inside the medium 1 in which 
the electron is moving, given by Eq. (18), and 
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is the loss due to reflect ion of the electron fi eld at the 
interface. Here, q\\ = ^ LtP- jv^ + q^, q^j = ^Jk'^^j — <?jj, 
and 



£2 qzi - £1 qz2 
£2 qzi + £1 qz2 



(26) 



and 



qzi - qz2 
qzi + qz2 

are Fresnel reflection coefficients for p and s waves, re- 
spectively. In the non-retarded limit, one has = 
and Tp = (e2 — £i)/(£2 + £1), from which Eq. (23) is eas- 
ily recovered. The advantage of this approach, based 
upon optical reflection coefficients, is that it can handle 
more complicated surfaces, like the periodically corru- 
gated boundary of a confined photonic crystal (Garcia 
de Abajo and Blanco, 2003; Pendry and Martm-Moreno, 
1994). For large impact parameters, Eq. (25) provides 
a retardation correction to Eq. (24), consisting in sub- 
stituting 2ujb/vj for 2ujb/v in the argument of the Kq 
Bessel function; the apparent impact parameter at very 
large velocity is contracted to 6/7. 



1. Excitation of surface plasmons and surface-plasmon 
polaritons 

Surface modes, signalled by the condition £1 + £2 = 0, 
are characteristic of the interface between a dielectric 
having Re{ei} > and a metal satisfying Re{e2} < 0. 
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FIG. 10 First observation of surface plasmons in the energy 

loss distribution of 2020 cV electrons specularly reflected on 
an Al surface under 45° incidence (adapted from Powell and 
Swan, 1959). 



These excitations were first identified and understood 
by Ritchie (1957) to explain anomalous valence losses 
that were previously observed at positions differing from 
bulk plasmons in the spectra of fast electrons transmit- 
ted through thin foils (see, for instance, feature D in Fig. 
2 of Watanabe, 1956). Ritchie's description in terms of 
a classical hydrodynamic plasma to represent the con- 
duction electron band led Stern and Ferrell (1960) to 
baptize these excitations as surface plasmons, which were 
later confirmed in electron energy-loss experiments (Pow- 
ell and Swan, 1959) (see Fig. 10). Bulk and surface modes 
were observed as fiLUp w 15 eV and hws ~ 10.6 eV energy 
losses in electrons reflected from an Al surface. [Inci- 
dently, the spectrum of Fig. 10 contains multiple plas- 
mon losses (i.e., two or more plasmons being excited by 
the same electron). The coherent aspects of double plas- 
mon excitation have deserved some attention in the past 
(Schattschneider et al., 1987).] 

Aluminum is a prototypical example of nearly-free elec- 
tron gas for which the dielectric function is well described 
by the Drude formula of Eq. (Al). The SP of the bare 
surface, satisfying the condition e + 1 = 0, has then a 
frequency lJs = LOp/V'2, in excellent agreement with ex- 
perimental observations. 

Surface plasmons were also observed for other mate- 
rials in subsequent reflection experiments performed on 
melted metals (Powell, 1968), and even the evolution of 
these modes through a solid-liquid phase transition was 
reported (Powell, 1965). Other pioneering electron trans- 
mission experiments managed to demonstrate the obser- 
vation of SPs as a characterization technique (Raether, 
1967) and to determine the dependence of SP energies on 
the thickness of thin films (Boersch et al., 1966a). Ev- 
idence of surface and bulk plasmon excitation has been 
also observed in the spectra of photoelectrons and Auger 
electrons emitted from metals (Osterwalder et al., 1990; 
Simonsen et al., 1997; van Attekum and Trooster, 1978). 

Similar to the bulk plasmons of Fig. 8(a), SPs are af- 
fected by nonlocal effects and exhibit significant disper- 
sion, which was already reported in early EELS studies 
(Krane and Raether, 1976). This has stimulated a rich 



literature intended to achieve a reasonable first-principles 
description of the dynamical response of crystal surfaces 
(see Pitarke et al., 2007, and references therein). Nonlo- 
cal effects involve comparatively short-range interactions 
extending a distance ~ l/g for typical wavevector trans- 
fers q in the A~-^ domain [see Fig. 8(a)]. However, we are 
mainly concerned here with longer-range optical excita- 
tions (e.g., SPPs), in which space dispersion can be gen- 
erally overlooked, except when sharp metallic structures 
or narrow gaps between metals are involved (Garcia de 
Abajo, 2008), as shown in Sec. V.E for electrons passing 
near metal-sphere dimers. Nonlocal effects can be also 
significant at small separations ^ 1 nm between the pass- 
ing electron and metal surfaces (Zabala and Echenique, 
1990). 

An interesting scenario is presented by the aloof con- 
figuration, in which the electron trajectory does not even 
intersect the sample. This leads to distant interaction, 
capable of providing surface-specific information, free of 
close-encounter events between the fast probe and the 
target atoms. The exponential decay of surface losses 
with 2u}b/v'y indicates that low-energy transfers are fa- 
vored in the aloof configuration. This has the additional 
benefit of minimizing sample damage. 

The aloof configuration was pioneered by Lecante et al. 
(1977), who forced the electrons to describe parabolic 
trajectories, defiected from the metal surface by a bias 
potential. Subsequent studies managed to aim a STEM 
aloof beam parallel to the planar surfaces of MgO cubes 
(Cowley, 1982a, b; Marks, 1982). Loss spectra were later 
recorded after transmission through perforated metallic 
channels (Warmack et al, 1984). The aloof configura- 
tion has been recently revisited and compared with near- 
field optical microscopy (Cohen et al., 2003, 1998, 1999; 
Echenique et al., 1999; Itskovsky et al, 2008). 

Retardation adds another source of plasmon disper- 
sion. Since SP modes are resonances in the surface re- 
sponse, they must involve a divergence in the refiectivity 
for incident evanescent waves with the right values of 
(g||,a;). More precisely, we obtain the plasmon disper- 
sion relation from the vanishing of the denominator of 
[Eq. (26)], leading to 

,f = fc.pg^ (27) 

" V ^1 + ^2 

under the condition that Re{ei} and Re{e2} have op- 
posite signs. The dispersion relation of Eq. (27) is rep- 
resented in Fig. 11 (left) for a silver- vacuum interface. 
When the retarded nature of SPs becomes important, 
they are usually referred to as SPPs, because they ac- 
quire a polaritonic character refiecting a mixture of col- 
lective electron motion and propagating electromagnetic 
fields. This tends to happen at low wavevectors (i.e., at 
long wavelengths according to Fig. 11), for which SPPs 
increasingly approach the light line (w = q\\c). It is im- 
portant to stress that this polaritonic regime occurs for 
momentum transfers that are well below the region ex- 
hibiting nonlocal effects [cf. horizontal scales in Figs. 8(a) 
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and 11]. This is shown in Fig. 8(b), where the hybridiza- 
tion of transversal bulk plasmons with light takes place 
at values of q that are two orders of magnitude smaller 
than the Fermi wavevector kp- More details on SPPs, in- 
cluding a comprehensive literature survey, can be found 
in the book of Raether (1988). 

The electric and magnetic fields associated to a SPP 
plane wave moving along x have similar strength in Gaus- 
sian units. For instance, in the vacuum region outside a 
semi-infinite metal of permittivity e, we have 

ESP(r,a;) = [(xgf - zqf)/k] exp[i(5f x + qfz)], 
HSP(r,u;) = y exp[i{qf x + qf z)], 



where gf^ = y fc^ — (g|P)2 = ^k/^/e+ 1. In particu- 
lar, the electric field satisfies V • E^^ = 0, and there- 
fore it is transversal.^*^ Close to the plasmon cutoff, the 
SPP momentum becomes increasingly large. This marks 
the transition towards the non-retarded regime, in which 
magnetic and electric fields are decoupled, so that we 
can approximate V x E^^ ~ and define a potential 
(j) = exp[i(a; — 2;)g|^] to write E^^ = — V^. The electric 
field appears to be longitudinal, and the plasmons are 
then known as SPs rather than SPPs, signalled by the 
condition e = — 1. 

An aloof electron like that of Fig. 9 (left) couples to 
SPPs of momentum component along the velocity vector 
given by Klu/v, according to Eq. (3). Consequently, the 
electron only excites parallel momenta above that value. 
This condition is represented in Fig. 11 (left) for two dif- 
ferent electron energies. The corresponding electric fields 
are shown in the same figure (right) within a plane that 
contains the velocity vector and the surface normal. The 
electron is passing at a distance of 10 nm above the sur- 
face. We can see several characteristic features in these 
plots: (1) the field shows a wake pattern characterized 
by oscillations of wavelength « 27r/g||, with g'n = w/v de- 
termined by the intersection points of Fig. 11 (left); (2) 
the electron that moves faster excites plasmons of lower 
energy, thus giving rise to oscillations of longer spatial 
period; (3) at variance with the continuity of the parallel 
electric field, the normal component changes sign across 
the surface in order to preserve the continuity of the nor- 
mal displacement, because the metal permittivity is neg- 
ative; (4) the normal field component takes large values 
compared to the parallel one, and this effect is more pro- 
nounced for faster electrons, which involve lower co's. 

In a different direction, the early 1990s witnessed re- 
markable experiments of energy losses in coincidence with 



Although the main E^^'-component of long-wavelength SPPs is 

along the surface normal, there is also a finite electric field along 
the propagation direction, but the normal components of both 
the wavevector and have a 7r/2 phase difference relative to 
the parallel components, leading to V • E^^ = 0. 



SEE, initially collected in amorphous carbon samples 
(Miillejans and Bleloch, 1992; Pijper and Kruit, 1991). 
Angle-resolved measurements were also carried out in 
this context to conclude that the SEE yield is larger in 
more localized excitations, involving larger-momentum 
transfers (Drucker and Scheinfein, 1993). Comparative 
studies for amorphous carbon and silicon suggested that 
bulk plasmons do not play a central role in SEE (Drucker 
et ai, 1993; Scheinfein et ai, 1993), possibly due to the 
noted begrenzungseffekt, which limits the strength of vol- 
ume plasmons in favor of SPs within the region accessible 
to the escape depth of SEs. Finally, coincidence experi- 
ments were conducted using the aloof beam geometry in 
diamond and MgO to compare EELS and SEE rates and 
to provide direct measurements of the probability that 
a surface excitation gives rise to a SE, which was found 
to be below 5% for low-energy losses (Miillejans et al., 
1993). 



2. Guided modes in thin films 

The SPPs on both sides of a thin film interact to pro- 
duced two hybridized modes. One of them can travel 
longer distances along the film due to exclusion of the 
electric field from the metal (Sarid, 1981). This mode 
is currently being applied to propagate electromagnetic 
signals in plasmonic devices (Barnes et al, 2003; Berini 
et ai, 2007). The dispersion relation of coupled SPPs in 
thin films was obtained from the EELS signal of transmit- 
ted electrons over 30 years ago by a series of outstanding 
experiments conducted by Silcox and co-workers in Al 
(Pettit et ai, 1975; Vincent and Silcox, 1973) and Si films 
(Chen et al., 1975). An example of these studies is given 
in Fig. 12 for a partially oxidized thin Al film. Both plas- 
mon branches follow the light line at low energies, and 
they converge to the non-retarded limit of the AI/AI2O3 
interface for large momentum transfer. This limit cor- 
responds to EAi + 4 = 1, or equivalently, TiLOg ~ 6.7 eV 
when Eq. (Al) is used to model Al, since the permittiv- 
ity of AI2O3 in this energy window is « 4. Large gy's 
involve fast plasmon oscillations along the film, so that 
alternating induced charges weakens the interaction be- 
tween both film sides and we recover the limit of the 
single interface separating two semi-infinite media. For 
small gy, the interaction between plasmons in both sides 
of the film gives rise to two SPP branches. 

Guided modes in graphite (Chen and Silcox, 1975a) 
and aluminum-oxide (Chen and Silcox, 1975b) films were 
also characterized using the same technique. In this case, 
the film behaved as a dielectric, and the electron coupled 
to optical Fabry-Perot resonances set up by successive 
total internal reflections at the film boundaries. The ex- 
periments were in excellent agreement with theory (Chen 
and Silcox, 1979). 

In a more recent development, a complex multilayer 
structure containing HfO, which is of interest in micro- 
electronics, has been characterized by Couillard et al. 
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FIG. 11 (Color online) Surface plasmon dispersion relation in Ag (left) and excitation of SPPs by swift electrons (right). 
The SPPs saturate at Ritchie's non-retarded frequency for large parallel wavevector and follow the light line u) = q\\c in the 
low-energy limit. An electron moving with velocity v parallel to the surface preferentially excites plasmons of frequency and 
wavevector related by a; < q\\v. This condition is indicated by solid dots in the dispersion relation for the two electron energies 
considered on the right. Faster electrons generate plasmons of lower q\\ and longer wavelength. The impact parameter is 10 nm 
in the field plots. The dielectric function of Ag is taken from Palik (1985). 
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FIG. 12 (Color online) Experimental determination of the 
SPP dispersion relation in a thin film (adapted from Pettit 
et al, 1975). Symbols: measured data obtained from energy 
and angle distributions of fast electrons traversing a partially 
oxidized Al thin film. Continuous curves: SPP branches cal- 
culated from optical data. 



(2007) . The agreement with theory is excellent after re- 
tardation effects are incorporated to correctly account for 
interface plasmons and Cherenkov guided modes. A sim- 
ilar study has been recently presented by Yurtsever et al. 

(2008) in slabs and stacks of Si and Si02. 
Additionally, a radiative plasmon branch inside the 

light cone, giving rise to losses near ujp, was predicted 
by Ferrell (1958) for thin metal films, and the resulting 
light emission was detected later by Steinmann (1960) in 
Ag. We discuss this issue in more detail in Sec. IV.F.l. 



D. Curved geometries 

The ability of TEMs to probe local response proper- 
ties has been extensively exploited to study non-planar 



sample geometries. Particular attention has been paid 
to spherical and cylindrical shapes, carbon nanotubes, 
composites, and other more complicated structures. The 
intense valence losses can provide spectral images, the in- 
terpretation of which is not as direct as when core losses 
are employed instead. Therefore, we need to strongly rely 
on theory, which is actually in a very advanced state, ca- 
pable of dealing with virtually any system. For example, 
small particles like fuUerenes or carbon nanotubes are 
reliably described by atomistic models, although first- 
principles theory is also possible (Marinopoulos et al., 
2003). For larger structures (e.g., nanoparticles) we have 
a suite of methods based upon classical dielectric the- 
ory, which gives excellent agreement with experiment in 
metallic objects of dimensions above ~ 10 nm. In this 
section, we explore some of these structures possessing 
spherical or cylindrical geometry. 



1. Cylinders 

Early work on cylindrical nanocavities relied on an 
electron beam to drill holes in AIF3 (Macaulay et al., 
1989), which were investigated by EELS (Scheinfein 
et al., 1985). The synthesis of such types of holes has 
considerably improved over the last few years in two 
different directions. (1) Bottom-up approach: chemical 
methods have been developed that are capable of sponta- 
neously forming arrays of self-organized nanoholes in ma- 
terials like alumina (Masuda and Fukuda, 1995); arrays 
of cylinders are also formed in eutectics (Pawlak et al, 
2008). (2) Top-down approach: advanced electron-beam 
and focused-ion-beam (FIB) lithographies have been ex- 
tensively employed to produce hole arrays for nanopho- 
tonics studies (Genet and Ebbesen, 2007). These and 
similar geometries have been investigated through TEM 
spectroscopic analysis (Degiron et al., 2004; Garcia de 
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FIG. 13 (Color online) Notation used for homogeneous 
spheres and cylinders or radius a, with an electron passing 
at a distance b from the center. 



Abajo et ai, 2003). 

Likewise, self-standing Si (Reed et ai, 1999) and Ge 
(Hanrath and Korgel, 2004) nanowires have been stud- 
ied to show clear evidence of the begrenzungseffekt in 
the former and confinement effects in the position of the 
bulk plasmons in the latter. Interface losses in Bi wires 
embedded in an alumina matrix have been also observed 
(Sander et ai, 2001) and theoretically explained (Garcia 
de Abajo and Howie, 2002). 

The axial symmetry of a cylinder enables us to clas- 
sify its modes according to their azimuthal dependence 
expiirmp), whereas a separation into expiiq^z) compo- 
nents of wavevector Qz is appropriate for the variation 
along the direction of translational invariance, chosen 
along z. This decomposition is employed in Sec. B.2 to 
write the non-retarded screened interaction, from which 
the cylinder modes are derived using the condition A,„ — 
[see Eq. (B4)]. A relatively simple expression can be 
obtained for these modes including retardation (Ashley 
and Emerson, 1974): 

xfxj [e2Xll'ra{x2) / Im{x2) - eiX2K'^{xi) / Km{xi)] 
X [Xil'^{x2) / Im{x2) - X2K'^{xi) / K^{xi)\ 

^m\e,-e2f{kqy)\ (28) 

where the notation of Fig. 13 has been adopted, the sub- 
scripts 1 and 2 refer to the media outside and inside the 
cylinder, respectively, and Xj = {u}a/v)yjl — e^u^/c^. 

The interaction of the electron with a cylinder has been 
described using both non-retarded (Tu et ai, 2006; Za- 
bala et ai, 1989) and fully-retarded (Walsh, 1991; Zabala 
et ai, 1989) dielectric theory. One can easily derive the 
loss probability using the screened interaction W in the 
non-retarded limit [see Eq. (B4)]. For electrons moving 
parallel to a cylinder, one just needs to insert W inside 
Eq. (10). The momentum transfer along the cylinder 
axis, is fully determined by the condition qz — w/v, and 
one obtains the loss probability 

Tll^y,{Ro,u;) ^ ^ lm{-W^{b,b,qz,u;)}, 

m— — oo 

where b is the distance from the electron beam to the 
cylinder axis, and different m's contribute separately. Re- 



tardation effects are however tedious to deal with analyt- 
ically (Walsh, 1991; Zabala et al., 1989), as compared to 
the planar interface. For instance, a beam moving along 
the center of the cylinder, which only couples to m = 
modes due to the symmetry of the electron-cylinder com- 
bined system, experiences a loss probability 

^ e2Xi/o(a;i)/i(x2) - eiX2li{xi)Io{x2) \ 
e2XiKo{xi)Ii{x2) + eiX2Ki{xi)Io{x2) j ' 

where we can immediately see the presence of the 
TO = mode predicted by Eq. (28) when we 
arrange the denominator of the last fraction as 

Ko{xi)Io{x2){e2XilQ{x2) / Io{x2) - '^1X2Kq{xi)/Ko{xi)). 

2. Spheres 

The interaction of electrons with spheres has attracted 
considerable attention since the publication of pioneer- 
ing experimental loss spectra in alkali-halide (Creuzburg, 
1966) and Al (Fujimoto et ai, 1967) particles, followed 
by their subsequent theoretical interpretation (Crowell 
and Ritchie, 1968; Fujimoto and Komaki, 1968). Further 
EELS measurements corroborated these results and fo- 
cused on the role of plasmon dispersion in small spheres 
(Batson, 1980; Ouyang et ai, 1992) and in the effect of 
the interaction between neighboring particles (Batson, 
1982b, 1985; Ugarte et al, 1992). More recently, UV 
surface exciton polaritons have been observed in gold 
nanospheres (Chu et al, 2008) and whispering gallery 
modes measured in silica beads (Hyun et ai, 2008) using 
EELS. 

A small metallic sphere embedded in a dielectric host 
exhibits plasmon modes at frequencies dictated by the 
condition (Kittle, 1987) 

{I + l)ei + 162 = 0, (29) 

where I labels the multipolar order of the excitation, and 
ei and £2 are the permittivities of the host and the metal, 
respectively (see Fig. 13). In particular, for an Al sphere 
in which 62 is described by the Drude dielectric function 
of Eq. ( Al) , the frequency of the non-retarded plasmon 
oscillations becomes 

For example, the dipole {I — 1) of an Al sphere in vac- 
uum has energy Tiujp/y/S ~ 8.7 eV. For very large /, the 
surface-charge oscillations are so dense that they cannot 
resolve the particle curvature, and therefore the mode 
frequency evolves to the plasmon of the planar interface, 

UJs — Uj/y/ei + 1. 

The remarkable works of Fujimoto and Komaki (1968) 
and Crowell and Ritchie (1968) established the basis 
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for the theoretical understanding of fast-electron inter- 
action with spherical particles. These authors described 
a metallic sphere using the hydrodynamical model of Ap- 
pendix A. 2 and derived analytical expressions for the loss 
probability and photon emission rate of a broad beam ne- 
glecting retardation. Their results showed spectral fea- 
tures at frequencies roughly described by Eq. (30). 

This theoretical development was followed by studies 
of voids in metals (Ashley and Ferrell, 1976), and later by 
the derivation of an analytical formula for the loss prob- 
ability of a STEM electron passing outside a sphere (Fer- 
rell and Echenique, 1985; Ferrell et al, 1987). Namely, 



1999c; Garcia de Abajo and Howie, 1998): 
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is the non-retarded multipolar polarizability of a sphere 
of radius a in vac;uum. and b is the distance from the 
trajectory to the sphere center, as illustrated in Fig. 13. 
Equation (31) can be derived upon insertion of the non- 
retarded screened interaction near a sphere [Eq. (B6)] 
into Eq. (10). This formula successfully explained pre- 
vious measurements for aloof electrons (Batson, 1982a) 
and was soon applied to explain EELS of spheres made 
of tin (Acheche et al., 1986), silicon (Ugarte et al., 1992), 
silver (Ouyang et al., 1992), and aluminum (Stockli et al, 
1997), as well as corrugated metal surfaces (Mamola 
et al., 1987), and diamond particles (Fehlhaber and Bur- 
sill, 1999, 2000; Peng et al, 2001). The effect of coatings, 
which can be important due to oxidation (Ferrell et al, 
1987; Ugarte et al, 1992), is described by Eq. (31) if we 
use the polarizability of a coated sphere instead of Eq. 
(32). 

Optical modes in nanoparticles (e.g., plasmons in 
metallic spheres) are well-known to be redshifted due 
to retardation when the particle size is not negligible in 
front of the wavelength (see Sec. IV. F. 2): the interaction 
between distant regions of a particle is mediated by elec- 
tromagnetic signals that propagate at the speed of light, 
but the signal simply takes longer to travel over larger 
distances as the particle dimensions increase, so that the 
period of the oscillation modes decreases, thus produc- 
ing the noted redshift. This effect was first described by 
Mie (1908) in his analysis of light scattering by metallic 
spheres. Similar phenomena occur in particles made of 
polaritonic materials (Fuchs and Kliewer, 1968). 

The extension of Eq. (31) to incorporate retardation 
also leads to an analytical formula (Garcia de Abajo, 
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where the Lorentz factor 7=1/ \/l — v^/c^ accounts for 
the contraction of the impact parameter b as seen by the 
moving electron. Here, and are coupling coeffi- 
cients that depend exclusively on the v/c ratio, and tf^ 
and tf are magnetic and electric Mie scattering coeffi- 
cients independent of beam energy and impact parame- 
ter. Analytical expressions for these coefficients are given 
in Appendix C for homogeneous spheres. Equation (33) 
results from a multipolar expansion of the electric field 
induced by the electron [Eq. (4)] after each multipole is 
scattered by the sphere to produce an outgoing wave act- 
ing back on the electron and giving rise to energy loss [see 
Eq. (8)]. 

It should be emphasized that Eq. (33) can be applied 
to any spherical object for external electron trajectories, 
provided tf^ and tf are known for that object. In partic- 
ular, this has been recently done for left-handed spheres 
(Xu and Zhang, 2008) and spherical metamaterial cloaks 
(Xu et al, 2008). 

Equation (33) predicts sizable redshifts even for rela- 
tively small Al spheres of radius a = 10 nm, as shown in 
Fig. 14(a) by comparison to spheres with a = 5 nm (solid 
curves). Spectral features associated to dipolar {I = 1) 
and quadrupolar (I = 2) excitations are clearly discern- 
able in the spectra. Higher-order multipoles produce a 
shoulder piling up near Us because their separation is 
smaller than the absorption width. Besides, the dipole 
feature is broader in the bigger sphere, which is again 
a retardation effect, since the coupling of this mode to 
emitted light increases with the particle dipole strength 
(i.e., with particle size), therefore limiting the lifetime of 
the dipolar plasmon due to radiative decay (Myroshny- 
chenko et al, 2008b). This point is further discussed in 
Sec. IV.F.2. 

For sufficiently small particles and aloof trajectories, 
we can just retain dipole terms {I = 1), in which case Eq. 
(33) reduces to 



r ( ^ ^ (2eu\\(u:b 

i EELS,dip(.WJ = — / 



where 



m = Ki{Q + Ki{c:)h\ 



Im{a}, (34) 



(35) 



Actually, this approximation works quite well even for 
relatively big spheres, as shown by the dotted curves of 



Loss and CL emission probabilities are plotted per unit of lost- 
energy range (i.e., in eV^"*"), so that the areas under the spectral 
curves of Fig. 14(a) yield dimensionless probabilities per electron. 
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FIG. 14 (Color online) Plasmon losses and spatial resolution 
in Al spheres, (a) The loss probability is calculated for two 
different particle sizes, exhibiting several multipole plasmons 
that redshift as the particle radius increases. Solid curves: 
result of Eq. (33) using local response. Dashed curves: same 
equation with nonlocal scattering matrices (see Appendix C). 
Dotted curves: prediction of Eq. (34). The Al is described 
by the dielectric function of Eq. (Al) with Tioj = 15 eV and 
fiT] = 0.5 eV. (b) Impact-parameter dependence of the dipolar 
mode maximum, (c) Apparent radius of the dipole excitation 
in an Al sphere, modeled within local response, and defined 
as the impact parameter 61/2 for which the EELS probability 
of a 200 keV electron is half the maximum as a function of h. 
The dipole-energy dependence on particle size (right scale) is 
taken into account. The V'y/2uj distance of exponential decay 
in Feels, dip [Eq. (34)] is shown for comparison. 



Fig. 14(a), obtained from Eq. (34) with the electric po- 
larizabihty 



(36) 



[see Eq. (C2); tf is a Mie scattering coefficient]. The 
dipole mode is given correctly, although higher multi- 
poles are obviously missed. In contrast to the particle- 
size-independent dipole-plasmon energy predicted by Eq. 
(32), we find that Eq. (36) incorporates retardation ef- 
fects leading to a sizeable dipole-plasmon redshift with 
increasing particle diameter. 

Spheres offer an excellent playground to discuss 
Coulomb delocalization effects (see Sec. III. A). As antic- 
ipated by Eq. (4) for the external electron field, Eq. (33) 
predicts an exponential decay of the EELS signal with 
b, and an attenuation distance given by v^/2uj. This is 
true when the argument of the Bessel functions is large 
compared to unity, a condition that is not satisfied near 
small spheres, in which those functions have a steeper de- 
pendence on b (see Sec. II. A). This is actually observed 
in the b dependence of the loss probability shown in Fig. 
14(b) (solid curves) for the dipole peaks of the spheres 
considered in Fig. 14(a). The probability falls down by 
a factor of two at a distance of ~ 2 nm away from the 
surface of the smaller sphere, so that its apparent radius 
bi/2 ^ 7nm is significantly smaller than v^jluj ^ 11 
nm. The sphere size dependence of b^j^ [Fig. 14(c)] fully 
corroborates these arguments. 

Small metallic particles exhibit further peak shifts orig- 
inating in quantum confinement and spatial dispersion, 
the effect of which is to break down the hypothesis of 
local response. Nonlocal effects in the performance of 
bounded metals can be easily included using the specular- 
reflection-model approximation (Garcia de Abajo, 2008; 
Ritchie and Marusak, 1966; Wagner, 1966). Then, Eq. 
(31) is still valid if we replace e by (Rojas et al, 1988) 



2a{2l+l) 



dq 



e(9,^) 



(37) 



in the polarizability of Eq. (32), where e{q,uj) is the non- 
local bulk response discussed in Sees. III.B and A.l, and 
ji is a spherical Bessel function. This level of description 
is similar to the work of Fujimoto and Komaki (1968). 
Under the conditions of Fig. 14(a), one needs however to 
simultaneously include nonlocal and retardation effects, 
for which the above procedure is not enough. 

The hydrodynamic model described in Sec. A. 2 incor- 
porates both spatial-dispersion and retardation effects. 
When applied to spherical particles, it still yields analyt- 
ical expressions for the resulting nonlocal Mie scattering 
coefficients, which we introduce in Appendix C. Using 
these coefficients in combination with Eq. (33), we pre- 
dict blue shifts in the plasmons of the Al particles of Fig. 
14(a) (dashed curves), even for relatively large spheres. 
This type of dispersion operates on a length scale com- 
parable to the wavelength of conduction electrons at the 
Fermi level, which is of the order of 1 nm for typical 
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metals. This is consistent with the large nonlocal effect 
observed for the quadrupole as compared to the dipole in 
the a — 5 nm sphere, since the I = 2 mode involves more 
sign changes in the induced surface charge that increase 
the relative importance of short distances. It should be 
noted that spatial dispersion has been found to be par- 
tially compensated by the smooth electron density profile 
at the particle surface (Ruppin, 1978). 

In this regard, nonlocality is important at distances 
vp/co (Larkin and Stockman, 2005; Larkin et al, 
2004) , as derived from the minimum momentum transfer 
in the e-h pair region for a given excitation frequency lu 
(see Fig. 8) . This is the so-called Landau damping, which 
has been phenomenologically accounted for through a 
size-dependent local dielectric function, with the dephas- 
ing rate increased by ~ a/vp (Batson, 1980; Hovel et al, 
1993; Kreibig and Vollmer, 1995). Nonlocal effects in no- 
ble metal particles are noticeable up to a size ^ 10—20 nm 
at optical frequencies, but it can be important in other 
materials for larger sizes ^ 100 nm (e.g., silver halide 
nanocrystals, studied with EELS by Oleshko, 2008). 

There are only scarce data on plasmon dispersion in 
small spheres (see, for example, Ouyang et al, 1992), and 
particularly on the spatial distribution of these effects, so 
there is a clear need for further EELS experiments that 
can shed more light into this subjcxrt, which is becoming 
increasingly important in many nanophotonics develop- 
ments that make use of metallic nanoparticles in the sub- 
20 nm range. We hope that these results, together with 
the dimer calculations of Sec. V.E, contribute to trigger 
future research in this direction. 



study was later extended to include retardation and to 
cope with more than two cylinders (Garcia de Abajo and 
Howie, 2002; Manjavacas and Garcia de Abajo, 2008). 

A great deal of work has been devoted in recent years 
to the interaction between metallic nanoparticles in the 
context of nanophotonics, due in part to the strong en- 
hancement of the light intensity in the inter-particle re- 
gion at close proximity (Xu et al., 1999). This has di- 
rect application to biosensing using surface-enhanced Ra- 
man scattering (Talley (d al., 2005). Furthermore, the 
plasmon modes of neighboring particles can be tuned to 
achieve on-demand frequencies by varying their separa- 
tion, as shown in optical measurements performed on 
lithographically printed metal dots (Atay et al., 2004). 
However, there are still many issues that remain unclear 
in the interaction between particles, specially when the 
gap distance is of the order of a few nanometers. Among 
them, the role of Landau damping and other dispersion 
effects. Energy loss spectroscopy turns out again to be 
ideally suited for exploring these situations, in which 
nanometer spatial resolution is needed to resolve optical 
excitations. For example, the singular transition in the 
response of touching and non-touching particle dimers 
(Romero (d al., 2006) has already prompted EELS mea- 
surements (Chu et al, 2009; N'Gom et al, 2008), al- 
though a comprehensive experimental study of this tran- 
sition is still pending. 



3. Coupled nanoparticles 

Coupling between plasmon modes in neighboring 
nanoparticles has attracted considerable interest in the 
electron microscopy community. Early energy-loss ex- 
periments focused on the interaction between neighbor- 
ing spheres (Batson, 1982b, 1985), leading to hybridiza- 
tion of individual particle plasmons, the energies of which 
were successfully explained through analytical deriva- 
tions. Further theoretical work addressed these modes in 
more detail (Ruppin, 1982; Schmeits and Dambly, 1991), 
and eventually the coupling of the electron to two neigh- 
boring spheres was derived in the non-retarded approxi- 
mation (Zabala et al, 1997). 

The interaction between more than two non-aligned 
particles is difficult to explore by analytical techniques. 
A particularly suitable method that yields highly conver- 
gent results, including retardation, consists in expanding 
the electromagnetic field in terms of multipoles around 
the particles, and then allowing these multipoles to prop- 
agate among particles, thus generating a self-consistent 
field (Garcia de Abajo, 1999b). More details on this 
method are given in Sec. III.E.4. 

Cylinder pairs have been theoretically addressed as 
well in the non-retarded limit by Schmeits (1989). This 



E. More complex shapes 

The real world is generally more complex than planar, 
cylindrical, or spherical surfaces. There are numerous ex- 
perimental studies of EELS involving more complicated 
geometries, like Si tips (Reed et al, 1999), metallic parti- 
cles supported by dielectric (Wang and Cowley, 1987a,d) 
or by metal (Wang and Cowley, 1987b, d) substrates, 
MgO wedges (Aizpurua et al, 1999; Garcia de Abajo 
and Aizpurua, 1997), M0S2 platelets (Cohen et al, 2003), 
2D photonic crystals (Garcia de Abajo et al, 2003), Ag 
spheroids (Little et al, 1984), Ag nanotriangles (Nelayah 
et al, 2007a, b), Ag nanoshells (Eggeman et al, 2007), 
and gold nanoparticles ranging from single nanospheres 
and nanorods (Bosman et al, 2007; Schaffer et al, 2009) 
to touching dimers and interacting nanorods (Chu et al, 
2009; N'Gom et al, 2008). The corresponding theoreti- 
cal efforts needed to explain these measurements involve 
more elaborate analytical developments as well as numer- 
ical simulations based upon either Poisson's equation in 
the non-retarded approximation or Maxwell's equations 
if full retardation effects are necessary. These methods 
use in general local-response functions and assume sharp 
boundaries separating different dielectric phases. 
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1. Analytical methods 

Several analj^tical studies of EELS have been carried 
out in the non-retarded approximation for geometries 
that can be naturally described using coordinate systems 
in which the Poisson equation is separable. 



a. Ellipsoids The electrostatic modes of ellipsoids have 
been known for a long time, including analytical expres- 
sions for the polarizability (Gans, 1912; Jones, 1945). In 
particular, the EELS probability has been worked out for 
oblate ellipsoids by lUman et al. (1988). 



h. Wedges Analytical solutions have been reported for 
the modes of sharp wedges (Davis, 1976; Dobrzynski 
and Maradudin, 1972), parabolic wedges (Eguiluz and 
Maradudin, 1976), and hyperbolic wedges (Davis, 1976). 
The latter display a discrete spectrum that suffers a sin- 
gular transition in the limit to a sharp edge, character- 
ized by a continuous spectrum. The EELS probability 
has been calculated for parabolic wedges (Garcfa-Molina 
et al., 1985b) with the beam running parallel to the edge. 



c. Supported particles A sphere or a hemisphere in- 
tersected by a planar surface can be described by 
mode matching, consisting in expanding the potential 
into spherical harmonics and plane waves within dif- 
ferent regions of the structure, thus generating a set 
of self-consistent linear equations upon application of 
the boundary conditions. This method, pioneered by 
Zaremba (1985) to explore the van der Waals interac- 
tion of an atom with a surface feature, was later applied 
by Wang and Cowley (1987a) to study plasmon modes of 
supported particles, which they measured in Al spheres 
sitting on AIF3 (Wang and Cowley, 1987a) and Al sur- 
faces (Wang and Cowley, 1987b). These authors incorpo- 
rated nonlocal effects within the hydrodynamic approx- 
imation (Wang and Cowley, 1987c). Further work on 
this geometry focused on calculating the loss probability 
(Aizpurua et ai, 1996; Rivacoba et al., 1992). 

Several other geometries in which Poisson's equation 
is separable are still unresolved, but the list of structures 
that can Ik^ accessed with analytical tools is short and we 
need in general to rely on numerical methods for more 
complicated shapes. We briefly discuss two of them in 
this section, and a third one is presented in Sec. IILG. 



the advantage that it requires a modest numerical effort, 
since it involves parametrization of boundaries rather 
than volumes. 

The BEM can be traced back to Maxwell, who utilized 
it to calculate capacitances (Maxwell, 1891). It was pio- 
neered by Fuchs (1975) in the computer age to simulate 
modes of ionic cubes, and later by Lu and Maradudin 
(1990) to find plasmons propagating inside channels cut 
into otherwise planar surfaces. Its debut in EELS was 
made by Ouyang and Isaacson (1989b), who obtained 
loss probabilities that they applied to study a particle 
supported on a substrate (Ouyang and Isaacson, 1989a). 
Subsequent developments focused on geometries charac- 
terized by either axial symmetry or translational invari- 
ance along a certain direction, thus permitting us to ac- 
cess a large number of complex morphologies with moder- 
ate computational effort (Garcia dc Abajo and Aizpurua, 
1997), since only curves rather than surfaces need to be 
parametrized in those cases. This approach gives excel- 
lent results as compared to experiment for the losses of 
electrons passing near Si tips (Reed et al, 1999) or close 
to the edge of an MgO cube with different orientations 
of the trajectory (Aizpurua et al., 1999). 

In the non-retarded BEM, the potential inside a given 
homogeneous region of space is separated into two terms, 

(f> = </-'=^* + ?!)"°'', 

the first of which describes the potential produced by 
the external charges present in that medium, treating it 
as if it were infinitely extended, while the second term 
accounts for reflection at the boundary. Assuming an 
inhomogeiieous structure composed of two materials with 
dielectric functions ei and €2, respectively, the interface 
S can be assimilated to an induced surface charge a, so 
that the reflected potential becomes 

r^(r) = [ds 

Js |r - s| 

Inserting these expressions into Poisson's equation, one 
finds the following self-consistent relation for the induced 
boundary charge (Garcia de Abajo and Aizpurua, 1997): 

A(t(s) = ns • V(/)''^*(s) + j ds' F{s, s') a{s'), (38) 

where 



2. Non-retarded boundary element method 

The boundary element method (BEM) consists in ex- 
pressing the fields inside each homogeneous region of a 
given structure in terms of interface sources that are then 
self-consistently calculated by imposing the continuity of 
the potential and the electric displacement. This has 



A = 27r(e2-hei)/(e2-ei), 

and Hg is the interface normal directed towards medium 
2 (see Fig. 15). The right-hand side of Eq. (38) is actu- 
ally the opposite of the electric field along the boundary 

normal, with the integral describing the field produced 
by induced interface charges. Equation (38) admits a 
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FIG. 15 (Color online) Scheme of parametrization used in the 
boundary element method. 



representation in terms of the eigenvectors of F , with 
eigenvalues denoted as 27rAj. This leads to the follow- 
ing general expression for the loss probability (Garcia de 
Abajo and Aizpurua, 1997): 



(39) 



(1 + Aj-)ei + (1 - Aj)e2 



where i runs over the materials 1 and 2, and Aij are 
coupling constants that depend on the actual sample ge- 
ometry and electron trajectory but not on t\ or 62. The 
second term inside the square brackets of Eq. (39) gives 
a negative loss probability corresponding to the begren- 
zungseffekt, the weight of which increases with Aij in the 
same proportion as the interface mode described by the 
first term. 



3. Retarded boundary element method 

We can follow a similar procedure to incorporate retar- 
dation by solving Maxwell's equations at the expense of 
a 6-fold increase in computational cost (Garcia de Abajo 
and Howie, 1998, 2002). First of all, the electric field can 
be expressed in terms of vector and scalar potentials as 

E = ifcA - v./). 

We find it convenient to work in the Lorentz gauge, in 
which V • A = \kt(\). Then, the potentials become 

0(r) = 0-*(r) + f dsG,{\r - s|) a,(s) (40) 

and 

A(r) = A-'(r)+ / ds G,(|r - s|) h,(s) (41) 

inside each homogeneous region j delimited by a bound- 
ary Sj, where 



exp (ikr^yej) 



is the Green function of Helmholz equation inside that 
medium, (j>°^^ and A'^'^* are the potentials produced by 
external sources in an infinite homogeneous medium de- 
scribed by a dielectric function ej , and aj and hj are aux- 
iliary interface charges and currents (see Fig. 15). The 
boundary conditions of the electromagnetic field are then 
imposed and one is left with a set of self-consistent surface 
integral equations that determine the boundary sources. 
For example, the continuity of the scalar potential leads 
to 



ds' [Gi(s - s')fTi(s') - G2(s - s')ct2(s') 



= 0r(s)-0r(s), 

where is the external potential in medium j, eval- 
uated at positions infinitesimally close to the interface. 
We proceed by discretizing the integrals in these equa- 
tions through a finite number of representative points 
distributed on the interfaces. This leads to a linear set 
of equations that we solve by linear algebra techniques. 
Further details on the method and its generalization to 
an arbitrary number of different media have been given 
by Garcia de Abajo and Howie (2002). 

The BEM is particularly advantageous to simulate 
axially-symmetric samples, for which different azimuthal 
components m having exp(im</3) dependence on az- 
imuthal angle are uncoupled. Then, the sample sur- 
face can be parametrized with a contour line, thus ef- 
fectively reducing the problem to one dimension. Simi- 
larly, only the ID contour defining the geometry needs 
to be parametrized for geometries with translational in- 
variance along a given direction. These two special cases 
have been discussed in detail by Garcia de Abajo and 
Howie (2002). However, the method is general and can 
be applied to arbitrary boundaries, as recently illustrated 
by Myroshnychenko et al. (2008a). 

A recent example of application of the retarded BEM is 
presented in Fig. 16, adapted from Nelayah et al. (2007a), 
in which measured loss spectra for Ag nanoprisms [Fig. 
16(a)] are compared to simulations [Fig. 16(b)] when the 
electron beam is aimed at the corner, the edge, or the cen- 
ter of the triangle. Different plasmon losses are observed 
at each of these positions. Moreover, one can obtain snap 
shots of these plasmons, as shown in the images of Fig. 
16(c-h), by scanning the beam over the sample area and 
collecting the rate of selected energy losses corresponding 
to the three dominant features observed in the spectra of 
Fig. 16(a-b). These findings have been corroborated by 
tip-scattering NSOM measurement (Rang et al., 2008). 

The fully retarded BEM has been used to investigate 
the optical response of nanoparticles (Aizpurua et al, 
2003; Grzelczak et al., 2007; Hohenester and Krenn, 
2005; Kumar et al., 2008; Myroshnychenko et al., 2008b; 
Sanchez-Iglesias et al, 2006) and extended nanostruc- 
tures (Blanco and Garcia de Abajo, 2004a; Cole et al., 
2007) , with particular emphasis placed on localized plas- 
mon resonances, the coupling between atoms/molecules 
and nanoparticles (Blanco and Garcia de Abajo, 2004b; 
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FIG. 16 (Color online) Mapping plasmon modes in an indi- 
vidual Ag triangular nanoprism supported on mica (adapted 
from Nelayah et al, 2007a). (a) Loss spectra acquired at 
different positions of the electron beam over the particle (see 
inset), (b) Parameter-free calculated loss spectra, given in 
probability per incident electron per eV of energy-loss range, 
(c-e) Simulated energy-filtered images at the positions of ab- 
solute maxima in the spectra A-C, respectively, (f-h) Mea- 
sured energy-filtered images. 



Triigler and Hohenester, 2008), and the interaction of 
fast electrons with nanoparticles to explain EELS (Chu 
et al, 2009; Garcia de Abajo et al., 2003; Garcia de Abajo 
and Howie, 1998, 2002; Nelayah et al, 2007a; Schaffer 
et al, 2009) and CL (Gomez-Medina et al, 2008; Hof- 
mann et al, 2007; Kuttge et al, 2009; Vesseur et al, 
2008; Yamamoto et al, 2001a) experiments. 



4. Multiple scattering approach 

For specimens containing several neighboring particles, 
multiple scattering techniques provide a natural way of 
solving Maxwell's equations in frequency space. The field 
induced by each of the particles adds up to the scatter- 
ing produced by the other particles in a self-consistent 
fashion. 

In practice, we need to express the electric field around 
each particle j in a complete basis set, for which vector 
spherical harmonics (Wang et al, 1993) are a convenient 
choice. The following specific steps are then followed to 
solve the electromagnetic problem in the presence of a 
passing electron (Garcia de Abajo, 1999b): (1) the exter- 
nal field is scattered by each of the particles, giving rise to 
spherical outgoing waves (single-scattering contribution) ; 
(2) every outgoing multipole centered around a given ob- 
ject is expanded into propagating multipoles centered 
around each of the remaining objects; (3) these prop- 
agating multipoles are scattered by the objects around 
which they are centered, giving rise to outgoing waves 
that add up to the ones obtained after first-order scat- 
tering of the electron field; (4) this process is repeated 
iteratively, giving rise to the self-consistent relation 

= er+^.E^^'4''' (42) 

where and ^j"*^ are vectors containing the coefficients 
of the outgoing-wave expansion of the single-scattering 
contribution and the self-consistent induced field, respec- 
tively, Hjji is the propagation matrix that expresses out- 
going waves centered around object j' in terms of prop- 
agating spherical waves centered around object j, and 
tj is the scattering matrix of object j, containing the 
scattering outgoing-wave coefficients after interaction of 
propagating spherical waves with that object. Analytical 
expressions of tj for homogeneous spheres are given in 
Appendix C. The multipolar expansions are truncated 
at some maximum orbital angular-momentum number 
^max ka for particles of radius a), so that the above 
equations contain L = 2A^(/i„ax + 1)^ variables, where 
iV is the number of particles, and the factor of 2 stands 
for the fact that we have electric and magnetic multi- 
poles, similar to the situation encountered in Eq. (33) 
for the single sphere. Direct inversion of the full sys- 
tem of equations takes a prohibitive computation time 
proportional to L^. Instead, we can utilize a highly con- 
vergent iterative procedure based upon Lanzcos' method 
(see Haydock, 1980) to solve these equations (Garcia de 
Abajo, 1999a, b). These computations demand an effort 
proportional to and yield highly converged spectra for 
electrons passing near clusters of up to several hundred 
particles in a matter of minutes with Imax ~ 15. 

An example of application of this method is presented 
in Fig. 17, taken from Cai et al (2009). A 100 keV elec- 
tron is considered to pass near the end of a linear ar- 
ray formed by 51 spherical gold particles [Fig. 17(a)]. 
The spectral distribution of the electric field [|E'"'^(r, uj)\. 
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plotted in Fig. 17(b) along a segment parallel to the ar- 
ray] exhibits prominent features at the localized plas- 
mon resonances of the spheres. The time-domain pic- 
ture [Fig. 17(c), obtained from Eq. (1)] shows plasmon 
propagation mediated by inter-particle interaction. The 
group velocity (dashed curve) is consistent with the slope 
at the intensity-maximum of the guided-mode disper- 
sion relation derived from the space Fourier transform 
of |E'"d(r,a;)| [Fig. 17(d)]. 



F. Composite materials 

A common situation arises when the sample contains 
composites made of a fine-scale mixture of different di- 
electric grains. Finding the optical response of the result- 
ing metamaterial is not easy, specially in cases for which 
the actual microscopic geometry is unknown. This is even 
more complicated in EELS due to the local nature of the 
electron probe. 

A first naive approach to the problem of explaining loss 
spectra in a composite consists in considering the appro- 
priately weighed bulk probabilities derived from Eq. (19) 
for the constituents. The result is a reasonable approxi- 
mation when the grains are large enough so that interface 
losses are small in front of volume losses. This type of 
analysis has been used to map water distributions in bi- 
ological specimens using valence excitations (Sun et al, 
1993), which are more intense than core losses and pro- 
vide similar information in this case, accompanied by 
lower sample damage. 

However, this approach breaks down for finer mixtures, 
as shown by Howie and Walsh (1991) through studies of 
losses in Al spheres created in AIF3 by the ionizing ac- 
tion of the electron beam. They found that the spectra 
were well reproduced by an expression like Eq. (39) with 
Xj = 1/3 (Howie, 2003; Howie and Walsh, 1991), which 
results from the dipolar mode of the spheres [Eq. (29) for 
1 = 1], now corrected by the bulk losses due to the be- 
grenzungseffekt. Similar conclusions were extracted for 
other systems, like Ni-Al mixed oxides (Cadete Santos 
Aires et al, 1993). 

The Maxwell-Garnett (1904, 1906) effective medium 
theory, which yields good results in many cases, is par- 
ticularly suited to describe small spheres of a given ma- 
terial 2 embedded in a matrix of material 1, and works 
extremely well for small volume fractions of the spheres, 
/ < 0.3. This theory predicts an effective dielectric func- 
tion 



(-1- 



2(1 - /)ei + (1 + 2/)e2 

(2 + /)6i + (l-/)e2 



(43) 



accounting for the spheres through their non-retarded 
dipolar response [af^ in Eq. (B5)]. Not surprisingly, the 
I = 1 mode of Eq. (29) (2ei + 62 = 0) is recovered from 
the zero of the denominator of Eq. (43) in the small / 
limit. The Maxwell-Garnett theory has been successfully 
employed by McComb and Howie (1995) to explain EELS 
measurements in Si02 polymorphs of different density. 



Corrections to the loss function of granular materials 
beyond the dipolar terms captured in Maxwell-Gaxnett 
theory have been explored by several authors (Fu et al., 
1993; Fu and Resca, 1993). In particular, Barrera and 
Fuchs (1995) obtained a nonlocal effective dielectric func- 
tion for a random system of spheres described their local 
bulk permittivity. This work was followed by the deriva- 
tion of a spectral representation of an effective local re- 
sponse for the same system (Fuchs et al, 1996). More 
complicated situations have been considered, such as an 
aloof electron moving above ordered and disordered ar- 
rays of spheres (Mendoza et al, 1998, 1999), and self- 
similar structures formed by a continuous distribution of 
sphere sizes (Sosa et al, 2001), all of them analyzed in 
the non-retarded limit. 

In this context, the remarkable works of McKenzie and 
McPhedran (1977) and McKenzie et al. (1978) provide 
analytical expansions for the effective local response of 
cubic arrangements of spheres with inclusion of an arbi- 
trarily large number of multipoles, which directly yield 
the loss function Im{— l/ceff}- However, a similar deriva- 
tion of EELS formulas for these structures is still pend- 
ing, although they could shed light into the open discus- 
sion of nonlocal effects in metamaterials. 

Retardation corrections to the effective dielectric func- 
tion can be obtained from the photonic bands of periodic 
composites in the long- wavelength limit. This approach 
was followed by Pendry and Martin-Moreno (1994) to 
explain the ~ 8 eV loss feature measured by Howie and 
Walsh (1991) in aggregates of small aluminum parti- 
cles. Their simulations were accompanied by a complex 
series of spectral features that subsequent analysis by 
Garcia de Abajo and Blanco (2003) demonstrated to be 
due to numerical inaccuracies. In fact, the loss spec- 
trum for small filling fraction (/ = 0.06) of an array of 
equally-spaced aluminum spheres already shows a single 
~ 8.5 eV feature (Garcia de Abajo and Blanco, 2003) 
[see Fig. 18(a), solid curves]. This peak splits into two 
for larger / = 0.3, due to multipolar interaction between 
neighboring spheres, quite different from the individual- 
plasmon excitation predicted by Maxwell-Garnett the- 
ory [Fig. 18(b), dashed curves, obtained from Eq. (43)]. 
Interestingly, the calculated EELS peaks in Fig. 18 are 
displaced to the red with respect to the maxima of the 
surface loss function Im{— 1/(1 -|- Ceff)} as a result of 
retardation, similar to the effect discussed in Fig. 11. 
Both the EELS probability and the effective dielectric 
function e^a are calculated here using a layer Korringa- 
Kohn-Rostoker (KKR) method, which allows solving the 
multiple-scattering equation (42) for parallel layers of 
spheres (Garcia de Abajo and Blanco, 2003; Stefanou 
et al, 2000). The aluminum is described by Eq. (Al) 
with Uup = 15 eV and fir] = 1 eV. A value of Zmax = 6 is 
sufficient to achieve convergence. We plot for comparison 
the loss probability computed in the absence of multiple- 
scattering interaction among spheres [Fig. 18(a), broken 
curves] , predicting a single peak that is roughly indepen- 
dent of /. 
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FIG. 17 (Color online) Electron-beam-induced plasmon generation in gold nanoparticle chains (taken from Cai et al, 2009). 
(a) Scheme of the structure and geometry under consideration, (b) Induced electric field |i?'"'^(a;, li;)| (in frequency space lo) 
produced by the electron along the points of the AB segment [see red double arrow in (a), || x]. (c) Same as (b), in the time 
domain: |i?'°'*(a;, f)| [see Eq. (1)]. (d) Fourier transform of E^^'^{x,lo) with respect to x, showing the dispersion relation of 
plasmons excited by the electron and propagating along the chain. The velocity of the launched pulse is found to be ~ 0.75c 
[dashed line in (c)], in good agreement with the group velocity at the maximum of excitation probability along the guided-mode 
dispersion curve [dashed line in (d)]. 
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G. Carbon molecules and low-dimensional structures: 
discrete-dipole approximation 
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FIG. 18 (Color online) (a) Energy loss spectra for an elec- 
tron moving parallel to the [100] direction of the (100) surface 
of a simple-cubic crystal formed by six layers of aluminum 
spheres in vacuum. The lattice period is 5 nm. We consider 
two different filling fractions as shown by labels. The elec- 
tron is moving at a distance of 1 nm from the sphere surfaces 
with velocity v = 0.4c. The dashed curves are obtained by 
neglecting the interaction among spheres, whereas the solid 
curves correspond to the full solution of Maxwell's equations 
using a layer KKR approach, (b) Surface loss function for 
the same crystal as in (a), using the effective dielectric func- 
tion obtained from Maxwell-Garnett theory [broken curves, 
Eq. (43)] and from the layer KKR method (solid curves). 



Carbon molecules and other low-dimensional struc- 
tures have fascinated chemists and physicists since the 
discovery of Cgo (Kroto et al., 1985) and carbon nan- 
otubes (lijima, 1991; Monthioux and Kuznetsov, 2006), 
followed by BN nanotubes (Chopra et ai, 1995) and other 
related molecules. 

Plasmon modes were first observed in solid Cgo as 
6.5 eV and 28 eV EELS features (Gensterblum et al, 
1991; Saito et al, 1991), corresponding to collective os- 
cillations of TT and a electrons, respectively. The optical 
properties of the individual molecules is not significantly 
altered in the solid, since their binding interaction is ba- 
sically due to weak van der Waals attraction. Actually, 
subsequent studies in gas phase Cgo molecules corrobo- 
rated these results (Keller and Coplan, 1992). 

The clear-cut distinction between loosely-bound 7r-like 
orbitals and tightly-bound cr-like orbitals persisted in the 
loss spectra of carbon nanotubes (Kuzuo et ai, 1992), in 
both single-wall and multishcU configurations (Stephan 
et ai, 2002), in multishcU fuUerenes (Henrard et ai, 
1999), in BN nanotubes (Kodak et al, 2000), and in 
other similar types of molecules (Kociak et al., 2001). 

These results were soon related to the dielectric proper- 
ties of graphite, the anisotropic permittivity of which was 
used as input to describe fuUcrcncs and nanotubes via 
continuous dielectric theory, with different e for field com- 
ponents parallel and perpendicular to the carbon shells 
(Lucas et ai, 1994). The dielectric approach yielded rea- 
sonable results compared to measurements of plasmons 
sustained by a diverse number of geometries and chemi- 
cal identities, including multishell BN nanotubes, C nan- 
otubes, and C fuUerenes (Henrard et ai, 1999; Kociak 
et al, 2000; Stephan et al, 2002), WS2 hollow nanoparti- 
cles (Kociak et ai, 2001), and WS2 nanotubes (Taverna 
et ai, 2002). Good agreement was obtained for nan- 
otubes of different thicknesses, down to the single-wall 
level (Stephan et al, 2002). 



26 



Also, first-principles calculations have been carried out 
to describe the interaction of fast electrons with carbon 
nanotubes (Nojeh et ai, 2006). The agreement with ex- 
perimental loss spectra of single-wall tubes is quite sat- 
isfactory (Marinopoulos et ai, 2003). However, this ap- 
proach is currently limited to relatively small structures, 
due to the heavy computational demand of spectroscopy 
ab initio methods. 

The discrete nature of these molecules suggests an 
atomistic approach to describe their response, with 
each atom represented by an induced dipole of differ- 
ent strength in the directions out-of-plane and along 
the tangent plane of the atomic shells. This is the so- 
called discrete-dipole approximation (DDA) , in which the 
atomic polarizability is obtained from the dielectric func- 
tion of the bulk material (e.g., graphite for carbon nan- 
otubes and fuUerenes) by inverting the Clausius-Mossotti 
relation. -'^^ A detailed comparison between the DDA 
and the continuous dielectric theory was carried out by 
Henrard and Lambin (1996), who showed that both ap- 
proaches are in excellent mutual agreement. 

The DDA was pioneered by Purcell and Pennypacker 
(1973) in the study of the optical response of dielectric 
grains by dividing them into fine volume elements that 
were described through induced dipoles. Actually, the 
DDA is a versatile method for dealing with arbitrary sam- 
ple morphologies and is extensively employed to predict 
optical properties of nanoparticles (Draine and Flatau, 
1994; Kelly et ai, 2003; Noguez, 2007). This versatility 
should find application in the simulation of EELS, and 
therefore, we offer a simple derivation of the method that 
generalizes previous approaches to fully include retarda- 
tion effects. 

In order to apply the DDA to carbon structures, each 
of the atoms (or alternatively, each of the bonds in a 
similar formulation) is assimilated to an induced point 
particle of polarizability tensor a^to). The dipoles in- 
duced on these points in response to an external electric 
field E°'^*(r) are then obtained from their self-consistent 
interaction, so that the dipole of the particle located at 
position Tj is given by 



experiment, 
Stephan etsi 
(2002] 
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FIG. 19 (Color online) Probing plasmons in a single-wall car- 
bon nanotube via EELS. The experiment shows prominent 
TT and a plasmons that are relatively well described both by 
the discrete-dipole approximation (DDA) and by a dielectric 
approach. The experimental curve and the dielectric theory 
have been taken from Stephan et al. (2002). The DDA calcu- 
lation performed here for a 16 x 16 nanotube is based upon 
the dielectric function of graphite, taken from Pahk (1991). 



pling with the dipole at r^'.^^ Equation (44) can be 
solved by direct matrix inversion to yield 



1 



P - 



(45) 



where p represents the vector formed by the components 
of all dipoles , the matrix is made of the components 
of Q'jji, and so on. The electric field corresponding to a 
passing electron is given by Eq. (4), that is, E°'^'(rj) = 
(2ecj/t;^7) gj, where gj = g(rj) [see Eq. (5)]. 

Using the same matrix notation as in Eq. (45), the 
fully-retarded loss probability in the presence of our po- 
larizable point particles representing carbon atoms reads 



EELS 



(C) 



1 



1 



(46) 



(44) 



where the 3x3 symmetric tensor Q^^ describes the cou- 



Notice that this expression reduces to Eq. (34) in the 
single-atom case {Q^ — 0). 

As an example of apphcation of Eq. (46), Fig. 19 shows 
the loss probability for an electron passing near a single- 
wall carbon tube under grazing incidence conditions. 



For carbon structures, the effective polarizability per atom is 
obtained from the dielectric function of graphite using = 
(47rnc/3) [{ei + 2)/{ei - I) -Si], where nc = 113.63 nm-^ is the 
atomic density of this material, i =||,-L refers to the response 
to an electric field either parallel or perpendicular to the carbon 
planes, and the constants S± = — 2S|| = 2.8216 originate in 
the lack of inversion symmetry of the graphite atomic lattice 
(Ashcroft and Mermin, 1976). 



The electric field created by dipole p,/ is given by 



X i [(fcu)2 + iku - 1] 



[{kuf + 3iku - 3] 



where u = r j — rjr . This quantity is related to the Green tensor 
G of Eq. (12) in vacuum through g° = -Attlj^G. 
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The experimental data and the results of the anisotropic 
dielectric theory have been taken from Stephan et al. 
(2002). The position of the plasmon features predicted 
by the above DDA formalism compares very well to both 
of them. It is remarkable that dielectric theory is still 
producing meaningful results for such thin nanotubes, al- 
though the DDA seems to attribute a relative weight to 
TT and a plasmons in better agreement with the observed 
spectrum. 



H. Relation to the photonic local density of states 

We can ask ourselves how exactly the loss probability 
is related to the photonic properties of the sample and 
whether it is possible to retrieve from measured data for 
bounded media a relevant quantity similar to the dielec- 
tric function. This question has positive answer in sys- 
tems with translational invariance along the direction of 
elec'tron motion: the EELS probability is proportional to 
the photonic local density of states (LDOS) (Garcia de 
Abajo and Kodak, 2008b)." 

The photonic LDOS is defined by analogy to its elec- 
tronic counterpart in solid state physics as the combined 
local intensity of all normalized eigenmodes of the sys- 
tem (Ashcroft and Mermin, 1976), but unlike electronic 
wave functions, the electric field is a vectorial quantity, 
which suggests defining separate Cartesian projections of 
the LDOS along unit vectors n as (Fussell et al., 2005) 

-2ui 

Pn(r,w) = Im{n- G(r,r,w) • n} (47) 

TT 

in terms of the Green tensor defined in Eq. (12).^^ The 
uniform LDOS in free space is known from black-body 
theory: 

p9(r,u;) = tuV37rV, 



^'^ There is a uiiivocal correspondence between the spatiaUy- 
dependcnt local dielectric function and the Green tensor, as can 
be clearly deduced upon inspection of Eq. (12). This correspon- 
dence works frequency by frequency. Furthermore, the photonic 
LDOS is directly obtained from the Green tensor following Eq. 
(47). The inverse relation is less direct, although one can ex- 
tend the equivalent of the Hohenberg and Kohn (1964) theorem 
to light, thus asserting that there is a univocal correspondence 
between dielectric functions and LDOS, with both quantities de- 
fined within a finite frequency range above a; = (i.e., the cor- 
respondence is not local in frequency). 

The definition of Eq. (47) assumes that the photon eigenmodes 
are well defined (that is, in the absence of absorption) (Colas des 
Francs et al., 2001). An alternative interpretation, which holds 
even in the presence of lossy materials (D'Aguanno et al, 2004), 
comes from the realization that {4n^ujD^ /h) p gives the sponta- 
neous decay rate for an excitation dipole strength D (Barnett 
and Loudon, 1996; Fussell et al., 2004). It should be also noted 
that a complete definition of the LDOS must include a magnetic 
part (Joulain et al., 2003), which is however uncoupled to our 
fast electrons. 



which contains an additional factor of 1/3 arising from 
the projection over a specific Cartesian vector n. 

If we choose the direction of translation invariance 
along the z axis, the Green tensor G{r,r',co) depends 
on z and z' only via z — z', and consequently we can 
write 



Gir,r',u) 



dqz 
27r 



iqz{z-z') 



(48) 



Using this in Eq. (13), the loss probability can be recast 
as 

rEELs(Ro,w) = ^^Im{-G22(Ro,Ro,g'2,w)} 



h 

2TTe^L 

TVjO 



Pz(R-0,9z,w), 



(49) 



where Qz = uj/v is determined by energy conservation 
[see Eq. (3)], L is the length of the trajectory, and 



Pft(R,g;,,w) 



-2uj 



Im{n • G(R, R, (/z, w) • fi} 



is the photonic density of states, local in real space along 
the R directions perpendicular to the beam and local 
in wave vector space qz along the direction of motion, z. 
This relation between the photonic LDOS and the EELS 
probability is further illustrated in Sec. IV.D.l for elec- 
trons moving inside 2D photonic crystals. 

For a non-lossy homogeneous dielectric of real per- 
mittivity e, one finds (R, , oj) = (a'/27rc^)(l — 
q^/k'^e) 9{v'^e — c^), which reflects the fact that only elec- 
trons moving faster than light in the material can lose 
energy [see Eq. (7)]. This is the Cherenkov effect that we 
discuss in Sec. IV.D. As expected, the vacuum density 
of states does not contribute to the EELS signal. 

Equation (49) constitutes a rigorous relation between 
the EELS probability and the LDOS valid for 2D systems 
with translational invariance along the direction of elec- 
tron motion. The LDOS enters this relation through its 
projection over specific components qz of light wavevector 
along that direction. A similar relation can be obtained 
for fully arbitrary geometries from Eq. (13) (Garcia de 
Abajo and Kociak, 2008b): 

4e2 

rEELs(Ro,w) = — Im{-G"z (Ro,Ro,g2,-q'z,w)}, 



27re^ 



Pz(Ro,gz,w), 



where G^^f {11,11' , qz, —q'^,u}) is the Fourier transform of 

the induced Green tensor with respect to z and z' , qz = 
uj/v, and we have defined, by analogy to Eq. (47), the 
quantity 

-2a; 

Pfl(R, q, uj) = Im{n • G(R, R, q, -q, uj) ■ n} 

TT 

as a generalized density of states that is local in real 
space along the R directions and local in momentum 
space along the remaining z direction. In 2D systems, 
we have pn = Lpn. 
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I. Electronic structure determination 

The publications containing EELS studies of local elec- 
tronic structures are so numerous that only a succinct 
review is possible here. A spectacular degree of spa- 
tial resolution is achieved by resorting to high-energy 
losses, for instance when mapping sp^ and sp^ carbon 
in a silicon-diamond interface with sub-nanometer detail 
(MuUer et ai, 1993), and also when imaging K atoms in 
doped fuUerenes (Guan et ai, 2005). However, valence 
losses can provide excellent resolution too with minimum 
sample damage, as illustrated by Daniels et al. (2003), 
who mapped carbon chemical states through character- 
istic bulk plasmon losses with nanometer detail. In gen- 
eral, comparison with ab initio theory greatly facilitates 
the identification of spectral features in the measured loss 
spectra (Keast and Bosman, 2008; Marinopoulos et ai, 
2003). 

The determination of band gaps in insulators is suc- 
cessfully addressed by valence EELS (Dorneich et ai, 
1998), although this requires special care to deal with 
the subtraction of the ZLP and to get rid of CR losses, 
which can be misleading in the assessment of the inelastic 
signal threshold (Gu et al, 2007; Stoger-PoUach, 2008; 
Stoger-PoUach et ai, 2006; Stoger-PoUach et ai, 2008; 
Stoger-PoUach and Schattschneider, 2007). 

The loss probability has been shown to be closely 
related to the so-called joint density of states (JDOS) 
(Pfliiger et ai, 1984, 1985), which is the convolution of 
density of states of occupied and unoccupied electronic 
bands. The JDOS gives an intuitive measure of the 
interband-transition probability when more subtle effects 
in the excitation matrix elements and lifetime effects are 
neglected. The JDOS predicts a dependence of the loss 
probability on the energy transfer ?iaj as ~ {hu; — E^)^^'^ 
and ^ (huj—Eg)'^/'^ at the edge of direct- and indirect-gap 
energies. Eg and Eg, respectively (Rafferty and Brown, 
1998). A good example of this is found in GaN, as shown 
in Fig. 20, adapted from Lazar et al. (2003). A direct gap 
involves a more abrupt decay of the loss probability near 
the gap [Fig. 20(a)-(d)]. Interestingly, both the smooth 
tail of an indirect gap and the sharper edge of a higher- 
energy direct gap are clearly resolved in the loss spectra 
of GaN [Fig. 20(e)]. 

In this context. Van Hove singularities (Van Hove, 
1953) appearing in the electronic density of states at ener- 
gies corresponding to maxima, minima, and saddle points 
of the bands have been also resolved using valence EELS 
(Dorneich et al, 1998). 

Finally, the spatial resolution of electron microscopes 
can be employed at its best to study electronic bands 
of nanostructures. For instance, this has been done to 
measure optical gaps in single-, double-, and triple-wall 
BN nanotubes (Arenal et ai, 2005), showing that the 
gap is rather independent of thickness and stays near 
the in-plane value of hexagonal bulk BN (^ 5.8 eV). In 
a separate development, quantum confinement effects in 
the gaps of individual CdSe quantum dots have been ob- 




FIG. 20 (Color online) Gap determination using EELS, (a-b) 
Direct and indirect gaps, (c-d) Loss probability near direct 
and indirect gaps (adapted from Rafferty and Brown, 1998). 
(e) Example of simultaneous fit of indirect and direct gaps in 
GaN for = 2 eV and = 3.2 eV (adapted from Lazar 
et ai, 2003). 



served using space-resolved EELS (Erni and Browning, 
2007). 

IV. CATHODOLUMINESCENCE: GENERATION OF 
LIGHT BY INCOMING ELECTRONS 

Cathodoluminescence (CL) was first observed in the 
mid nineteenth century as light emission produced when 
cathode rays (electrons) struck the glass of evacuated 
discharge tubes. Since then, a whole industry has been 
generated in the search for outstanding phosphorescent 
materials (Ozawa, 1990), and in particular, CL has been 
extensively used in TV screens, and even TEMs have ben- 
efited from it by projecting sample images in phosphor 
screens. 

Cathodoluminescence offers a method for characteriz- 
ing microstructures causing minimum sample damage. 
It can be regarded as a contactless technique under the 
aloof configuration (see Sec. III.C). CL is widely uti- 
lized in mineralogy (Pennycook, 2008; Pratesi et ai, 
2003; Schieber et ai, 2000) to resolve the composition, 
contamination, and defects of natural structures with 
sub-micron detail. It is also a common tool to study 
the electronic bands of insulators (Balberg and Pankove, 
1971) and semiconductors (Shubina et al., 2004; Yacobi 
and Holt, 1986), and particularly the effect of dopants 
(Holt and Yacobi, 1990). Recently, CL has been used 
as a means to read stored information by tracing phase 
transitions in Ga nanoparticles (Denisyuk et ai, 2008), 
which can be in turn controlled by electron excitation 
(Pochon et ai, 2004). Cathodoluminescence is typically 
performed in adapted SEMs using 1-50 keV beams, and 
it is now possible to record spectra for every pixel of a 
given image (Galloway et al, 2003). In an separate devel- 
opment, low-energy electrons (150-300 eV) emitted from 
scanning tunneling microscope tips have been employed 
to produce surface CL, which benefits from the small in- 
elastic mean free path of electrons at those energies (Ma 
et al, 2006). 
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Here, we are interested in the application of CL to 
investigate the optical properties of nanostructured ma- 
terials. Besides e-h pair recombination, which dominates 
CL in semiconductors, we are interested in other mecha- 
nisms of light emission such as excitation of leaky plas- 
mons in metal nanoparticles. Overall, CL constitutes 
an excellent technique to access optical properties with 
nanometer resolution and without spurious effects due to 
modifications produced by the presence of a probe in the 
region near the sample, as is the case in NSOM. A close 
relative of CL is found in light emission accompanying 
inelastic electron tunneling, which can be assisted by lo- 
calized plasmon excitations (Persson and Baratoff, 1992), 
although the photon yield is very low in that case. 



continuu m of 




FIG. 21 (Color online) Some of the processes taking place 
during energy transfer from a swift electron to a sample, (a) 
Excitation from the specimen ground state g to an excited 
state e, and subsequent radiative decay of the latter (see Sec. 
IV.A.l). (b) Excitation followed by incoherent radiative de- 
cay to an intermediate state i (see Sec. IV. A. 2). (c) Excitation 
followed by non-radiative decay. 



A. Mechanisms of light emission 

We find it convenient to classify the mechanisms of 
CL emission according to their degree of coherence with 
respect to the external electric field of the bombarding 
electrons. 



1. Coherent electron-induced radiation emission 

The interaction of the electron with a material gives 
rise to electromagnetic far-field components that are 
coherent with respect to the external evanescent field 
accompanying the moving charge, described by Eq. 
(4). These far-field components produce what we call 
electron-induced radiation emission (EIRE). There are 
several mechanisms of coherent light emission, the most 
notorious of which are transition radiation (Sec. IV. C), 
Cherenkov radiation (Sec. IV. D), and diffraction radia- 
tion (Sec. IV. E). The electric field associated to each of 
these mechanisms is subject to the same set of Maxwell 
equations, so that we must observe interference effects 
between them. 

The interference of coherently emitted light is actually 
due to the fact that the quantum-mechanical state of 
the sample does not change after the emission process 
takes place. All the energy released by the electron in a 
scattering event giving rise to EIRE is liberated in the 
form of radiation, as schematically shown in Fig. 21(a). 

EIRE sources have a relatively delocalized origin, since 
they involve regions of the sample that are large from 
the electronic point of view, although they can be small 
compared to the wavelength of the emitted light. As 
a consequence, the coherent photon yield increases with 
electron energy: the faster the charge, the more extended 
the external field (see Sec. II. A) and the more efficient the 
coupling to delocalized excitations. 

When heavy encounters with sample atoms take place, 
the sudden change in the electron trajectory gives rise to 
radiation. This is the so-called bremsstrahlung (Jackson, 
1999), which can maintain certain degree of coherence 
with the sources of EIRE discussed so far. In particular. 



theoretical analysis has shown that bremsstrahlung pro- 
duces interference with transition radiation (TR) upon 
electron bombardment (Ritchie et ai, 1964; Shieh and 
Ritchie, 1970). Experimental evidence of this effect has 
been reported under grazing incidence conditions (Cram 
and Arakawa, 1967), so that deflections in the electron 
trajectory occur close to the surface. However, the con- 
tribution of bremsstrahlung to the total photon yield un- 
der usual TEM and SEM operation conditions is minor, 
unless very grazing trajectories are considered (Tomas 
et ai, 1972, 1974). 



2. Incoherent cathodoluminescence 

The incoherent part of the emission bears close rela- 
tion to photoluminescence: the sample is first excited and 
then it decays inelastically by emitting light of a differ- 
ent color, displaced to the red with respect to the exter- 
nal illumination. The only difference lies in the source 
of excitation: while external light preferentially couples 
to strong-dipole modes, electrons can efficiently excite 
metastable levels because they are capable of transferring 
relatively large momenta to the sample (see, for example. 
Fig. 8) that enhances the probability of non-dipole tran- 
sitions. Furthermore, the electron field can be regarded 
as a super-continuum source of evanescent light, suitable 
for creating energetic excitations (for instance, in large 
band-gap insulators) that would be difficult to produce 
using a laser instead. 

Figure 21(b) shows a typical incoherent emission pro- 
cess, whereby the excited state produced by the electron 
decays into an intermediate state, giving rise to light 
of frequency lu' below the transferred uj. The sample 
changes its quantum-mechanical state after the interac- 
tion with the electron and the optical emission take place, 
and consequently, the resulting light cannot produce in- 
terference with coherent EIRE. This contribution to CL 
is similar to the outcome of excited-atom decay (i.e., pro- 
portional to the radiative part of the photonic LDOS). 

Incoherent emission is generally associated to the cre- 
ation of e-h pairs and excitons. Actually, CL is routinely 
used to investigate these types of excitations (Koizumi 
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et al, 2001; Rodt et al, 2005; Saucr et al, 2000), and 
particularly when one is interested in resolving them in 
space. For instance, to study quantum dots (Grundmann 
et al., 1995; Rodt et al., 2005), quantum-confinement in 
semi-conductors (Chen et al., 2006; Yamamoto et al., 
2006a), and optical eS'ects driven by surface atomic steps 
(Minoda and Yamamoto, 2006). 

We should expect the opposite dependence of the in- 
coherent emission yield on electron velocity as compared 
to coherent EIRE. The excitations in this case are highly 
localized, so that the strength of the coupling to the 
more delocalized field of faster electrons turns out to 
be weaker. This effect is actually refiected in the en- 
ergy dependence of the electron inelastic-mean-free-path 
(e-imfp) (Tanuma et al., 1994), represented in Fig. 1: 
this quantity passes through a minimum at ^ 50 eV, in 
which inelastic excitation is more probable, and exhibits 
a monotonic increase with energy due to weaker Coulomb 
interaction when the electron moves faster. Incoherent 
CL can take very large values in suitably designed in- 
sulators [for example, > 5 photons per incident 4keV 
electron in LaFs doped with rare earths (Wang et al, 
2008b)]. 

There are other situations in which the excited states 
involved in incoherent CL are not so localized. For in- 
stance, in emission assisted by nanoparticle plasmons. 
More precisely, d-band holes created by photo-excitation 
in Au nanocrystals have been shown to recombine non- 
radiatively with sp electrons that give rise to excitation 
of intermediate particle plasmons, the decay of which re- 
sults in photoluminescence intensities that are four orders 
of magnitude higher than in metal films (Dulkeith et al, 
2004), but still insignificant compared with electronic de- 
cay channels [see Fig. 21(c)]. Nevertheless, these phe- 
nomena are relevant to understand CL assisted by par- 
ticle plasmons in metal clusters embedded in dielectrics 
(Drachscl et al, 2002; Karah et al, 2005). 

Incoherent CL can be regarded as a way to deliver 
a localized source of light (the inelastic excitations) with 
which one can investigate local photonic properties. This 
is possible thanks to the effect predicted by Purcell 
(1946): an increase in the radiative decay rate of excited 
nuclei coupled to a resonant circuit. The Purcell effect 
is quantified by a linear dependence of the light emission 
rate of excited atoms on the photonic LDOS (see Sec. 
III.H) that has been observed in numerous experiments 
(see Andrew and Barnes, 2001, and references therein). 
But of course, this decay c;lianiiel must compete with non- 
radiative processes and their relative magnitudes depend 
on the dielectric environment. As a thumb rule, atomic 
radiative decay dominates in insulators and is negligible 
in metals. 

A microscopic quantitative analysis of incoherent CL 
is made more complicated by the excitation of hot elec- 
trons, which are in turn a source of CL. Some of these 
electrons escape from the sample surface, giving rise to 
SEE, the yield of which can reach values above 20 elec- 
trons per primary electron in some insulators at ~ 1 keV 



incidence energy (Cazaux, 2006). This high emission 
is due to hot charge carriers right above the insulator 
gap, characterized by long mean free paths and there- 
fore larger probability of escaping from the solid (Ortega 
et al, 1998). 

Secondary electrons of significant energy above the 
Fermi level are less abundant in metals because they are 
rapidly quenched by strong dynamical screening. More- 
over, the CL emission comes from the skin- depth sub- 
surface layer (~ 20 nm in the visible and near-IR for Al 
and noble metals), which drastically erases any effects 
arising from the electron cascade that is produced by the 
energetic impinging probe in deeper regions. 

As noted above, the incoherent CL emission must 
compete with other relaxation channels, like Auger de- 
excitation [see Fig. 21(c)], which are relatively inefficient 
in insulators. However, radiative inelastic decay gives 
rise to just a minor contribution to CL in metals be- 
cause electronic relaxation channels are several orders of 
magnitude faster. Therefore, the CL signal arising from 
metallic nanostructures, such as the ones employed in 
plasmonics, is dominated by EIRE, the details of which 
we discuss in the next sections. 



B. Calculation of coherent light emission 

The intensity of light emission produced by swift elec- 
trons interacting with a sample can be calculated from 
the induced far field. Using the retarded BEM to solve 
Maxwell's equations (Sec. III.E.3), we obtain the far field 
from the integrals of Eqs. (40) and (41) in the kr ^ oo 
limit. The radiation is commonly produced and mea- 
sured in the vacuum medium j = 1, in which ei = 1 and 
Gi(|r — s|) « exp[i(A;r — fcr • s)]/r (for kr :s> 1), with 
f = r/r. Then, the induced field reduces to 

pifer 
fcr— ^oo r 

where denotes the orientation of r, and the far-field 
amplitude f = f — (f • f)f is expressed in terms of the 
boundary current on the vacuum side as 

f{n) = ik [ dse-"=*-"hi(s). (50) 
Jsi 

Here, we are relying on the fact that the contribution 

of — Vc6 to E'"'^ only produces longitudinal components. 
However, the far field is transversal, and consequently, 
entirely coming from the ifcA contribution to E. This 
gives rise to f [Eq. (50)], from which the longitudinal 
part has been subtracted to obtain f . 



A word of caution should be added here; while EIRE is dom- 
inant at relatively-high electron-beam energies (for example, 
above lOkeV), it decays and is overcome by incoherent CL at 
sufficiently-low electron energies. 
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The EIRE intensity can be obtained from the Poynting 
vector integrated over emission directions. The emitted 
energy per incoming electron is then given by 

AE = ^ j dtr^ j f • [E(r,t) x H(r,i)] . 

Now, expressing the fields in time Fourier transform [Eq. 
(1)] and using the far-field limit just discussed, AE can 
be decomposed into photon-energy components Tiuj as 

AE = hujduj j rfnrcL(f^,w), 

where 

is the number of photons emitted per incoming electron, 
per unit of solid angle of emission and per unit of 
photon frequency range u). 



electrons impinging on Ag, Au, Bi, Cu. and W surfaces, 
and showed a dense scries of features that might origi- 
nate in incoherent CL, since coherent EIRE is strongly 
suppressed at such low energies (see Sec. IV.A.l). More 
recently, relativistic beams (6.1MeV) have been utilized 
to produce millimeter TR that is subsequently focussed 
with a parabolic mirror (Naumenko et al, 2008). 

The theory of TR for an electron normally incident on 
a planar surface is covered by several textbooks (Jackson, 
1999; Landau et al, 1984). We outline it here in a form 
that is suitable for obtaining plasmon launching proba- 
bilities in Sec. V.D. We consider a swift electron moving 
along the z axis and crossing the free surface of a metal 
occupying the 2; < region (Fig. 9, right). The electric 
field produced by the moving charge in each medium j 
(with j — I and 2 representing the materials above and 
below the interface, respectively) can be separated into 
the contribution of the external electron field for an infi- 
nite bulk material plus the field reflected at the surface. 



ref 



C. Transition radiation 

When a fast charge impinges on a metal surface, the 
sudden annihilation of its image acts like an induced 
dipole that produces radiation, as illustrated in Fig. 22. 
This efi^ect was predicted by Ginzburg and Frank (1946) 
and observed by Goldsmith and Jelley (1959) in the vis- 
ible range for protons incident on Ag, Al, and Au sur- 
faces. In fact, transition radiation occurs whenever a 
swift electron crosses the interface between different di- 
electric phases. 

A detailed exposition of TR was given by Ter- 
Mikaelian (1972), and a particularly useful result was 
produced by Kroger (1968), who reported a closed- 
form expression for the EELS probability of an elec- 
tron traversing a thin film under arbitrary oblique inci- 
dence, including the contributions of Cherenkov modes 
and TR. This formula has been valuable in devising 
methods to eliminate retardation effects from measured 
spectra and to retrieve bulk dielectric functions (Stoger- 
PoUach, 2008). In a related development, the theoretical 
analysis of the radiation produced by an electron passing 
by the center of a sphere showed an interesting evolution 
from CR (see Sec. I V.D) to TR when considering large or 
small spheres compared to the wavelength (Pogorzelski 
and Ych, 1973). 

On the experimental side, the results of Yamamoto 
et al. (1996a,b) for CR and TR in mica, Ag, Al, and Si 
films, including detailed comparison with theory, consti- 
tute an outstanding reference that illustrates the com- 
plex interplay between both types of radiation emission 
channels in dielectric films as compared to metals (for 
example, Ag and Al show relatively featureless backward 
emission spectra for wavelengths above the bulk plas- 
mons, in contrast to mica or Si). In a separate work, 
Borziak et al. (1976) reported CL spectra for 5-500 eV 



In particular, the electric field inside an infinite medium j 
obtained from Eq. (2) after integrating the z component 
of q reads 



Ef*(r, 



.) = / 



Ef*(q||,2,i^)e'' 



q, R 



with 



\(l»,z,uj) 



47rie 



where q = (q||,w/?;) and <fq = z x q||, so that the 
set {q||,(^q,z} forms a positively oriented 3D reference 
frame. Faraday's law permits writing the corresponding 
magnetic field as 



tij [m,z,uj) = e ' 



q'^ — k'^Cj 



It should be stressed that the bulk fields are evanescent 

away from the trajectory, unless ej is a positive real num- 
ber and v'^y/e] > c^, in which case CR can be produced. 

The reflected component of the field finds its sources 
in the induced surface charges and currents. The charges 
can be in turn expressed in terms of the currents by using 
the continuity equation. After some algebra, we obtain 

E-f(q||,z,a;) = 



The electron velocity varies dramatically inside the metal due to 
both heavy encounters with target atoms and inelastic collisions, 
until it is eventually stopped. However, no relevant changes 
in the trajectory occur with significant probability for energetic 
electrons 50 keV) within a distance of the order of the skin 
depth underneath the surface. 
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FIG. 22 (Color online) Intuitive explanation of transition radiation, (a) An electron cancels its image charge when it crosses 
a metal surface, (b) This cancelation creates an effective dipole capable of emitting light and generating plasmons. (c) The 
resulting transition radiation has indeed a dipolar angular pattern (solid curve) similar to that arising from a dipole normal to 
the surface (dashed curve), and vanishing for grazing emission, in contrast to the pattern of the dipole in the absence of the 
surface (dotted curve) . The angular pattern is calculated for 200 keV electrons impinging on a planar Au surface and for 600 
nm light wavelength. 
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where qzj = y/k-^^j — q\\'^- Here, h^ are the noted surface 
currents (defined on either side of the interface), which 
are determined from the continuity of the parallel com- 
ponents of E and H. Finally, the reflected magnetic field 
reduces to 



Hf (q,|,z,w) = -2^e'«-l^l i?//,sign(z) ^q, (51) 



where 
Ml 

and 



-iuj/v)e2 + qz2ei -iuj/v)ei + qz2<^i 



— k^ei — k^e2 

{uj/v)e2 + g2i£2 _ iuj/v)ei + ^^162 
q^ — k^ei q^ — k'^t2 ' 

2ieg||/c 



D = 



qzie2 + qz2e-i 



(52) 
(53) 

(54) 



Performing the integral over the azimuthal angle of q|| , 
one finds 

Hf (r,c.) = (55) 



-\Sj if 



/ (7||<i(7|| i:>Atje'«-^l^l Ji((7||i?), 
^0 



where (p is the azimuthal unit vector of r. The TR emis- 
sion is obviously p polarized for normal incidence, with 
the magnetic field pointing along the azimuthal direction. 



parallel to the surface. This is in agreement with experi- 
mental observations, in which residual s-polarized emis- 
sion is claimed to originate in surface roughness (Hatten- 
dorff, 1977). 

The homogeneous field H™* decays evanescently away 
from the electron trajectory (see Sec. II. A). Therefore, 
any TR emission must arise from H"^"^, which behaves as 

H-f(r,u;) ^H{e,uj) — 

kr — *oo T 

in the far field. Asymptotic analysis of Eq. (55) shows 
that 

f/f (6*, (jj) = ifc cos 6 D III (p. 

This expression must be evaluated at q\\ = k sin 9, where 
9 is the emission angle with respect to the surface normal. 
Finally, proceeding in a way similar to Sec. IV. B, the TR 
probability is given by the angular integral 



TR 



(co) 



1 

2TThk 



Tr/2 



d9\{{9,uj)f' 



(56) 



performed over the upper hemisphere (vacuum). 

Typical TR spectra are shown in Fig. 23(a) (broken 
curves) for 200 keV electrons incident on Ag and Au sur- 
faces. The spectral emission intensity exhibits an ap- 
proximately linear decrease with wavelength above the 
SP threshold. The integrated emission in the 600-1200 
eV region yields ~ 10~^ photons per incident electron 
and shows a monotonic increase with probe energy [Fig. 
23(b)]. 



D. Cherenkov radiation 

The electric field set up by a swift electron in a ho- 
mogeneous dielectric can exhibit far-field oscillatory be- 
havior according to Eq. (5) whenever the charged probe 
moves faster than light in the medium under the condi- 
tion (7). This is reflected in the fact that the argument 
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FIG. 23 (Color online) (a) Spectral dependence of the TR and SPP emission probabilities upon electron beam bombardment 
under normal incidence. A 200 keV electron impinges on a planar metal surface, giving rise to TR [dashed curves, calculated 
from Eq. (56)] and SPPs [solid curves, obtained from Eq. (71)]. The metal dielectric function has been taken from Johnson 
and Christy (1972). The probability is given per unit of light wavelength range, (b) Dependence of the TR and SPP emission 
probabilities on electron kinetic energy for two different free-space light wavelengths, A = 800 nm and 1500 nm. (c) Measured 
(symbols) and simulated (dashed curve) SPP emission probability spectra for normally-incident 30 keV electrons. The SPP 
emission rate produced by an oscillating dipole placed at the surface is shown for comparison (solid curve), normalized to the 
spectral integral of TR and far-field dipole radiation, (d) Computed TR spectrum for 30 keV electrons normally impinging on 
a gold surface, (c) and (d) have been taken from Kuttge et al. (2009), to which the reader is referred for more details on the 
microscope setup, experimental methods, and theoretical analysis. 



of the Bessel functions in Eq. (5) becomes imaginary un- 
der such conditions. Actually, these oscillations signal 
the emission of light, originally observed by Cherenkov 
(1934) and successfully explained soon after its discov- 
ery (Frank and Tamm, 1937; Tamm, 1939). The theory 
of the Cherenkov effect in homogeneous media was later 
generalized to accommodate frequency dispersion (Fermi, 
1940; Tamm, 1939), and more sophisticated descriptions 
have been produced since then (Lalor and Wolf, 1971; 
Patro, 1982). 

Cherenkov radiation has been extensively discussed 
(see Ginzburg, 1996, and references therein), but a 
particularly interesting demonstration was reported by 
Ohkuma et al. (1991), who observed the effect for 
bunched ultra-relativistic electrons moving faster than 
light in air. 

The probability of CR emission in a homogeneous ma- 
terial was first obtained by Frank and Tamm (1937). 
It can be easily derived by integrating Eq. (20) over 
wavevector transfers, assuming real and local e{uj). One 
finds 



rche(w) 



1 



1 



(57) 



where the dependence on emission frequency comes ex- 
clusively from the dielectric function. Equation (57) is 
a pillar in the design of Cherenkov detectors for high- 
energy physics, in which CR signals the passage of 
charged particles resulting from energetic subatomic col- 
lisions (Julley, 1958). 

Non-trivial finite-trajectory effects have been success- 
fully described using the field of Eq. (4) in structures of 
large dimensions compared to the emission wavelength, 
for which the contribution of reflection at the boundaries 



can be neglected (Ciljak et al, 2003; Ruzicka and Zrelov, 
1993; Zrelov and Ruzicka, 1989). However, this approx- 
imation is not sufficient for small samples like the thin 
films theoretically addressed by Kroger ( 1968) . Actually, 
the experiments reported by Yamamoto et al. (1996a) on 
mica and Si films show a complex evolution of the CR 
intensity with film thickness (particularly for backward 
emission) that is correctly explained by Kroger's formula. 

The Cherenkov effect has been theoretically investi- 
gated for a charge traversing planar interfaces (Kroger, 
1968; Lucas and Kartheuser, 1970; Ter-Mikaelian, 1972; 
Von Festenberg and Kroger, 1968) or moving parallel to a 
surface (Schieber and Schachter, 1998), for layered struc- 
tures (Bolton and Chen, 1995; Lastdrager et al, 2000), 
for an electron moving along the axis of a cylinder (Olsen 
and Kolbenstvedt , 1980), and for beams passing near a 
number of different geometries (Garcia de Abajo et al, 
2004). Other related effects have been explored, like CR 
due to bunched electrons moving along waveguides (On- 
ishchenko et al, 2002), solitons producing CR in optical 
fibers (Akhmediev and Karlsson, 1995; Skryabin et al, 
2003), a second narrower Cherenkov cone in strongly 
dispersive active media (Carusotto et al., 2001), THz 
SPPs emitted by browsing a light spot at superlumi- 
nal velocities on a planar surface (Bakunov et al. , 2005) , 
and inverted Cherenkov cones in materials of negative 
refraction-index (Averkov and Yakovenko, 2005; Matloob 
and Ghaffari, 2004; Veselago, 1968). 

Furthermore, exotic production of CR has been ob- 
served due to superluminal vortices in Josephson junc- 
tions (Goldobin et al., 1998), whereas subluminal CR 
emission has been reported assisted by phonons that sup- 
ply the required extra momentum (Stevens et al, 2001). 
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Homogeneous dielectric 




FIG. 24 (Color online) Tutorial description of the Cherenkov 
effect in a photonic crystal. A charge moving in a homoge- 
neous dielectric with velocity below the Cherenkov threshold 
V < c/n (a) cannot transfer energy and momentum [shaded 
region in (b)] overlapping the light dispersion line [upper line 
in (b)]. However, a periodic modulation of the index of reflec- 
tion (c) allows light to undergo umklapp processes [i.e., ab- 
sorption or emission of lattice momentum, as shown in (d)] , 
so that band folding over the first Brillouin zone of the result- 
ing photonic crystal (e) produces mixing and gap openings in 
the photonic band structure (f). The charge can now couple 
to propagating optical modes. 



1. Cherenkov effect in photonic crystals 

An interesting situation is presented when an electron 
moves inside a photonic crystal, consisting of a periodic 
modulation of the index of refraction along certain(s) 
spatial direction(s). The resulting CR is then affected 
by the photonic band structure, which can considerably 
differ from the vacuum light cone, as schematically il- 
lustrated in Fig. 24 for a ID crystal. The loss of trans- 
lational invariance allows subluminal electrons to couple 
to photonic bands of the crystals, thus resulting in the 
emission of CR. 

This effect has been theoretically analyzed by Luo et al. 
(2003) and Kremers et al. (2009) for electrons moving 
perpendicularly with respect to the pores of a 2D pho- 
tonic crystal. These authors found a complex evolution of 
the electric field, which is molded by the photonic bands. 

An experimental investigation of the interaction of fast 
electrons with a 2D photonic crystal has been reported 
(Garcia de Abajo et al., 2003) and is illustrated in Fig. 
25. The photonic crystal consisted of a porous alumina 
film, grown by electrochemical methods under the appro- 
priate conditions to form a locally regular array of pores 
(Masuda and Fukuda, 1995). The precise control of elec- 
tron optics in a STEM allows focusing the probes so that 
they pass along one of the pores of the crystal in aloof 
configuration, with the beam positioned tens of nanome- 
ters away from the pore wall. The resulting loss spectra 
display prominent peaks around 7-9 eV, depending on 




Energy loss (eV) 

FIG. 25 (Color online) Experimental observation of CR losses 
in a photonic crystal (adapted from Garcia de Abajo et al, 
2003). A 1.4^m thick self-standing porous-alumina film, as 
shown in the TEM image of (a), is used to measure the en- 
ergy loss experienced by electrons moving inside one of the 
holes (b). The measured energy-loss spectrum is represented 
in (c) for 120 keV and 200 keV electrons (symbols) alongside 
theoretical calculations based upon full solution of Maxwell's 
equations for an electron moving inside the 2D photonic crys- 
tal (dashed and solid curves, respectively). The theory has 
been convoluted with a measured ZLP of 0.86 eV FWHM. 

the electron energy, which were attributed to radiative 
losses produced by coupling to photonic modes of the 
crystal. The spectral features are of course related to 
the photonic LDOS, as discussed in Sec. III.H, and they 
exhibit singularities similar to those explained by Van 
Hove (1953), but now arising from characteristic points 
in photonic rather than electronic bands. These results 
are indirect evidence of CR in photonic crystals. A di- 
rect experimental detection of CL in these systems has 
not been yet reported. 

Incidentally, similar effects could be responsible for un- 
usual spectral features in porous silicon at large beam- 
surface separation (Williams et al, 2000). 

E. DifFraction radiation 

We have seen that radiation can be produced when an 
electron crosses a dielectric interface or when it moves 
faster than light in a homogeneous medium. Another 
common source of light emission originates in the passage 
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Energy loss (eV) 

FIG. 26 (Color online) Losses in aloof trajectories due to 
diffraction radiation. The figure shows measured and sim- 
ulated energy-loss spectra for 200 keV electrons passing just 
outside alumina spherical particles. Symbols: measured data 
taken from Abe et al. (2000). Thick curves: retarded EELS 
calculations using Eq. (33). Dotted curve: non-retarded 
EELS calculation for the larger sphere, obtained from Eq. 
(31). Thin solid curves: CL calculations using Eq. (58). The 
computed probability is given per eV of energy-loss range. 
The dielectric function of alumina is taken from Palik (1985). 

of the electron near a structured surface (without actu- 
ally crossing it), giving rise to the so-called diffraction 
radiation (Potylitsyn, 1998). For instance, this effect oc- 
curs if a swift charge passes near and perpendicularly to 
metallic slits, a geometry that has been actually used for 
beam diagnosis in electron accelerators (Karataev et al, 
2004). 

We present in Fig. 26 the effect of diffraction radiation 
on the loss spectrum of electrons grazingly passing with 
respect to alumina spheres of different size. The largest 
sphere produces a sizable inelastic signal in the region 
below ~ 7 eV, in which the material is basically trans- 
parent for the dimensions under consideration. This is 
actually originating in the polarization charges induced 
by the passing electron in the sphere, which oscillate giv- 
ing rise to light emission. Understandably, the smaller 
sphere is poorly polarized, and therefore, couples weakly 
to propagating light, so that it does not produce observ- 
able losses within the low-energy end of the spectrum. 

The emission probability admits a closed-form expres- 
sion similar to Eq. (33) (Garcia de Abajo, 1999c): 

TcL^sphM (58) 

2 oo I / ujb\ 

= ^ E E ^™ ( ^ ) F'- + ^^rn \tf\\ ■ 

i=l m=-l ^ ' ^ 

One can easily see from Eqs. (CI) and (C2) that the 
scattering matrices satisfy the relation \V^\ — ImjtJ'} in 
spheres made of real dielectric function, implying that the 
loss probability equals the CL rate in that case, because 
the material is unable to absorb any energy. The CL rate 
calculated from Eq. (58) coincides with the loss proba- 
bility in the region below 7 eV of Fig. 26 (thin and thick 
continuous curves, respectively), which further supports 
the dominant effect of the diffraction radiation mecha- 



nism in that part of the spectrum. 

Incidentally, the EELS spectra of Fig. 26 show absolute 
maxima in the 15-20 eV range, corresponding to plas- 
mons of the alumina spheres (for comparison, the bulk 
plasmon occurs at 25 eV and the planar SP at 21 eV). 

1. Smith-Purcell emission 

The effect discovered by Smith and Purcell (1953) 
(S-P) is a well-known example of diffraction radiation. 
These authors recorded light emitted when an electron 
beam passes parallel and above a metallic grating. The 
direction and wavelength of emission followed a simple 
relation derived from a Huygens construction. More pre- 
cisely, 

At = dcos9/c^ d/v + Xn/c, (59) 

where n is an integer number, A is the wavelength, and 
the rest of the parameters is defined in the following in- 
set: 




Equation (59) arises from the condition of constructive 
interference in the emission originating in contiguous 
spheres. 

Successive experimental demonstrations of this effect 
were carried out using 30-300 keV electrons to produce 
light of wavelength lying in the fim region (Bachheimer, 
1972; Salisbury, 1970). Also, relativistic beams have been 
employed for generating millimeter radiation (Blackmore 
et al, 2008; Doucas et al, 1992; Shibata et al, 1998), 
which turns out to be very intense along the forward 
direction (Woods et al, 1995). Furthermore, 20-40 keV 
dense electron beams coming from conventional electron 
microscopes and coupled to metallic gratings of periods 
in the sub-millimeter domain have been used to produce 
super-generation of far-infrared radiation (Temkin, 1998; 
Urata et al, 1998). 

Early theoretical analyses of the S-P effect relied on 
diffraction of the evanescent waves accompanying the 
moving electron (Haeberle et al., 1994; Toraldo di Fran- 
cia, 1960; van den Berg, 1973; van den Berg and Tan, 



Within linear response theory, valid for most situations consid- 
ered in this review, the EELS and CL probabilities are propor- 
tional to the square of the particle charge, and they scale linearly 
with the beam intensity for moderate currents. However, bunch- 
ing can produce nonlinear scaling when the electrons move closer 
to each other than the emission wavelength (Urata et al, 1998), 
so that they behave like a composite charge, the square of which 
is obviously larger than the sum of the squares. This leads to 
super-radiant Smith-Purcell emission. 
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FIG. 27 (Color online) Energy loss probability and Smith- 
Purcell radiation emission for an electron moving parallel to 
an infinite periodic string of aligned Al nanoparticles (see up- 
per inset). The probability is normalized per sphere and per 
eV of energy-loss range. The loss probability for an individual 
sphere is also shown for comparison. Adapted from Garcia de 
Abajo (1999a). 



1974). It was also shown that an equivalent descrip- 
tion results from considering the currents induced by the 
passage of the electron on the grating surface (Brownell 
tt al, 1998). Moreover, a large deal of work has been de- 
voted to studying the conditions upon which this effect 
could be utilized for efficient x-ray generation (Moran, 
1992), for instance in free-electron lasers (FELs) (Li 
et al, 2007a,b; Liu et al, 2007). 

The excellent spatial resolution of electron microscopes 
could be used to direct a beam parallel to an array of 
nanoparticles and produce S-P emission. This situation 
is considered in Fig. 27, which predicts intense super- 
continuum emission in the UV domain. This type of sub- 
micron light emitters could find eventual application in 
integrated photonic devices. Interestingly, Fig. 27 shows 
the presence of an absolute threshold below which no 
emission is possible because Eq. (59) cannot be satisfied. 

In recent works, the S-P emission from photonic crys- 
tals has been extensively investigated (Garcia de Abajo 
and Blanco, 2003; Ochiai and Ohtaka, 2004a,b, 2005; 
Pendry and Martm-Moreno, 1994) and compared with 
gratings (Ohtaka and Yamaguti, 2001; Yamaguti et al, 
2002) using theoretical tools. Besides, experimental ob- 
servations have been reported for millimeter waves emit- 
ted by electrons passing near photonic-crystal opals (Ya- 
mamoto et al, 2004). Finally, interesting effects due to 
finite size of the periodic structure along the direction 
of the beam have been noted (Ochiai and Ohtaka, 2006) 
and theoretically analyzed (Horiuchi et al, 2006). 



F. Cathodoluminescence and plasmons 

Plasmons in metallic nanoparticles have finite lifetime 
due to radiative decay. Landau damping (e-h pair cre- 
ation), and other absorption mechanisms. The former is 
rather efficient when the particle has non-negligible size 
compared to the corresponding light wavelength, as we 
discuss below. Actually, the CL signal resulting from 
plasmon decay offers an excellent opportunity to spec- 
trally and spatially resolve plasmons. 



1. Plasmons in metallic films and gratings 

An electron traversing a thin film can give rise to a 
peak in the CL emission at the bulk plasmon frequency 
ujp, as predicted by Ferrell (1958), who found that the 
emitted light was p polarized (magnetic field parallel to 
the film) for normal incidence, and the emission intensity 
vanished at normal and grazing directions for a Drude 
metal. This prediction was soon confirmed by Steinmann 
(1960) and Brown et al (1960) in Ag films, and later by 
Arakawa et al (1964) in Ag and Mg films. A study by 
Ritchie and Eldridge (1962) generalized Ferrell's results 
to arbitrary metals and v ^ c, while the comprehen- 
sive work of Cram and Arakawa (1967) showed excellent 
agreement between theoretical and experimental results 
in Ag films, compared on the same absolute-intensity 
scale. Additionally, the Ferrell mode was observed using 
optical excitation rather than electron beams (Brambring 
and Raether, 1965), and its presence was shown to affect 
the response properties of thin films near the bulk plas- 
mon frequency (Kliewer and Fuchs, 1967). Subsequent 
TEM experiments managed to resolve the dispersion rela- 
tion of the Ferrell mode in Al (Vincent and Silcox, 1973), 
basically following the transversal bulk plasmon curve of 
Fig. 8(b). 

In contrast, SPPs in a thin film cannot directly couple 
to external light because their parallel momentum lies 
outside the light cone, according to Eq. (27), with ei and 
£2 having opposite sign. However, SPPs can be made to 
radiate if they find obstacles in their propagation along 
an otherwise planar surface. In particular, gratings pro- 
vide a source of parallel momentum that plasmons can 
use to couple to external light. The optical emission re- 
sulting from excitation of SPPs by electrons impinging on 
a grating was first observed by Teng and Stern (1967). 
These authors managed to correlate SPP-related peaks 
in the angular distribution of the emission with dips ob- 
served in the specular reflectivity of periodically corru- 
gated surfaces. 

In a remarkable development, Heitmann (1977) used 
CL to map the SP dispersion in sinusoidal Ag gratings 
of different height. The gratings assisted SPPs to couple 
to freely propagating light. His analysis of extrapolated 
results in the limit of shallow corrugation were in ex- 
cellent accord with both the dispersion relation and the 
momentum-energy width predicted by Eq. (27). Further- 
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more, his results for stronger modulation of the grating 
profile resulted in dispersion relations further apart from 
the light line, so that the grating contributed to increase 
the binding of SPPs to the surface. A similar increase 
in surface-mode binding driven by periodic corrugations 
has been also observed in microwaves trapped by drilled 
metal films (Ulrich and Tacke, 1973). This behavior is 
accompanied by gaps in the boundaries of the Brillouin 
zones, as first discussed by Ritchie et al. (1968) for p- 
polarized light incident on gratings, and later studied in 
2D plasmon crystals (Kitson et al., 1996) capable of guid- 
ing plasmons through lines of defects (Bozhevolnyi et al, 
2001). 

Two recent outstanding experiments have indepen- 
dently supplemented Heitmann's work by out-coupling 
electron-generated SPPs in planar surfaces using distant 
gratings (Bashevoy et al, 2006; van Wijngaarden et al, 
2006), thus suggesting the possibility of using fast elec- 
trons as plasmon sources. We address this issue in more 
detail in Sec. V.D. 



2. Plasmon mapping 

Cathodoluminescence has been used to map plasmon 
excitations in scanning mode, by analogy to the EELS 
plasmon-mapping study presented in Fig. 16. The CL 
signal is generally much weaker that the loss intensity, 
which is problematic when good signal-to-noise ratio is 
needed to resolve neighboring spectral features. How- 
ever, CL has the advantage that it can be performed in 
SEMs, which are more widely available than TEMs, and 
it works for thick samples, since there is no need to collect 
transmitted electrons. 

Several examples of application of this technique to 
spatially resolve plasmons in nanostructures have been 
presented by Yamomoto and c;o-workers over the last 
few years, including studies of plasmons in spherical Ag 
nanoparticles (Yamamoto et al, 2001a), gratings (Ya- 
mamoto et al, 2001b), holes in thin films (Degiron et al, 
2004), and particle dimers and nanorods (Yamamoto 
et al, 2006b). 

A characteristic example of these studies is shown in 
Fig. 28 for Ag particles. The spectra of Fig. 28(a) unveil 
the excitation of a broad dipole mode and a narrower 
quadrupole. The measurement is in excellent agreement 
with the theory of Eq. (58). By scanning the electron 
beam over the particle and collecting polarized photons 
emitted at the peak wavelength, one can form the image 
of the dipole plasmon shown in Fig. 28(b). The inten- 
sity of eac;h pixel represents the photon rate measured 
when the electron is positioned at that spot of the sam- 
ple. This plasmon map displays two characteristic max- 
ima on either side of the particle along the direction of 
the polarizer. This study was carried out for particles 
of various sizes, the results of which are summarized in 
Fig. 28(c). The plasmon modes show a distinctive red- 
shift with increasing particle size, in good agreement with 



the theory of Mie (1908), who provided the basis for un- 
derstanding light scattering by spherical particles. The 
origin of these shifts is retardation, as we discuss in Sec. 
IILD.2, and they occur when the diameter of the particle 
is a sizable fraction of the wavelength. As an empirical 
rule, we can assert that retardation effects show up for 
radius a > AZ/27r.^^ These effects are accompanied by 
an increase in the width of the multipole plasmons with 
particle size, which arises as a result of stronger coupling 
to propagating light [see Fig. 28(d)] (Myroshnychenko 
et al, 2008b). 

Figures 28(a)- (b) clearly demonstrate that aloof ge- 
ometries (external trajectories) provide significant infor- 
mation on the particle plasmons. This idea is further 
illustrated in Fig. 29, which presents a quantitative com- 
parison of EELS and CL probabilities for Au disks. The 
difference between these two probabilities corresponds to 
events in which the lost energy is absorbed by the parti- 
cle rather than being re-emitted. These calculations have 
been carried out using the boundary element method de- 
scribed in Sec. in.E.3. 

The CL spectra of Fig. 29 (dashed curves) show just a 
subset of the features observed in the loss spectra (solid 
curves), as an indication that some plasmon excitations 
are dark modes, in the sense that their main decay mech- 
anism is via absorption by the particle rather than cou- 
pling to radiation (Chu et al, 2009). More precisely, 
the feature marked by an arrow in the loss spectrum of 
trajectory A is not found in the light cross section of ran- 
domly oriented disks (not shown), but it can be clearly 
resolved by EELS. 

The spatial resolution of CL should be analogous to 
EELS for mapping plasmons, as shown in Fig. 14(b). 
This figure proves that aloof beams can yield similar 
emission intensities as penetrating electrons, in agree- 
ment with the images of Fig. 28(b). Finally, an addi- 
tional thought on the aloof configuration: the incoherent 
CL signal, which we argued in Sec. IV. A. 2 to be small 
compared to the coherent EIRE contribution reported in 
these figures, should be actually smaller when the elec- 
tron does not traverse the metal and cannot create local- 
ized excitations, which are the main source of incoherent 
CL. 



G. Ultrafast cathodoluminescence 

In a recent development, the CL technique has evolved 
to allow picosecond resolution using a microscope in 



The multipolar order I intervenes in determining the relevant dis- 
tance to which the wavelength has to be compared for assessing 
the relative importance of retardation. This must be of the order 

of the separation between nodes in the induced surface charge, 
which is ^ a/ 1, since the interaction between charges at larger 
distances is partially washed out by cancelation of contributions 
with different sign. 
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FIG. 28 (Color online) Plasmon mapping of single metallic nanoparticles (adapted from Yamamoto et al, 2001a). (a) Measured 
and calculated CL spectra for a 122 nm Ag nanoparticle under grazing incidence, showing dipole {I — 1) and quadrupole {I = 2) 
plasmons. The dielectric function of Ag is taken from Johnson and Christy (1972). (b) Energy-filtered CL scan corresponding 
to the dipole- mode maximum of the spectra shown in (a). The collected light is passing through a polarizer with the orientation 
shown by the arrow, (c) Symbols: measured dispersion of multipole plasmons with particle size. Curves: Mie theory, (d) 
Plasmon width, obtained from Mie theory. 
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count the time of flight delay of the electron in its way 
to the target. This technique has been recently utilized 
to resolve charge-carrier dynamics in quantum dots and 
quantum wells with 50 nm spatial resolution and 10 ps 
time resolution (Merano et ai, 2005; Sonderegger et al., 
2006). Similarly, picosecond electron pulses have been 
used to develop ultrafast electron microscopy (Lobastov 
et al, 2005), so far applied just for simultaneous dynami- 
cal and structural characterization (Barwick et al, 2008; 
Grinolds et ai, 2006). These studies open an avenue 
towards a whole range of interesting phenomena in the 
domain of the ultrafast nanoworld. 



FIG. 29 (Color online) Comparison of EELS and CL inten- 
sities in a gold nanodisk. Two different electron beam tra- 
jectories are considered (see inset) separated 10 nm from the 
disk in both cases. The disk edges have been rounded with a 
curvature radius of 3 nm. These results have been obtained 
using the retarded BEM with the dielectric function of gold 
taken from Johnson and Christy (1972). 



which the electron source consists of a photocathode be- 
ing excited by external light pulses, so that the time- 
resolved CL emission can be assigned an excitation time 
origin determined for each of the pulses, taking into ac- 



V. RELATED PHENOMENA AND SUGGESTED 
EXPERIMENTS 

Electron microscopy is a rather versatile discipline that 
flnds in the nanoworld its natural domain. We can cer- 
tainly think of many new phenomena to be explored with 
focused electron beams, particularly when low-energy, 
valence excitations are involved. In this section, we 
briefly go through some of them. 
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A. Mechanical momentum transfer 

A fast electron beam has the potential to act as a nan- 
otool that manipulates molecules and nanoparticles by 
transferring linear and angular momentum. This should 
become an excellent complement to optical tweezers, cur- 
rently employed to trap small particles ranging in size 
from nanometers to several microns (Meiners and Quake, 
2000), and to study their mechanical behavior by exerting 
electromagnetic forces (Grier, 2003; Korda et ai, 2002). 
Studies based on optical tweezers range from elastic prop- 
erties of deformable living cells (Guck et al., 2000) to 
fluctuation forces in biomolecules (Meiners and Quake, 
2000). Likewise, an electron microscope can modify (and 
simultaneously image) small structures with far greater 
resolution than what is achieved by optical imaging. A 
more sophisticated setup would involve in-vacuo optical 
tweezers integrated in an electron microscope for observa- 
tion of trapped particles. This system is experimentally 
challenging because it requires to trap particles in vac- 
uum, for which pioneering results on optical levitation 
were reported three decades ago by Ashkin and Dziedzic 
(1976). A solution that is perhaps more compatible 
with trapping in vacuum consists in relying on evanes- 
cent fields (e.g., localized SPs), as proposed by Novotny 
et al. (1997) and recently demonstrated by Righini et al. 
(2008). 

Here, we estimate the magnitude of the momentum 
transfer imparted on small objects using the methods 
elaborated in previous sections (Garcia de Abajo, 2004). 
The electromagnetic force is obtained by integrating 
Maxwell's stress tensor over a surface surrounding the 
object (Jackson, 1999). Performing the time- integral of 
this force, we find 



Momentum transfer 



F{t) dt 



¥{Lu)duj, (60) 



where the spectral decomposition in the last identity is 
similar to the one used for the EELS probability in Sec. 
II.B. 

For simplicity, we first limit our discussion to a small 
isotropic particle described by its electric polarizability 
a{Lo). The force exerted by an external field E(r,cij) act- 
ing on such particle is given by the expression (Gordon 
and Ashkin, 1980) 

FH = Re |a ^£;,(r,^) V[ii;,(r,^)]*| , (61) 

where the sum runs over Cartesian components and the 
external field supplied by the electron is given by Eq. (4) . 
Using this in Eq. (61), one finds 



F(^) 
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27lm{a}/ ( — I z 



(62) 
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FIG. 30 (Color online) Momentum transfer from an electron 
to an spherical alumina particle (taken from Garcia de Abajo, 
2004). (a) Sketch of a 200 keV electron passing at a distance 
of 10 nm from the surface of an alumina particle, (b) Lin- 
ear momentum transfer, quantified by the velocity change of 
the particle along directions parallel and perpendicular to the 
beam, and normalized per incident electron. Solid curves; full 
electromagnetic calculation. Dashed curves: dipole approxi- 
mation [Eq. (62), with the prescription of Eq. (36) for a\. 
(c) Change in rotational velocity produced by angular mo- 
mentum transfer from a single electron (full electromagnetic 
simulation) . 



where the function / is the same as in Eq. (35) and the 
particle is located at a distance h from the electron beam. 

Symmetry considerations lead to the conclusion that 
Rayleigh scattering by the particle produces a radiation 
pattern independent of the direction of v. This means 
that the overall transfer of momentum in the form of in- 
duced radiation is zero in the small-particle limit, so that 
Eq. (60) accounts for all momentum transfer emanating 
from the moving electron along the z axis, parallel to 
the trajectory. Then, the contribution of each oj com- 
ponent to the electron energy-loss probability is vFz{uj). 
In fact, we find from Eqs. (62) and (34) the identity 
vFz{u}) — ?iij-'FEELS,dip(w): it is reassuring that the me- 
chanical energy exchanged through momentum transfer 
equals the loss probability multiplied by the transferred 
energy for small targets. 

The mechanical effect can be significant in nanopar- 
ticles, as shown in Fig. 30(b), which compares Eq. (62) 
(dashed curves) with the results obtained when all mul- 
tipoles are included (solid curves). For instance, a typi- 
cal STEM beam with a current of 0.1 nA passing 10 nm 
away from the surface of a 40 nm alumina particle exerts 
a macroscopic attractive acceleration equal to 0.13 times 
the gravity. 

The transfer of angular momentum A£ can be written 



r(w) duj, 



where the torque r is also given by the integral of the 
Maxwell stress tensor (Jackson, 1999) over a surface sur- 
rounding the target. For the small particle considered 
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above, one obtains 



r(a;) 



(63) 
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Interestingly, the expression inside the square brackets 
of this equation vanishes for a non-dissipative particle.^'' 
Equation (63) shows that the torque produced on a small 
particle drops exponentially with impact parameter [see 
Fig. 30(c)], but it is still sufficient to cause an astonishing 
3.9 MHz/s angular acceleration when a 0.1 nA electron 
beam passes 20 nm away from a 40 nm alumina parti- 
cle. The type of experiment needed to demonstrate this 
prediction should be attainable using the combination of 
optical tweezers and TEM suggested above. Perhaps the 
reader shares the curiosity of the author to see what hap- 
pens to the particle as it spins increasingly faster until 
centrifugal forces compete with atomic binding. 



B. Vicinage effects 

At half the speed of light, electrons in a typical 
100 keV, 0.1 nA STEM beam are separated by an aver- 
age distance of 26 cm along the propagation direction and 
they impact the specimen with a time delay of 1.6 ns. The 
interaction among them is therefore negligible in vacuum, 
but a leading electron produces excitations in the sample 
that a second trailing projectile can pick up, whenever 
they last more than the delay time. This situation is 
discussed in Sec. V.B.I. 

Another type of vicinage effect involves the electron 
self-interaction, taking place when its wave function ex- 
tends over different regions of the sample. Current tech- 
nology makes the observation of this effect feasible, and 
it should find direct application to measure nonlocal 
screened interactions, as we show in Sec. V.B.2. 



1. Interaction between two electrons 

This effect bears close relation to the amply discussed 
topic of swift charged-molecule dissociation after travers- 
ing thin foils, in which the interaction among the result- 
ing ionic fragments affects their stopping (Gemmell et al. , 
1975) and the SE yield (Tomita et al, 2006). Likewise, 
mutual interaction between bunched electrons is at the 
heart of super-radiant Smith-Purcell emission (Li et al., 
2006; Urata et al, 1998), as discussed in Sec. IV.E.l. 

We can easily find from Eq. (8) the overall loss prob- 
ability for two electrons moving a distance a apart from 
each other along the same straight-line trajectory. The 



ratio of this probability to that of a single electron is ex- 
actly given by 4cos(a;a/2t;) with independence of the na- 
ture of the sampled excitations. When the electrons move 
in close proximity as compared to u/w, a two- fold increase 
in the combined loss probability is produced. An inter- 
esting situation occurs under the condition a = 'kv/oj: 
the energy released from the leading probe is picked up 
by the trailing one and the overall probability vanishes. 
This may find application to shadow intense energy-loss 
or light-emission features in order to make weaker spec- 
tral details more visible, in the spirit of the masking in- 
terferometry techniques utilized in astronomy (Ireland 
et al, 2008). 



2. Electron self-interaction 

The quantum nature of the electron probe can be ex- 
ploited to retrieve nonlocal information on the dielec- 
tric response if we separate the incident beam into two 
or more coherent spots at the position of the sample. 
This is routinely done in TEM holography by inserting 
a prism in the electron optics, basically exploiting par- 
tial coherence of the energy-loss signal (Herring, 2008; 
Schattschneider and Werner, 2005; Verbeeck et al, 2008). 
Inelastic holography has been performed both in the real 
image plane (Lichte and Freitag, 2000) and in the diffrac- 
tion plane (Herring, 2005), although the current theoret- 
ical understanding of the obtained results suggests that 
further development of the technique is required in order 
to extract new information on the coherence of inelas- 
tic excitations (Verbeeck, 2006). This technique permits 
addressing profound questions on quantum mechanics, 
such as the Feynman thought experiment (Potapov et al, 
2007). 

Here, we discuss a similar scenario, consisting of an 

electron beam split into two coherent beams that are fo- 
cused on identical neighboring spots. This situation is 
well described by Eq. (17) if one assumes that the lateral 
wave function of the incident electron has the form 



V;i±(R) = [F(R - Ri) + F(R - R2) e'^] , 



(64) 



where F gives the spot profile (e.g., a Gaussian of width 
Ai?), Rj are the centers of the two separate spots (|Ri — 
R2I ^ Ai?), and the phase 5 w m^va/fi accounts for 
the difference in path length between the two beams (a), 
which can be controlled through the electron optics of 
the microscope. 

We now make two assumptions regarding the spot size: 
(1) the induced part of the screened interaction W™'^ is 
considered to be a smooth function that does not vary 
significantly over the extension Ai?; and (2) the lateral 



This is a consequence of the optical theorem for small particles, 
stating that Im{ — > 2fe'^/3, with the equal sign only occur- 
ring for non-absorbing particles (van de Hulst, 1981). 



Nonlocality in this paragraph refers to W^(r, r', a;) for r ^ r', and 
not to spatial dispersion in e. 
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variation of the detected plane wave components across 
a distance AR is negligible, that is, pf±AR ^ 1 in Eq. 
(17). The first assumption limits our analysis to excita- 
tions extended over large regions as compared to the size 
of each beam spot. This is quite reasonable for a realistic 
value of Ai? - 2 A, since most valence excitations includ- 
ing plasmons are extended over much larger distances 
(see, for example. Fig. 16). The second assumption de- 
fines the angular region of transmitted electrons in which 
we are interested (for instance, pf± <C 1/AR ~ 0.5 A~^, 
or equivalently, collection angles <C ISmrad for 100 keV 
electrons) . 

Inserting Eq. (64) into Eq. (17), one finds 



dPf± 



cx rgtLs(Ri.^) + rS^Ls(R2,c.) + rf„f (Ri,R2,c.) 



NR./ 



where 



r 

rf„f(Ri,R2,w) = ^ / dzidz2 Im{-VF(ri,r2,w)} 



^ + P/j- ■ (R-i - R-2) — Z2) 



contains the interference between the two spots. 

At variance with the EELS probabihty r^fLg(R,w) 
of Eq. (10), obtained when the entire inelastic signal is 
collected, the quantity dT^^{uj)/dpf± offers direct ac- 
cess to the nonlocal dependence of the screened interac- 
tion, Im{VF(ri, r2, w)}. This requires to have exquisite 
wavevector resolution, Apf± < l/|Ri — R2I, in order to 
prevent cancelations of the interference signal produced 
by the oscillations of the cosine function in Tf^^ . Using 
100 keV electron with an achievable angular resolution of 
font ~ 1 /irad, one should be able to retrieve information 
for spots separated by a distance of ~ 50 nm, which is 
sufficient to study excitations of interest in nanoparticles 
and nanostructured materials. 

We illustrate this formalism by considering two coher- 
ent beams of impact parameters Ri = 5x and R2 = — 6x 
relative to the center of a sphere of radius a, described 
by Eq. (B6). We obtain 



EE 
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(65) 



lm{af^} [1 + (-1)" cos((5 + 2pf^b)] , 



where af^ is defined in Eq. (32). Figure 31 represents 
the results obtained from Eq. (65) for a 200 keV elec- 
tron split into two beams that are passing at a distance 
of 2 nm from a 20 nm gold particle, as shown in the in- 
set. Interestingly, dTf^/dpf± is periodic in pf± with 
period ir/b and is dominated by a dipolar feature . Fur- 
thermore, the interference disappears under the condition 
d + 2pf±b — (2n -|- l)7r/2 for integers n, leading to a loss 
probability proportional to Eq. (31). 




FIG. 31 (Color online) Loss probability according to Eq. (65) 
for two coherent beams passing close to opposite sides of a 
20 nm gold sphere (see inset). The gold is described with 
the dielectric function of Johnson and Christy (1972). The 
electron energy is 200 keV. 



C. Electron energy-gain spectroscopy 

The width of the ZLP limits the resolution of EELS 
and masks low-energy features sitting in the elastic peak 
tail. These problems are partially solved as we show be- 
low if one is capable of reversing the process, so that 
energy gain rather loss is recorded. Electron acceleration 
events were actually reported over 40 years ago when 
Boersch et al. (1966b) observed absorption of thermal 
phonons by electrons transmitted through thin LiF films, 
as illustrated in Fig. 32. The gain signal was later shown 
to depended on temperature as 

r^"M = «MrEELs(^), 

where Feels ('^) is the EELS probability discussed in Sec. 
II and n{uj) = l/[exp{huj/kBT) — 1] is the Bose-Einstein 
distribution function at temperature T (Garcia de Abajo 
and Kodak, 2008a). Actually, Boersch et al. (1966b) 
managed to suppress the measured F^gg^''^ by coohng 
their samples. Unfortunately, thermal transfers are lim- 
ited to modes of energies ^ ksT k, 26 meV at room tem- 
perature. 

A more flexible strategy consists in exciting through 
external illumination only those modes that we want 
to probe, as already suggested by Howie (1999). The 
proposed technique is illustrated in Fig. 33 and can 
be termed electron energy-gain spectroscopy (EEGS) 
(Garcia de Abajo and Kociak, 2008a), in which the elec- 
trons absorb photons from an external light source. How- 
ever, the energy-momentum mismatch between freely 
propagating charges and photons prevents the linear cou- 
pling of fast electrons and light plane waves in vacuum 
(for example, an electron moving in vacuum cannot emit 
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FIG. 32 (Color online) Pioneering evidence of energy gain af- 
ter transmission of fast electrons through a thin foil (adapted 
from Boersch et ah, 1966b). Energy gain events involve 
room-temperature phonon absorption by the electron and is 
strongly reduced at lower temperature. 



CR).^^ Therefore, we must rely on the interaction with a 
material that is structured on subwavelength scales (the 
specimen under scrutiny) to produce evanescent optical 
fields carrying the necessary complex momentum to make 
this coupling feasible. This is already the case in the in- 
verse Cherenkov effect, whereby an electron moving un- 
der the condition given by (7) can absorb energy from a 
laser aligned close to the Cherenkov cone, as proposed by 
More (1966) and observed fifteen years later by Edighof- 
fer et al. (1981). The inverse Cherenkov effect can find 
application to accelerate charged particles (Kimura et al. , 
1995). In a related context, the observation of the inverse 
Smith-Purcell effect (Mizuno et al., 1975, 1987) demon- 
strated the coupling of free electrons to evanescent com- 
ponents of momentum produced by interaction of exter- 
nal light with a grating. 

Here, we consider a more suitable scenario for spec- 
troscopy. Light scattered from a nanostructure, and in 
particular evanescent components in its vicinity, posses 
the necessary momentum to couple to the electron. Con- 
sequently, the energy-gain process can occur near the 
specimen and should yield information on its spectral re- 
sponse, directly mapping the evanescent light strength. 
In such setup, the full width of the gain peak must cor- 
respond to electrons that have absorbed one photon out 
of the external light, and therefore the energy resolu- 
tion of this technique is only limited by the width of 
the external illumination, typically in the sub-meV do- 
main (Rapoport and Khattak, 1988). This is a 100-fold 



In-vacuum optical acceleration of an electron moving along the 
axis of a Bessel beam has been recently proposed (Zhao and Lii, 
2008). 
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FIG. 33 (Color online) Schematic description of electron 
energy-gain spectroscopy (BEGS), (a) An electron passes 
near a nanostructure (a nanoparticle in the illustration) that 
is being illuminated by light of frequency ui. (b) A regular 
EELS spectrum acquired in a scanning transmission electron 
microscope (STEM) shows a gain peak on the left-hand side 
with respect to the zero-loss peak, corresponding to electrons 
that have absorbed one photon out of the external light, (c) 
The electron energy-gain probability for light of frequency 
u is proportional to the area of the energy-gain peak. The 
gain process is assisted by localized excitations of the nanos- 
tructure. They supply the necessary momentum to couple 
evanescent, scattered light to the electron, so that the EEGS 
probability exhibits features as a function of photon frequency 
mimicking the excitation modes of the sample. 



increase in spectral resolution with respect to currently- 
available EELS technology. [Notice that low-energy ex- 
citations (e.g., sub-meV modes) are difficult to resolve 
due to the ZLP background. Here we refer instead to the 
ability to resolve higher-energy excitation features (e.g., 
1 eV) separated by minute energy differences (e.g., 1 eV 
versus 1.001 eV).] 

In order to evaluate EEGS probabilities, we consider 
an energetic electron focused near the sample and de- 
scribe it through the wave function of Eq. (16) (Sec. 
II. C. 2). We adopt the non-recoil approximation and 
assume that the normalized transversal wave function 
ipi± (R) is focused within a narrow spot centered around 
R = Rq, so that only the motion along the beam direc- 
tion z needs to be considered. The electron wave function 
is perturbed by interaction with the external light field 
and picks up components of higher energy to become 



V'ii(z) = 



1 



ipiz -icit 



where we now have Sf = Ei + uj and uj is the photon 
frequency. Using first-order perturbation theory, we find 



- (5'0|| (z) I 



-iCft 



Sijj\\{z)= / dz'Gf{z~ z')H^{z')^P\\{z'), 



(66) 



where H^{z) = {—efi/m(.uj)Ez{^o,z)dz is the energy- 
gain part of the electron-photon coupling Hamiltonian, 
Ez is the z component of the light electric field (parallel 
to the trajectory) , dz denotes differentiation with respect 
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to z, and 



Gf{z - z') 



-I 



dp e'p(^-^') 

27r (me£//7i-pV2 + iO+) 



niePf 



Jpf\z-z'\ 



is the non-rclativistic Green function of the free electron 
propagating along the z direction with final energy TiSf. 
The integral in Eq. (66) must be performed in the z ^ oo 
limit in order to ensure that the electron has already 
abandoned the region of interaction with the nanostruc- 
ture. Then, Sip\\ can be expressed as a coefficient multi- 
plying a plane wave function of energy TiSf. Finally, the 
gain probability is given by the square of that coefficient 
as the dimensionless quantity 



Pi 



EEGS 



{CO) 



J dze-''^'^^E,{Iio,z) 



(67) 



where we have used the non-recoil approximation Pf — 
Pi ~ uj / V Pi for the electron moving with high veloc- 
ity V = Upi/nie- If the incident field is a plane wave and 
the electron moves in vacuum, only scattered components 
contribute to the EEGS probability, because the momen- 
tum gained by the electron {hu/v) is always larger then 
the free- space momentum of light. Therefore, we can 
substitute the induced field for the total field in Eq. (67). 

Useful analytical formulas can be derived for small par- 
ticles, in which the dipolar response is dominant, via their 
electric polarizability a. The EEGS probability is given 
by integration of the dipole scattered field, 

4"^* = a[(wVc^)£^r* + (E°^* • V)5^] exp(iwr/c)/r. 
Direct application of Eq. (67) yields 



-PEEGS,dip(w) = 



2euia 



cob 



(68) 



u)b 



where b is the impact parameter of the electron relative 
to the particle (b _L z) and h = h/b. Interestingly, for 

an external field that is perpendicularly oriented with re- 
spect to the electron trajectory, the scattered field picks 
up z components that couple to the electron, giving rise 
to the contribution in Eq. (68). This expression needs 
to be compared with the EELS probability for a small 
particle [Feels. dip, Eq- (34)]. For a lossless particle, the 
optical theorem (Imja} = {2u}^/3c^) |ap) allows us to 
write TEELS.dip = (pEEGS.dip)/?"/, where p° = u?I'k^(? is 
the free-space local density of photonic states, the factor 
/ = I'K'hiO corrects for the number of photons contained 
in a normalized external field E*"^* (i.e., for lE**^*] = 1), 
and the average over orientations of this field has been 
performed. In other words, the EELS probability, or 



equivalently the CL probability (Fcl) in our lossless par- 
ticle, coincides with the EEGS probability for an incoher- 
ent source containing one photon per photon state. For 
a dissipative particle, we have Feels > Tcl, but the 
relation Fcl = (fEEGS,dip)p°/ still holds. 

We are now ready to make a quantitative comparison 
of EEGS and EELS. Incidentally, Feels has units of in- 
verse energy, whereas Peegs is dimensionless, so that 
the comparison requires integrating Feels over a finite 
energy-loss range ?iAa;. We have -Peegs/(^AwFeels) ~ 
7r|i;^''*|V(?iAwIm{-l/Q;}). For instance, for a 100 nm 
gold nanoparticle illuminated with a laser of moderate 
intensity ~ lO^W/cm^, tuned to the particle plasmon 
energy ?ia; = 2.5 eV, and taking TiAw = 0.1 eV (the plas- 
mon width), we find PEEGs/(^A^rEELs) ~ 40. This 
must be compared with Fcl/I^eels < 1- 

We conclude that EEGS should provide an excellent 
probe of local photonic properties, combining the spatial 
resolution of TEMs and the energy resolution of optical 
sources. The EEGS signal is free from bulk features that 
dominate the loss spectrum as soon as the electron tra- 
jectory intersects the sample. Control over the intensity 
of the pumped light adds an extra handle to facilitate 
the detection of the gain signal, in contrast to CL. Fi- 
nally, EEGS is sensitive to light polarization and it can 
be applied to investigate nonlinear optical response with 
nanometer resolution. 



D. Surface plasmon launching 

As microchip features continue to shrink, 
lithographically-patterned metal structures are be- 
coming natural candidates to replace electronic mi- 
crocircuits. The new devices are expected to operate 
at frequencies above the THz, rather than GHz, and 
to encode information in plasmon signals rather than 
electric currents. The field of plasmonics is actually 
experiencing a tremendous expansion in preparation for 
this emerging technology (Barnes et ai, 2003) but is still 
facing the lack of suitable plasmon sources that provide 
the versatility demanded by highly integrated plasmonic 
circuits, which conventional techniques like frustrated 
total-internal refiection and NSOM cannot satisfy. 

Fast electrons are a promising alternative for generat- 
ing plasmons with exquisite spatial control and relatively 
high efficiency. Plasmon launching by swift electrons was 
already demonstrated several decades ago by Teng and 
Stern (1967) and Heitmann (1977) in metallic gratings. 
These authors utilized the periodic corrugation to cou- 
ple plasmons to propagating light, which was analyzed in 
energy and angle of emission. 

These pioneering studies have been recently supple- 
mented by monitoring SPP generation in metallic planar 
surfaces. The SPPs were eventually decoupled by means 
of a grating placed at some distance from the position 
in which they were generated (Bashevoy et al, 2006; van 
Wijngaarden et al., 2006), thus yielding information on 
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SPP propagation lengths through the dependence of the 
resulting emission intensity on the distance between the 
beam and the grating, and conclusively demonstrating 
the excitation of SPPs by fast electrons. Subsequent 
studies have added more insight into the spectrally- and 
spatially-resolved propagation of SPPs (Bashevoy et al, 
2007; Kuttge et al, 2009; Yamamoto and Suzuki, 2008) 
and the increase in plasmon losses produced by grain 
boundaries in polycrystalline gold compared to single- 
crystal gold (Kuttge et al, 2008). The excitation near 
the plasmon threshold has been also addressed (Peale 
et al, 2008). 

Plasmons can be identified in Eq. (51) with the pole 
in the denominator of D [Eq. (54)], which yields the dis- 
persion relation given by Eq. (27). The q^^ integral of Eq. 
(55) can be then performed in the i? 3> A limit by only 
retaining the contribution of the g'n = q^^ pole. This 
is the so-called plasmon-pole approximation (Ford and 
Weber, 1984), which enables us to approximate 



D 



C 



1\\ 



with 



C = 



2\eq^ 



^SP 



£2 



(69) 



[This comes from the Taylor expansion of £2 qzi + ei qz2 
around the q\\ = q^^ pole.] Now, we proceed by insert- 
ing the plasmon-pole expression back into the reflected 
components of the fields [Eq. (55)] and by approximating 
the integral over using the q^^ = pole of Eq. (69). 
Finally, Eq. (55) reduces to 



H-f(r,u;) « A,-e'«-l^l H['\qfR)ip, (70) 



where 



Aj = wqf sign{z) C1J.J [qf^ , 

and iJ,j is defined in Eqs. (52) and (53). Here, h[^^ is a 
Hankel function, dying off as w (2/7rgy^ii)^/^ exp(ig|^i?) 

at large distance R from the electron trajectory^^ 
(Abramowitz and Stegun, 1972). Equation (70) indicates 
that the induced plasmon field is a cylindrical wave of 

m = symmetry and parallel wavevector q^^ . 

Figure 34(a) shows the real part of the refiected mag- 
netic induction {W^^[r,u>)), calculated from Eqs. (54) 
and (55) as a function of R and cj for r = {R, 0, 10 nm) 



along a silver surface bombarded by a 50 keV electron. 
Surface plasmons produce strong modulations of the field 
at energies below Ritchie's non-retarded limit (hLUg « 3.7 
eV for Ag). The period of these oscillations increases 
with decreasing oj in accordance with the h[^^ behavior 
noted above. The surface field exhibits faster oscillations 
for Lo > LUs, which rapidly decay away from the trajectory. 
The full spatial dependence of B''°' (r, co) is represented in 
Fig. 34(b) for hu = 1.5 eV, clearly showing the excitation 
of a cylindrical SPP wave (see intensity maxima near the 
surface) and TR (oscillatory wave components evolving 
away from the surfac;e). The TR region is dominated by 
oscillations of short wavelength {2n/k) compared to the 
SPP wavelength (27r/Re{g|^}), in consonance with the 

larger momentum of SPPs (Relgiy^} > k). 

We are ultimately interested in determining how many 
plasmons are launched per incident electron, which can 
be calculated from the integral of the Poynting vector 
derived from Eq. (70) over a cylindrical surface centered 
around the electron trajectory. This leads to 



rsp(w) 



E 

i=i,2 



l^j-r J 911 1911 1 



(71) 



for the probability normalized per incident electron and 
per unit of plasmon frequency range. ^'^ 

The prediction of Eq. (71) is represented in Fig. 23(a) 
(solid curves) for Au and Ag surfaces bombarded by 
200 keV electrons, and normalized per unit of emission 
wavelength A [i.e., the plotted quantity is (a;/A)rgp(a;)]. 
Plasmon production is boosted near iVg in both cases, and 
it actually surpasses the production of TR in that region. 
Roughly speaking, ^ lO^'' gold plasmons are produced 
for each normally- incident electron within the 600-1200 
nm spectral range, a figure that reduces to ^ 1.7 x 10~* 
for 20 keV incident energy and monotonically grows with 
electron velocity until it levels up at a value of^2xl0~^ 
around 2 MeV [Fig. 23(b)]. This behavior is qualitatively 
consistent with the spatial extension of the external elec- 
tron field ~ vj/u} [see Eq. (5)], and one should expect 
maximum coupling to the plasmon field when that exten- 
sion is of the order of the SPP wavelength. The depen- 
dence of the emission on primary electron-energy, plotted 
in Fig. 23(b), confirms this explanation. 

Good agreement between the above theory and experi- 
ment has been recently obtained by Kuttge et al. (2009), 
who inferred the absolute SPP yield by first normalizing 
the measured CL signal to the calculated TR in a flat sur- 
face, and then analyzing the interference between SPPs 
decoupled by a grating and TR. The results are shown 
in Fig. 23 (cj. 



The l/y/R dependence of the field is the signature of a surface 
mode, since tlic resulting Poynting vector scales as l/R, so that 
its integral over a circle of radius R centered at the electron 
source stays independent of R, except for the exponential atten- 
uation due to non-vanishing SPP absorption, exp(— 2Im{g?^}i?). 



It should be noted that the approximation of Eq. (69) is fully 
justified for 2> A, that is, in the limit that we have examined 
to derive the plasmon energy flux far away from the electron 
trajectory (Garcia de Abajo, 2007). 
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FIG. 34 (Color online) Spatial visualization of plasmon launching, (a) The magnetic induction B(r, lj) set up by an electron 
impinging on a semi-infinite silver sample is represented as a function of distance to the trajectory along the surface and is 
resolved in energy components hto. Plasmon generation leads to oscillations of increasing period with decreasing plasmon 
energy, in agreement with the dispersion relation of Fig. 11. The induction is calculated at a distance of 10 nm from the surface, 
(b) Full spatial dependence of the Tiuj = 1.5 eV component of the magnetic induction under the same conditions as in (a). 
The left inset shows a log-scale overview stretching 1.2 mm along the surface, in which SPP and TR components are clearly 
separated. The dielectric function of Ag has been taken from Johnson and Christy (1972). 



Propagating plasmons can be generally excited with 
electron beams in any plasmon-supporting system, in- 
cluding metal films, nanowires, and particle arrays (Cai 
et al., 2009). An example of plasmon excitation in a 
linear particle array is discussed in Fig. 17, whereas ex- 
perimental observation of electron coupling to plasmons 
in nanowires was obtained by Yamamoto et al. (2006b) 
and Vesseur et al. (2007). 

The emission mechanism so far considered is coherent 
in the sense explained in Sec. IV. A. 1. However, there 
are at least two other incoherent mechanisms of emission 
involving plasmons: (1) plasmon-assisted decay of local- 
ized excitations produced by the impinging electron close 
to its trajectory and (2) SEs, which can in turn couple to 
the plasmon, acting as primary electrons. These contri- 
butions are only effective within a distance of the order of 
the skin depth inside the metal. We can invoke Bethe's 
theory of charge stopping (Bethe, 1930) to conclude that 
the energy deposited within that region scales as ^ 1/E 
with the primary electron energy E. Consequently, one 
should observe relatively more incoherent emission at low 
electron energies and an increase in the degree of coher- 
ence with faster electrons. 



E. Nonlocal effects in nanostructured metals 

The interaction between closely spaced particles has 
recently attracted much interest because they produce 
large electric field enhancement in the gap region, which 
is beneficial for molecular sensing applications (Xu et al. , 
1999) . Gap distances in the range of one to a few nanome- 
ters have actually been used to measure surface-enhanced 



Raman scattering (SERS) signals, but nonlocal effects 
have been largely neglected in this context until very 
recently (Garcia de Abajo, 2008; Zuloaga et al, 2009). 
Spatially-resolved EELS is well suited to experimentally 
study these unexplored nonlocal effects in dimers, and 
therefore we discuss this problem in this section. 

Spatial dispersion in the material response generally 
produces blue shifts and decreasing intensity of plasmon 
features. Since the early work of Fujimoto and Komaki 
(1968) on spheres, nonlocal effects in EELS have been 
theoretically addressed for several geometries in the non- 
retarded limit, and they have been predicted to show up 
for electron impact parameters in the nanometer range 
(Aizpurua and Rivacoba, 2008; Zabala and Echenique, 
1990). 

Nonlocal effects in the optical response of small metal 
spheres (< 10 nm) have been known for a long time [see, 
for example, Ouyang et al. (1992) for an EELS experi- 
mental study of small silver particles] , and they are well 
described by a phenomenological model incorporating 
conduction-electron scattering at the particle boundary 
as a source of plasmon decoherence (Kreibig and VoUmer, 
1995). This model consists in supplementing the damp- 
ing rate rj in Eq. (Al) with a size-dependent term, vp/a. 
However, the vp/a quantum-confinement model is incon- 
sistent with the fact that an electron state trapped in- 
side a finite cavity (mimicking a conduction electron in a 
nanoparticle) has infinite lifetime. Therefore, many-body 
effects rather than reflection at the boundaries seem to 
be the source of dephasing in small particles, facilitated 
by spatial dispersion, which is responsible for the break- 
down of the local response approximation at distances 
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below a few nanometers (Garcia de Abajo, 2008). 

The hydrodynamic model offers an alternative way of 
accounting for nonlocal effects, as discussed in Appendix. 
A. 2 for a general case and in Appendix C for spheres. 
The results of this model for single aluminum spheres 
are shown in Fig. 14(a) (dashed curves), compared with 
a local description (solid curves). As expected, nonlo- 
cal calculations are blue shifted with respect to the local 
theory. 

For small gold dimers at optical frequencies, it is rea- 
sonable to describe the particles using Eq. (37), and 
to plug the resulting multipolar polarizability into the 
multiple-scattering formalism of Sec. III.E.4. We have 
used Eq. (21) for e{q,u!). The calculated loss probability 
is represented in Fig. 35 for 20 nm particles (see Garcia 
de Abajo, 2008, for more details on the model). The loss 
spectra show a dominant feature that corresponds to the 
lowest-order longitudinal dipole mode of the structure 
(i.e., this mode corresponds to the co-axial dipole-dipole 
mode characterized by a charge pattern " (-1 — ) • (-1 — )" at 
large separations). Significant enhancement of the elec- 
tric field is produced near the gap region at the plasmon 
mode frequency (Fig. 35, lower-left inset), so that the 
details of the dielectric response at small distances be- 
come critical. The local description of the sphere (dashed 
curves) is clearly insufficient at distances below 2nm, 
since they produce redshift and increase of intensity as 
compared to the nonlocal results (solid curves). The re- 
sponse of the dimer is extremely sensitive to the gap dis- 
tance, but electron microscopes are ideally suited to si- 
multaneously determine this parameter and the plasmon 
energy. 
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FIG. 35 (Color online) Loss probability of a 200 keV electron 
passing near dimers formed by spherical gold particles of 20 
nm in diameter. The loss spectra are dominated by the lon- 
gitudinal dipole mode. Local (dashed curves) and nonlocal 
(solid curves) calculations are compared for several separa- 
tions between the particle surfaces, d. Consecutive curves 
are vertically offset by 0.003 eV~^ for clarity. The electron 
trajectory is perpendicular to the dimer axis and located as 
indicated in the upper-right inset. The lower-left inset shows 
the hu) — 2.26 eV component of the total- field intensity in 
logarithmic scale for d — 0.6 nm, calculated from the BEM; 
the field diverges along the electron trajectory (vertical bright 
line) and it has a local maximum in the dimer-gap region. 



VL PROSPECTS FOR PLASMONICS 

Browsing the figures of this review, we observe that 
the historical development of plasmons has been inti- 
mately correlated with advances in electron microscopy 
and spectroscopy. Since the early days of SPs, which 
Ritchie (1957) predicted to show up as characteristic en- 
ergy losses in electrons traversing metallic thin foils (see 
Fig. 10), up until now that particle plasmons can be 
mapped with nanometer precision (see Fig. 16), these 
collective excitations have captured the imagination of 
electron microscopists and material scientists alike. 

Plasmons have now taken off as excellent tools in ul- 
trafast technology for processing information encoded as 
plasmonic signals (Barnes et ai, 2003; Ozbay, 2006; Zia 
et ai, 2006) and for apphcations in biology (Schuster 
et al, 1993) and biomedicine (Burda et ai, 2005). In par- 
ticular, localized plasmon resonances (LPRs) play a rele- 
vant role because they are typically accompanied by large 
enhancement of the electromagnetic field that is useful in 
molecular sensing [for example, using SERS (Xu et al, 



1999)].^^ It is of course important to have experimen- 



Localized plasmons can be confined down to a few nanometers 
in small particles, for which retardation effects are negligible, so 
that LPRs are safely described by Poisson's equation. Actually, 
the minimum size of well-defined LPRs is only limited by broad- 
ening due to nonlocal effects (see Sec. V.E). Now, the question 
arises, are metals (in the optical sense) necessary to produce lo- 
calized electromagnetic resonances? In other words, do we need 
to have negative e materials? A positive answer to this question 
comes from examining the eigenvalues A = 27r(e2 + ei)/ (e2 — ei) 
that appear in the solution of Poisson's equation, recast into Eq. 
(38) (i.e., the solution corresponding to vanishing external field): 
one finds that the eigenvalues are real and they must satisfy 
|A| < 27r, which is only possible if ei and 62 have opposite sign. 
An alternative explanation based upon physical intuition can be 
found in the formal equivalence between Poisson's equation and 
the heat transport equation under stationary conditions, with 
e and (j> playing the role of the thermal conductivity and the 
temperature, respectively: real materials have positive thermal 
conductivity, which guarantees that the second law of thermody- 
namics if fulfilled (heat does fiow towards regions of lower tem- 
perature); but we know that heat cannot be spatially confined 
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tal tools that permit us to visualize LPRs in patterned 
metallic surfaces and nanoparticles. One of these tools is 
NSOM (see Kim and Song, 2007, and references therein), 
which has been successfully applied to map plasmons 
(Bozhevolnyi et al, 2006; Drezet et al, 2007), although 
it is limited to around several tens of nanometers in spa- 
tial resolution. In this context, absorption and scattering 
spectroscopy has been extensively used to characterize 
the optical properties of metallic nanoparticles in colloids 
(Burda et al, 2005), but it requires narrow distributions 
of particle size and morphology, and it can only access the 
average properties of particle assemblies. Similarly, dark- 
field microscopy (DFM) is becoming a popular technique 
in the colloidal chemistry community, allowing localized 
plasmons to be studied in individual nanoparticles (Mock 
et al, 2002). However, unUke EELS, the DFM technique 
cannot provide details on the local distribution of these 
modes. 

Spatially-resolved EELS and CL offer alternative 
methods to characterize plasmons with unmatched res- 
olution. Actually, no other available technique has 
achieved the degree of detail obtained by the plasmon 
studies that we review in Figs. 16, 28, and 36. Fur- 
thermore, the theoretical tools that are necessary for in- 
terpreting EELS and CL spectral images have already 
reached a mature status that allows addressing arbitrar- 
ily complex geometries (see Sec. IILE). 

The plasmon energies and widths obtained from EELS 
and optical extinction are in excellent agreement [Fig. 
36]. In particular, the spectra of dumbbells that we 
show in Fig. 36(c) provide similar information in both 
techniques. However, EELS adds the spatial resolution 
that allows mapping plasmons in real space [see Figs. 
16 and 36(d)]. The resulting EELS images compare 
reasonably well with the near-field enhancement that 
is expected upon conventional optical illumination [Fig. 
36(e)]. Moreover, EELS can access plasmon modes that 
are not resolvable with light (Chu et al., 2009). For ex- 
ample, multipoles of order ? > 1 in small metal spheres 
[Fig. 36(f)]. 

Like far-field optical spectroscopy, the EELS and CL 
techniques benefit from the use of relatively non-invasive 
probes (the electrons) that produce negligible changes 
in the optical properties of what is being observed. In 
contrast, the NSOM tip exerts a strong perturbation on 
the sample. The spatial resolution of EELS and CL is 
now pushed down to the sub-A level, but even a mod- 
est SEM setup equipped with CL can focus an electron 
beam down to < 2 nm spots. The photon-energy resolu- 
tion of CL is similar to that of other optical techniques, 
although the signal-to-noise ratio becomes annoying, par- 
ticularly when the electron does not penetrate the speci- 
men but simply interacts through electromagnetic fields 



for ever using passive materials, so that e > media arc unable 
to trap the nonuniform ip associated to a localized resonance. 



at a short distance from it (aloof geometry). Regarding 
EELS, the 0.1 eV resolution achieved by forefront com- 
mercial STEMs is already sufficient for most plasmonic 
structures, and this situation should improve in coming 
years with further advances in electron monochromators 
and spectrometers. 

Finally, it would be unfair to conclude this section 
without a summary of the impressive amount of infor- 
mation that has been gathered on both localized and ex- 
tended plasmons using electron microscopes. We cannot 
review the vast literature here, but we present a succinct 
overview with some relevant references in Table I. 



VII. CONCLUSION 

Overall, valence excitations, probed by electron mi- 
croscopes, have been a continuous source of information 
on electronic structure and optical response of materi- 
als on the nanometer scale over the last six decades. 
The energetic electrons employed in these microscopes 
permit obtaining structural and spectroscopic details of 
the investigated samples with sub-nanometer resolution, 
and they can be regarded as relatively contactless probes 
compared to other techniques such as STM of NSOM, 
particularly when the aloof configuration is employed, 
in which the electron beam passes outside the targeted 
materials and interacts with them via electromagnetic 
coupling. Optical excitations are the most intense in the 
energy-loss spectra, so that they can be collected using 
low beam-currents, and therefore, sample damage is re- 
duced to small levels, even for penetrating trajectories. 
However, electron microscopes are versatile instruments 
in which the beam current can reach doses that cause 
substantial changes in the specimen, including reshap- 
ing of the morphology of metallic nanoparticles (Batson, 
2008). 

In the context of electron microscopy. CL constitutes 
an inexpensive addition to electron microscopes that can 
be very valuable when combined with appropriate the- 
oretical tools, the versatility of which we have tried to 
illustrate in this review. Electron microscopes in scan- 
ning mode of operation, combined with either EELS or 
CL analysis, are suitable tools for studying the optical 
response of nanostructures, and no other existing tech- 
nique can match them in combined spac;e-eiiergy reso- 
lution. Furthermore, these instruments should facilitate 
the eventual application of swift electrons to light and 
plasmon sources for nanoscale devices. Finally, there is 
a huge background of both experimental and theoreti- 
cal expertise that constitutes a valuable guide for future 
research making use of electron microscopes in nanopho- 
tonics studies. 
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FIG. 36 (Color online) Comparison between EELS, CL, and optical extinction spectroscopy, (a) Electron energy-loss spectrum 
(EELS) of the dumbbell gold nanoparticle shown in the dark-field TEM-image inset, (b) Plasmon maps obtained from the 
intensity of electrons that have lost an energy corresponding to the dipolar longitudinal mode (1.7 eV) and the dipolar transversal 
mode (2.4 eV). (c) Calculated EELS probability (left scale) and extinction cross section (right scale) of a nanoparticle mimicking 
the dumbbell of (a). The cross section is obtained for external electric field orientations both parallel and perpendicular to 
the particle axis, in order to excite longitudinal {!) and transversal (t) plasmons, respectively. The EELS spectra correspond 
to positions of the electron beam as shown in (d) by white circles. The longitudinal mode is red-shifted in the experiment, as 
compared to the calculation, presumably due to the effect of a mica substrate, which is not accounted for here, (d) Calculated 
EELS plasmon maps, with the particle contour shown as a shaded area. Notice that the experimental images of (b) exhibit 
depletion of electron intensity for trajectories passing through the particle, due to elastic scattering of the electrons with gold 
atoms, (e) Near-field plots for longitudinal and transversal polarization of the incident electric field under external optical 
illumination, (f) Comparison of EELS and CL probabilities (solid and dashed curves, respectively) with the optical extinction 
cross section (Text (dash-dotted curve) for a metal sphere described by the dielectric function of Eq. (Al). The damping rate is 
77 = 0.04ajp (similar to Al), the sphere radius is a tenth of the free-space wavelength at frequency uip, and the impact parameter 
is 6 = 1.1a. The experimental data of (a) and (b) are taken from Bosman et al. (2007). The calculations of (c)-(e) are adapted 
from Myroshnychenko et al. (2008b). 
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The local response approximation yields excellent re- 
sults when comparing theory and experiment both in the 
context of nanophotonics and in numerous EELS and 
CL studies that focus on the role of sample geometry. 
The different materials forming a structure are then de- 
scribed by frequency-dependent dielectric functions, usu- 
ally taken from tabulated optical data (see, for example, 
Palik, 1985). 



Several metals, including some alkali elements and alu- 
minum, are remarkably well described by a free-electron 
gas model, in which the local dielectric function takes the 
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plasmon maps 


LPR 


Al nanospheres 


EELS 


2 


-40 


StockU et al. (1997) 




SP 


AI2O3 nanospheres 


EELS 


2 


-30 


Abe et al. (2000) 


diffraction radiation losses 


SPP 


Ag gratings 


CL 


1 


-4 


Yamamoto et al. (2001b) 


plasmon standing waves 


LPR 


Ag nanospheres 


CL 


2 


-4 


Yamamoto et al. (2001a) 


plasmon maps of Mie modes 


LPR 


WS3 nanotubes 


EELS 


3 


- 30 


Kodak et al. (2001) 




LPR 


C nanotubes 


EELS 


3 


-25 


Stephan et al. (2002) 


TT and a plasmons 


LPR 


Ag rough surfaces 


PEEM 




3 


Kubo et al. (2005) 


femtosecond dynamics 


SPP 


Au planar surfaces 


CL 


1 


-3 


Bashevoy et al. (2006) 
van Wijngaarden et al. (2006) 


SPP propagation 


LPR 


Ag sphere dimer 


CL 


2 


-4 


Yamamoto et al. (2006b) 


plasmon maps 


LPR 


Ag nanowires 


CL 


1 


-4 


Yamamoto et al. (2006b) 


plasmon maps 


SPP 


Au atom wire array 


LEEM 


1 


-3 


Nagao et al. (2006) 


ID plasmon dispersion 


LPR 


Ag triangular nanoprisms 


EELS 


1 


-4 


Nelayah et al. (2007a) 


plasmon maps 


LPR 


Ag nanowires 


CL 


1 


-4 


Vesseur et al. (2007) 


plasmon maps 


SPP 


Au planar surfaces 


CL 


1 


-3 


Bashevoy et al. (2007) 


hyper-spcetral imaging 


SPP 


Ag planar surfaces 


PEEM 




3 


Kubo et al. (2007) 


femtosecond dynamics 


LPR 


Au annular grating 


CL 


1 


-3 


Hofmami et al. (2007) 


plasmon maps 


SP 


Iff02 in stacks 


EELS 


2 


-30 


Couillard et al. (2007) 




SP 


Si/Si02/Si stacks 


EELS 


2 


-30 


Couillard et al. (2008) 




SPP 


Au planar surfaces 


CL 


1 


-3 


Kuttge et al. (2008) 
Peale et al. (2008) 
Yamamoto and Suzuki (2008) 


SPP propagation 


LPR 


Au ridge 


CL 


1 


-3 


Vesseur et al. (2008) 


plasmon maps 


LPR 


Au nauorods, nanospheres 
ellipsoids, dimers, and 
touching dimers 


EELS 


1 


-4 


Bosman et al. (2007) 
N'Gom et al. (2008) 
Chu et al. (2009) 
Schaffer et al. (2009) 


plasmon maps and spectra 



TABLE I Historical overview of experimental surface plasmon studies making use of electron beams. 
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form 



where 



e°(a;) 



u{u) + it]) ' 



(jp = (47re^no/me)^^^ 



(Al) 



(A2) 



is the bulk plasmon frequency, determined by the conduc- 
tion electron density no,^^ and 77 is a small phenomeno- 



^® Equation (Al) describes the correct way of incorporating damp- 
ing, via u!{ijj+iri) rather than (w + i??)^. The latter is inconsistent 
with conservation of the number of electrons in the metal (Mer- 
min, 1970). 

^'^ The dependence of uip on no has been experimentally monitored 
by Imbusch and Niedrig (1970) using EELS and changing no 
through thermal expansion in the 3 — 295 K range. 



logical relaxation rate that accounts for electromagnetic 
damping in the material (e.g., firj = 0.6 eV in Al). This 
model gives hUp = 15.8 eV for Al, in good agreement 
with the value of 15 eV obtained from EELS measure- 
ments (sec Fig. 8 and Tabic II). 

Some ionic materials, and alkali-halide crystals in par- 
ticular, are better described by assuming localized ionic 
oscillation modes, which result in the dielectric function 
(Ashcroft and Mermin, 1976) 



e{(j) = Co 



uj^oj + iri) — ojI 



(A3) 



where is the frequency of transverse vibrations and 
t^ion signals a plasmon-likc longitudinal mode character- 
ized by e « 0. Retardation produces dispersion of these 
modes for momentum transfers comparable to %k, as 
shown in Fig. 5. 

Neglecting the influence of the atomic lattice, the ratio 
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between components of the external potential and the to- 
tal potential (resolved in momentum and frequency) de- 
fines in general a dispersive dielectric function e{q,uj). 
Spatial dispersion becomes important for wavevector 
transfers q of the order of oj/vp, where the Fermi ve- 
locity vp is typically two or three orders of magnitude 
smaller than the speed of light c. Consequently, the de- 
pendence on q is typically important for q ^ k, away 
from the region showing retardation [see Fig. 8(b)]. 



1. Lindhard and Mermin dielectric functions 

Spatial dispersion effects in good metals such as alu- 
minum and silver can be accounted for by assimilating 
their conduction band to a gas of non-interacting free 
electrons (Ashcroft and Mermin, 1976). An analytical 
formula for e{q,u}) in the free-electron gas was derived 
by Lindhard (1954). The Lindhard formula, which pre- 
dicts a continuum of e-h pair excitations and collective 
plasmon modes as illustrated in Fig. 8(a), is obtained 
by perturbing the plane-wave representing each electron 
within the Fermi sphere by means of an external poten- 
tial of the form (j)'^^*^ = exp(iq • r — iojt). The sum of these 
perturbations for all conduction electrons produces an 
induced charge p™'^ proportional to 0'''^*. The system is 
then assumed to be perturbed by the self-consistent total 
potential rather than the external potential in what has 
been termed the random-phase approximation (RPA). 
This leads to (Ashcroft and Mermin, 1976) 



e{q,uj) 



1 + 



47r p" 



1 



ir'^hq'^ 



,3, /k-q/2 - /k-Fq/2 



?ik • q/rrie -|- u) ' 



where the integral extends over 3D electron wavevectors 
k, fi; is the Fermi distribution (approximately 1 for fc < 
kp and otherwise at room temperature), and kp = 
nieVp/h, is the Fermi wavevector. This integral admits 
the closed-form expression (Lindhard, 1954) 



eL(g,a;) = 1 + 



2mee^kp 



[l + R{x,y) + R{x,-y)], 



where x = q/kp, y = ui/ep, hep = h'^kp/2me is the 
Fermi energy, and 



R{x,y) 



2x 



2x 



In 



2a;- 



2a;- 



Here, w is assumed to have a positive imaginary part. 
If Imjo;} is infinitesimally small then lm{e^{q,uj)} takes 
nonzero values only for x(x — 2) < y < x(x + 2) (i.e., in 
the so-called e-h region) , which corresponds to the shaded 
area of Fig. 8(a). 

A naive way of introducing a damping rate rj in the 
electron-gas motion consists in adding a finite imaginary 



part to uj, but this results in a net loss of the number of 
electrons with time. The Lindhard formula was extended 
by Mermin (1970) to include finite damping while main- 
taining the number of electrons constant. He found the 
expression 



,M 



(a,w) 



(A4) 



= 1 + 



jio + iTl) [e^{q,u> + iTj) - 1] 
uj + iri [e^{q, a; + It?) - 1] / [e^{q, 0) - 1] ' 



where 



eL(g,0) = l + 



2me'^kp 



1 + 



4a; 



In 



a;-2 



is the static limit of f^. 
2. IHydrodynamic model 



A simpler approach, which can be extended to nonlocal 
structured metals, consists in describing the conduction 
band as a classical electron plasma (Bloch, 1933). We 
generalize this approach here to incorporate retardation 
effects. The force equation for the hydrodynamic flow 
takes the form^^ 

{d/dt + r])r'^ = (eVme)noE (A5) 
+ {n/mef [P^ - (1/4)V2] Vp'"^ 

where p"^'* and j™'* are the charge density and cur- 
rent density perturbations induced in the electron gas, 
E is the self-consistent electric field, Uq is the unper- 
turbed electron density, and the hydrodynamic pres- 
sure is related to the Fermi wavevector kp through 
the coeflicient (3 = {5/5y/^kp (Ritchie, 1957). We 
have supplemented this equation with the cutoff term 
— (/i/me)^(l/4)V^Vp'"'^, which accounts for e-h pair ex- 
citations in an effective way. The response of the system 
is then obtained by combining Eq. (A5) with Maxwell's 
equations for the total charge and current (p'^'^'-l-p'"'^ and 
jext -|- jind^ respectively) in the presence of the ionic back- 
ground, the response of which is accounted for through a 
permittivity ej,.^^'^'^ This formalism is a direct generaliza- 
tion of the common non-retarded approach that is based 



The linearized Bloch hydrodynamic equations constitute a suit- 
able approximation within linear response under the assumption 

^® Interband transitions in aluminum give the main contribution 
to 77, but €(, = 1 works well for that material, indicating that 
they are not strongly polarized near the bulk plasmon frequency 
(Raether, 1980). In contrast, interband transitions play a domi- 
nant role in noble metals such as Ag, Au, and Cu (Johnson and 
Christy, 1972), where ei, differs considerably from 1 (see Table 
II). 

Notice that the continuity equation V • j = iu>p is implied by 
Maxwell's equations. 
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Material 


Equation 


Tiu)p or /iwion (eV) 


hr) (eV) 


ib or Coo 


hjJtr (eV) 


Energy range (eV) 


Ag (a) 


(A7) 


9.17±0.19 


0.025±0.01 


«4 


- 


< 1.8 


Al (b) 


(Al) 


«15 


«0.6 






8-25 


Au (a) 


(A7) 


9.06±0.18 


0.071±0.007 


^9.5 




< 1.8 


Cu (a) 


(A7) 


8.88±0.18 


0.096±0.01 


«8 




< 1.8 


KCl (c) 


(A3) 


«0.0274 


wO.0014 


«2.19 


wO.0185 


0.01 - 0.025 


SiC (b) 


(A3) 


wO.075 


wO.012 


W6.15 


wO.099 


0.01-0.2 



TABLE II Dielectric parameters of selected materials, described by Eqs. (Al), (A3), and (A7), and chosen to fit the measured 
dielectric functions reported in the following references within the specified light-energy range: (a) Johnson and Christy (1972); 
(b) Palik (1985); (c) Ashcroft and Mermin (1976). The reader is referred to the book of Raether (1980) for comprehensive 
tables of bulk and surface plasmon energies and widths. 



upon Poisson's equation and expresses the current den- 
sity through a velocity potential ^ as j'"'^ = noV^ (see, for 
instance, Ritchie, 1957, for appHcation to planar surfaces, 
Fujimoto and Komaki, 1968, Crowell and Ritchie, 1968, 
and Barberan and BauscUs, 1985, for spheres, and Wang 
and Cowley, 1987c, for supported particles). The cus- 
tomary boundary conditions of the electromagnetic field 
at the metal surface must be now supplemented by the 
condition that the normal c;oinponciit of j'"*^ is zero. The 
results of this formalism for planar surfaces and spheres 
are presented in Appendices B.l and C, respectively. 

In homogeneous media, it is convenient to project the 
hydrodynamic and Maxwell equations in the space of 
wavevector q and frequency w, from which we can work 
out the longitudinal and transversal dielectric functions 
eion and etr using the definition of electric displacement 

D = eionElon + CtrEtr 

in terms of electric field components parallel and perpen- 
dicular to q, respectively. We obtain 



eion(5,w) = £6 



(A6) 



and 



u){u} + vq) 



Eq. (A6) was derived by Ritchie (1957) from an expan- 
sion of Lindhard's function for large uj. The dispersion 
relation of longitudinal and transversal bulk plasmons is 
represented in Fig. 8(b) for aluminum. 



APPENDIX B: Non-retarded screened interaction in simple 
geometries 

The electrostatic screened interaction W{v, r', u) is de- 
fined as the amplitude of the potential created at r by a 



charge of magnitude cxp(— icjf) placed at r'. For mate- 
rials described by a local, frequency-dependent dielectric 
function e(r, w), Poisson's equation leads to 



V • e(r,w)Viy(r,r',w) = -Att5{y - r'). (Bl) 
In free space, one has W°(r — r') = |r — r'|~^. While Eq. 
(Bl) is valid for arbitrary, smooth spatial dependence of 
the permittivity, most situations encountered in practice 
are reliably modeled by sharp boundaries. We can then 
obtain W from the Coulomb interaction in eax;h homo- 
geneous region after imposing the continuity of both the 
potential and the normal displacement at the interfaces. 
Next, we offer analj^ical expressions for W in the simple 
cases of planes, cylinders, and spheres. 



1. Planar surface 

A planar interface at 2; = is best described by ex- 
pressing W in the space of parallel wavevector qy. In 
particular, the bare Coulomb interaction can be written 



(B2) 



where r = (R, z) and 



w^im,z,z') = — c-«iii^-^ I. 

911 



The full screened interaction W admits an expansion sim- 
ilar to Eq. (B2), in which W^{q\\,z,z') must be replaced 
by 
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W{qii,z,z',uj) 



2tt 



otherwise, 



(B3) 



. ei+e2 



using the notation of Fig. 9. 

Nonlocal effects can be incorporated through the hydrodjmamic model outlined in Appendix A. 2. The result in the 
non-retarded approximation reduces to 



27r 



'e-9iil--'l + £iLW_ic-9ii(-+-'), 



£3(0)+! 



^MV)_p-g||z 
e»(0) + l*^ ' 

U3(0) + l^ 



z,z' > 0, 
z, z' < 0, 
z>0, z' <0, 
z <0, z' >0, 



where 



Csiz) = 



dqz 



(p eion(g,t^)' 



and = q\\'^ + q^. An analytical expression for es{z) can be obtained using Eq. (A6) (Garcia de Abajo and Echenique, 
1992). This approach is formally equivalent to the specular reflection model (Ritchie and Marusak, 1966; Wagner, 
1966), which applies to any eion(9, <^)- 



2. Cylinder 

For axially-symmetric problems, it is convenient to write the bare Coulomb interaction as 



where 



-1 -1 /"oo 

^\ = lo:;^ d9.e'«^(^-^') e™('^-'^')t^°(i?,i?',g.,a;), 



W° (i?, i?', q,) = AnI„MzR<) Km{qzR>), 



i?< = min{i?, i?'}, and i?> = max{i?, i?'}. In a homogeneous cylinder centered around the z axis (see Fig. 13), the 
screened interaction is found to be 



Wm{R,R',qz,u))=47rx { 1 

^^Im{qzR<)Km{qzR>), 

where A„ = e2H^{x)K^{x) - eiIm{x)K'^{x) and a; = g^a. 



Im{qzR<) + ^(ei - ^2)Ini(.x)r^{x)K„i{qzR<) K„,{qzR>), if R,R' > a, 
K^{qzR>) + s^(ei - e2)K^{x)K'^{x)I^{qzR>)\ I^{qzR<), if R,R' < a, (B4) 

otherwise, 



3. Sphere 

The bare Coulomb interaction can be expanded in spherical harmonics as 



' Im 



where 
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r< = niin{r, r'}, r> = max{r, r'}, and il denotes the direction of r. Using the notation of Fig. 13, the interaction 
near a homogeneous sphere is given by the above expansion with 

i-W^ir,r') 
Wi{r,r',Lo) = <j iVKf(r,r' 




substituted for W;°(r,r'). Here, rji = {I + l)ei + le2. It 
is also useful to write the r,r' > a part in terms of the 
multipolar polarizability 



and 



Ke2 - ei) 



lei + + l)ei 



as 



W;(r, r'.w) 



<(r,r') 



(rr') 



/V+l 



(B5) 



(B6) 



tf (C2) 
^ -eij;(a:i)[x2i;(3;2)]' + e2[a;ij;(a:i)]'i;(a;2) 
e^h!i^\xi)[x2ji{x2)]' - e2[xxh^i^\xx)]'ji{x2) ' 



APPENDIX C: Retarded scattering matrices and coupling 
coefficients in a sphere including nonlocal effects 

The coupling coefficients of Eqs. (33) and (58) are di- 
mensionless functions oi v/c that admit the following 
closed- form expressions: 



and 



M 



1 



/(/ + !) 



1 

I(ITi) 



2mv 



-Ml, 



where xi = ka^/Ii, X2 = ka^/ei, the prime denotes dif- 
ferentiation with respect to the argument, and j; and 
h^^^ = iji — yi are spherical Bessel and Hankel functions, 
respectively. 

We can incorporate nonlocal effects in the sphere re- 
sponse using the formalism of Sec. A. 2. This allows us to 
generalize the non-retarded study of dispersive spheres 
reported by Fujimoto and Komaki (1968). Considering a 
metallic sphere in vacuum [i.e., ei = 1, and with medium 
2 inside the sphere described by eion(9,w) and etr(w)], a 
detailed analysis leads to 



where 



_ 2; + l(;-|m|)! (2|m|-l)!! (im\+i/2)fc 



-jljx) {[X2jl{x2)]' + gi} + e2[xjl{x)]'jl{x2) 

{[X2jiix2)y + gi} - e2[xh'i+\x)yji{x2) 



Mim = Nirr,+iy^{l + m~-l){l 
+ Ni 



m+m), 



and C'm ^ are Gegenbauer polynomials (Abramowitz and 
Stegun, 1972). Similar formulas are encountered to de- 
scribe the Coulomb excitation of atomic nuclei (Winther 
and Alder, 1979). 

The scattering matrix elements of Eqs. (33) and (58) 
depend on the composition and structure of the spheres. 
For the homogeneous particle of Fig. 13, Mie theory leads 
to (Jackson, 1999; Mie, 1908) 



.M 



-jl{Xl)x2j'i{x2) + Xij'i{Xi)jl{X2) 



hl+\xi)x2fi{x2) - Xi[h\-^>{Xi)]'jl{X2) 



(CI) 



where 



^ + i)jK^2)ji(go«) _i 



X = ka, X2 = kuy/e^, and go is given by the condition 
£\on{qo,<^) = [see Eq. (A6)]. The matrix elements tf 
are still given by Eq. (CI) with ctr [Eq. (A7)] substituted 

for 62- 



54 



APPENDIX D: List of selected symbols and acronyms 



A 

a 
B 

b 

c 

/~iE /-iM 

e 
E 
F 
f 

G 

go 

H 
h 

n 

j 

ji 

k = Lu/c 
kp 

Km 

L 

I, m 

^max 
Pi, P/ 

q 

qii = {(ix,(iv) 
q_L 

R 

R = (x, y) 
R 

r 

r= (R,^) 

ro = (Ro.^o) 
t 



t 



Vf 

W 



vector potential 
sphere or cylinder radius 
magnetic induction 
electron impact parameter 

speed of light 

electron-sphere coupling coefficients 
minus the electron charge 
electric field 
electromagnetic force 
far field amplitude 
electromagnetic Green tensor 
dipole-dipole interaction in vacuum 
magnetic field 
boundary current 
Planck's constant 

spherical Hankel function of order / 
modified Bessel function of order m 
electron current density 
spherical Bessel function of order I 
free-space light wavevector 
Boltzmann's constant 
Fermi wavevector 

modified Bessel function of order m 

length of the electron trajectory 

orbital and azimuthal numbers 

maximum orbital angular-momentum number 

electron mass 

initial and final electron wavevectors 
wavevector transferred from the electron 
wavevector parallel to a planar surface 
wavevector transfer perpendicular to v 
radial distance in cylindrical coordinates 
in-plane or _L z component of r 
cylindrical radial unit vector 
radial distance in spherical coordinates 
position coordinate vector 
electron trajectory 

electron impact parameter at time t = 
time 

sphere scattering matrices 

electron velocity vector 

electron velocity 

Fermi velocity 

scalar screened interaction 

Cartesian unit vectors 



a 
P 

7 

7e 



S 

V 
A 

p 

Pa 
a 

<t> 
V> 

V 

Vin 

<^out 

X 

OJ 
U)p 



particle polarizability 

hydrodynamic wavevector 

probability per unit frequency w 

7^ in vacuum (e = 1) 

Lorentz's contraction factor in a medium 

of permittivity e 

dielectric function 

Drude, Lindhard, and Mermin e 

longitudinal and transversal e 

Tie is an electron energy 

bulk plasmon damping rate 

free-space light wavelength 

coefiicients for TR and SPP emission 

electron charge density 

photonic local density of states (LDOS) 

boundary charge 

scalar potential 

azimuthal unit vector 

azimuthal unit vector in q space, -L qy , z 

azimuthal angle 

electron-beam divergence half-angle 

electron-spectrometer aperture half-angle 

electric susceptibility 

electron wave function 

solid angle 

frequency 

classical bulk plasma frequency 
Ritchie's surface plasma frequency 
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ir> 9r> '^r> 


Ullc-, IWL)-, llllcc-tlllllcllblUllcll 


RF,M 


UkJ LlliLldi y clcllltJlit lllc tliUU. 




CcllllUQUlUlillllcbCcliCfc; 




v^llcicllKUV IcU-lldllUll 


TATA A 


discrete- dip ole approximation 


Ur iVi 


dark-field microscopy 


T7T7PQ 


electron energy-gain spectroscopy 




electron energy-loss spectroscopy 


rjf irjiVi 


energy-filtered TEM 


e-h 


electron-hole 




electron-induced radiation emission 


r VV xliVl 


full width, at half maximum 


1 1 1 \ 1 . 1 . 1 -D 


nign- resolution rL/H/ijo 




joiiii cieiioiiy oi oiaieb 


l\ l\ 


Kramers-Kronig 


r\ r\ I V 


Korringa-Kohn-Rostoker 






T PPA/T 


low-enei^y eiecuion iiiicioocopy 


T VT) 


localized plasmon resonance 


NTT? 




IN OvJlVi 


near-field scanning optical microscopy 








Y■ofloJ^^-i^^n TTTTT Q 

renection HjUjijo 


oH/ 


secondary electron 




secondary electron emission 


oH/iVl 


scanning electron microscopy /microscope 




Gil VTQ r'O—on ri Q Ti r'oH R qtti q ti ar'Q'H'OT'i n or 
oUl id'^c-clilidliCcU. JTLdlildli oljdttcillig 




tl Lli idjL-L- JJidoiiiUii 


kjir r 


bui idce-pidbiiioli puidiiuoii 




Smith- Pur cell 


STEM 


Qr'fi TTmnc i'T'sn qtti i qqi c\Ti o or*!" tt^ti th i r'VOQrr^TiiT' 

O*^ diili-l-ii ti LjI CljiiolilloolUli tyityl^ Li Uli lllll^l \Jo\j\J\jy 


O i M 


scanning tunneling microscopy 


O ± O 


scanning tunneling spectroscopy 


TF,M 


■f T" Q n GTTl 1 GGl i~\n ol Or''i" T'r~»T1 




TTnr'T*r^ar'r^T»i7' / Tmr'Tr^arTM^o 
iiiUjiUoou^y / iiiUjUjotjU^t; 


TR 


transition radiation 


UV 


ultraviolet 


XRM 


x-ray microscopy 


ZLP 


zero-loss peak 
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